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FOREWARD 


This  report  was  prepared  by  the  Grumman  Aircraft 
Engineering  Corporation,  Bethpage,  New  York,  for  the 
Bureau  of  Naval  Weapons,  Washington,  D.C.,  under  contract 
NOw  60-0449-c,  dated  8  June  1960. 

The  work  was  performed  under  the  cognizance  of 
Grumman  Advanced  Development  by  the  Dynamic  Analysis 
and  Aerodynamics  Sections  of  the  Grumman  Engineering 
Department.  The  project  was  monitored  initially  by 
Mr.  E.  J.  Griffin  and  subsequently  by  Mr.  J.  H.  Walls, 
RAAD-223,  Loads  and  Dynamics  Branch,  Airframe  Design 
Division,  Bureau  of  Naval  Weapons. 

This  report  is  the  final  report  on  the  project,  and 
supersedes  previous  progress  reports.  It  is  presented  in 
three  complementary  volumes  as  follows: 

Volume  I  -  Theoretical  Development  (Unclassified) 

Volume  II  -  Calculation  Procedure  (Confidential) 

Volume  III  -  Results  (Confidential) 

All  three  volumes  have  the  same  basic  title.  An  Engineer¬ 
ing  Evaluation  of  Airplane  Gust  Load  Analysis  Methods. 

The  theoretical  development  and  calculation  procedures  have 
general  application  to  subsonic  aircraft.  The  illustrative 
examples,  data,  and  specific  results  are  for  the  A-6A  air¬ 
plane  . 
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SUMMARY 


Airplane  longitudinal  motion  is  described  by  a 
linear,  multi-degree-of-freedom  system  of  forced  re¬ 
sponse  equations  in  terms  of  normal  modes,  where  the 
basic  inputs  to  the  system  are  the  generalized  aero¬ 
dynamic  forces  due  to  motion  and  gust.  The  modes 
considered  in  the  analysis  are  two  symmetric  rigid- 
body  modes  (vertical  translation  and  pitch)  and  an 
arbitrary  number  of  symmetric  vibration  modes  of  the 
wing  and  fuselage  with  the  stabilizer  assumed  rigid. 

In  formulating  the  generalized  forces,  strip  theory, 
appropriately  modified  to  yield  experimental  complete- 
airplane  steady-state  stability  derivatives,  is  used. 
The  lifting  elements  (strips)  of  the  airplane  are  an 
arbitrary  number  of  wing  and  tail  panels  per  side. 
Aerodynamic  loading  of  the  fuselage  is  assumed  to  be 
acting  on  the  wing  sections  obtained  by  extending  the 
wing  to  the  centerline  of  the  airplane. 

Development  of  the  generalized  forraNi^riations 
is  carried  out  through  use  of  Dtihamel'sllJ^eghal^ 
''With  the  unsteady  aerodynamics  due  to  mtttiVi,  gusfc^^ 
and  downwash  being  described  by  indicia jL  11^ 
functions  which  include  both  compressibi\ity^and 
aspect  ratio  effects ;  The  unsteady  aerodynamics  of 
swept  planforms  is  considered  by  developing  indie i a 1 
generalized  force  functions  which  account  for  tm^ 
sweep  of  the  leading  edge  of  the  wing  and  tail  by  a 
transport  time  delay  to  the  leading  edge  of  each 
panel.  Auxiliary  differential  equations  derived  from 
exponential  and/or  damped  trigonometric  approximations 
to  the  indicial  functions  are  used  to  represent  the 
aerodynamics  due  to  arbitrary  motion  and  gusts  in 
either  the  time  or  the  frequency  domain. 

Control  system  dynamics,  including  stick  motion, 
stabilizer  rotation,  and  autopilot  control,  is  also 
provided. 

Incremental  structural  loads  (shear,  bending 
moment,  and  torque)  are  determined  for  several  wing 
spanwise  stations  and  one  tail  station.  The  wing 
loads  are  determined  using  the  mode  acceleration 
method,  with  the  aerodynamic  forces  being  approxi¬ 
mated  by  the  rigid-body  airload  distribution  which 
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would  produce  the  previously  calculated  rigid-body 
vertical  acceleration.  The  tail  incremental  shear 
is  based  on  the  force  summation  method,  with  the 
aerodynamic  forces  being  summed  over  the  panels. 

The  associated  tail  torque  and  bending  moment 
due  to  the  aerodynamic  forces  are  obtained  by  multi¬ 
plying  the  lift  by  appropriate  moment  arms. 

The  response  in  each  of  the  variables  considered 
may  be  obtained  either  in  the  form  of  time  histories, 
frequency  responses,  or  in  the  statistical  forms  of 
power  spectra  and/or  response  peaks  per  unit  time. 
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INTRODUCTION 


Early  airplane  gust  response  investigations  were 
based  on  relatively  simple  mathematical  models  of 
the  disturbance  and  the  disturbed  vehicle.  The  gust 
models  consisted  of  simplified  representative  shapes, 
such  as  the  sharp-edged  gust  and  the  ramp- type  gust. 

The  resultant  airplane  motion  was  computed  primarily 
on  the  basis  of  the  vertical  response,  with  correction 
factors  added  to  account  for  pitching  motion,  flexi¬ 
bility,  and  unsteady  aerodynamics. 

With  advances  in  the  state  of  the  art  of  aircraft 
design,  gust  response  analysis  methods  were  of  course 
gradually  improved.  And,  in  the  late  1950's,  with  the 
introduction  of  power- spectral  methods  in  the  gust  re¬ 
sponse  field  and  the  subsequent  measurement  of  actual 
atmospheric  turbulence  characteristics,  it  became 
possible  to  base  gust  response  analyses  on  a  more  real¬ 
istic  description  of  the  gust  input. 

To  complement  this  advance,  it  was  also  desir¬ 
able  to  obtain  additional  information  on  analytical 
techniques  required  for  a  realistic  description  of  the 
airplane  response.  To  provide  such  information  for 
subsonic  attack-  and  fighter-type  aircraft,  the  project 
on  which  this  report  is  based  was  established. 

A  set  of  equations  of  motion  was  formulated 
which  includes  control  system  dynamics,  a  relatively 
large  number  of  structural  modes,  and  a  fairly  com¬ 
prehensive  treatment  of  the  unsteady  aerodynamics. 
Thereafter,  successive  simplifications  in  these  equa¬ 
tions  were  introduced  to  obtain  information  on  the 
degree  of  analytical  sophistication  required  to  obtain 
gust  response  data  of  reasonable  engineering  accuracy. 
Finally,  a  comparison  was  made  between  results  from  the 
power- spectral  approach  and  those  from  the  one-minus- 
cosine  discrete-gust  method. 

This  volume,  Volume  I,  presents  the  development 
of  the  equations  of  motion  used  in  the  study.  The 
translation  of  these  equations  into  working  digital 
computer  programs  is  presented  in  Volume  II.  Finally, 
in  Volume  III,  the  numerical  results  of  the  project 
are  presented  and  discussed. 
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DISCUSSION  AND  SUMMARY  OF  THEORETICAL  DEVELOPMENT 


I .  Introduction 

The  detailed  theoretical  development  of  the  equations 
for  the  solution  of  the  flexible-airplane  longitudinal 
gust  response  problem  is  carried  out  in  Appendices  A  to  F. 
The  development  is  so  arranged  that  the  reader  interested 
in  acquiring  a  complete  understanding  of  the  details  in¬ 
volved  in  the  development  may  continue  from  Appendix  A 
right  through  Appendix  F.  Symbols  are  listed  in  Appendix  G. 

Appendix  A  presents  the  development  of  the  autopilot 
and  control  system  equations  which  consist  of  the  auto¬ 
pilot  signal,  the  stick  and  sprashpot  signals,  and  the 
stabilizer  rotation.  In  Appendix  B,  the  concept  of  the 
indicial  lift  function  is  discussed,  and  the  method  by 
which  these  functions  are  evaluated  including  compressi¬ 
bility  and  aspect-ratio  effects  is  described.  The  indicial 
lift  functions  presented  herein  account  for  the  unsteady 
aerodynamics  of  subsonic  rectangular  planforms  due  to  a 
stationary  step  gust,  a  step  change  in  angle  of  attack 
caused  by  motion,  and  the  downwash  at  the  tail  caused  by 
wing  motion  and  gust  on  the  wing.  In  Appendix  C,  the  un¬ 
steady  aerodynamics  of  swept  planforms  is  considered  by 
developing  indicial  generalized  force  functions  which 
account  for  the  sweep  of  the  leading  edge  of  the  wing  and 
tail.  The  development  of  these  functions  is  dependent 
upon  the  indicial  lift  functions  due  to  gust  and  downwash 
developed  in  Appendix  B.  Appendices  D  and  E  develop  the 
generalized  forces  which  are  the  inputs  to  the  equations 
of  motion.  These  generalized  forces  are  divided  into  the 
following  six  groupings:  (1)  motion  of  the  wing,  (2) 
motion  of  the  tail  and  rotation  of  the  stabilizer,  (3) 
downwash  at  the  tail  due  to  motion  of  the  wing,  (4)  down- 
wash  at  the  tail  due  to  gust  on  the  wing,  (5)  gust  on  the 
tail,  and  (6)  gust  on  the  wing.  Finally  in  Appendix  F, 
expressions  for  the  airplane  structural  loads  are  devel¬ 
oped.  These  expressions  are  specifically  developed  for 
several  wing  spanwise  stations  and  for  one  tail  station 
defined  by  the  juncture  of  the  fuselage  side  and  the 
stabilizer  torque  tube. 

The  final  form  of  the  equations  summarized  in  the 
subsequent  sections  of  the  main  part  of  the  report  are 
reproduced  from  these  appendices  after  a  certain  amount 
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of  algebraic  manipulation.  Also,  the  time-history  differ¬ 
ential-equations  representation  of  the  variables  is  con¬ 
verted  to  a  frequency- dependent  algebraic-equations 
representation  through  the  assumption  of  a  sinusoidal  form 
of  the  variables  as  given  by  the  following  expression: 


£(3-Std) 


fa(^)e 


i(xs  +  fa(a))-KSTD) 


(1) 


The  function  f  (ai)  represents  the  steady-state  harmonic 

a 

amplitude  of  the  assumed  sinusoidal  variable  on  the  left- 
hand  side  of  equation  (1)  .  The  constant  f  a(jo)  is  the  phase 
angle  of  the  variable  f(s)  with  respect  to  a  reference 
forcing  function  and  k.  is  the  reduced  frequency  in  radians 
per  wing  reference  semichord.  The  time  delay  s,^  is  in¬ 
cluded  for  those  special  cases  where  a  physical  time  delay 
is  associated  with  a  particular  function,  such  as  the  gust 
on  the  tail.  Equation  (1)  may  be  rearranged  as  follows: 


f  (S“Sipjj) 


-iKS^  iKS 

f  (cd)  [  e  TD][e  ] 


(2) 


where 

f(cD)  =  f  (aa)eif&(a^ 
a 


-  fa(o>)[cos  (fa(oo))+  i  sin(fa('i))>]  V 


(3) 


and 


K  *  aob^/V  =  CJDt/s 


J 


Now  let  the  Greek  letter  rho,  p,  be  used  as  a  subscript 
to  represent  the  real  part  and  the  Greek  letter  iota,  i, 
be  used  to  represent  the  imaginary  part.  Equation  (3) 
is  now  written  as 

f(^)  -  fp(<u)  +  if  t  (u>)  (4) 
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where 


fpfa)  -  fa(u)  cos<f0(o>)) 
ft(Q))  "  fa^  sin(fa(oi)) 

The  magnitude  of  the  function  is  given  by 

f,(®)  -  t£p(ffi)  +  (6) 

and  the  phase  relationship  is 

fa(a>)  »•  tan_1[ft(co)/fp(a>)]  (7) 

Equation  (2)  may  be  rearranged  further  as 

f (s~sTjj)  «  f  (d>)  [cos  KS^  -  i  sin  KSTD]eiKS  (8) 

Equations  (4)  and  (8)  are  the  primary  equations  which  will 

be  used  in  converting  the  time  history  expressions  to  sim¬ 
ple  harmonic  motion.  Furthermore,  the  required  first  and 
second  derivatives  may  be  evaluated  from  equation  (8)  as 
follows : 

f  (s-sTD>  -  ixf  (cd) [ cos  KS-f1p  -  i  sin  KSTD]eilcs  (9) 

f(s”sTD)  *  -K  f(o))[cos  KS.jp  -  i  sin  KS^le1  (10) 


where  a  circle  over  the  variable  indicates  differentiation 
with  respect  to  the  non-dimensional  time,s. 


II. 


Autopilot  and  Control  System  Equations 


The  A-6A  (A2F-1)  automatic  flight  control  system  has 
three  basic  modes  of  operation:  manual  with  stability 
augmentation,  pilot  relief,  and  command.  In  the  present 
study,  which  is  restricted  to  longitudinal  perturbations 
at  constant  velocity,  only  the  manual  mode  with  stability 
augmentation  and  the  pilot  relief  mode  in  attitude  hold 
will  be  considered.  For  these  cases,  autopilot  signals 
are  integrated  into  the  longitudinal  control  system  so  as 
to  effect  stabilizer  motion  independently  of  stick  motion. 
The  autopilot  signal  equation  given  herein  is  written  in 
general  form  to  accommodate  either  one  of  the  above  cases. 
Signals  from  the  motion  sensing  devices  are  written  in 
terms  of  the  normal  mode  components  of  airplane  motion. 

The  equation  describing  the  autopilot  signal  is  reproduced 
from  equation  (A5)  as 


p-(s)  +  Ll0p.(s) 


N 

-  2  [L 
n=l 


OO 

ln?n<s> 


+  L2A 


(s)  +  L3nfn(s) ] 


(ID 


Simple  harmonic  motion  substitution  yields  for  the  auto¬ 
pilot  equation 


(12) 


where 


P-(“)  *  M-p(o))  +  int(o>) 


(13) 


and 


Npn  *  <L10  '  iK) [ (L3n  ‘  K  Lln^  +  ixL2n] 


2  2 

Dp.  -Li0  +  K 


(14) 


In  the  stick-free  condition,  the  stick  deflection  re¬ 
sulting  from  aircraft  motion  is  determined  by  an  artificial 
feel  system  consisting  of  a  feel  spring,  a  sprashpot,  and 
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two  bobweights.  The  feel  spring  and  the  stick  inertia  are 
represented  by  the  left  side  of  the  stick  equation,  while 
the  reaction  force  from  the  sprashpot  and  the  inertial 
forces  from  the  bobweights  (in  terms  of  normal  mode  com¬ 
ponents)  are  contained  on  the  right  side.  The  sprashpot 
equation  represents  a  spring  and  dashpot  in  series  connected 
to  the  stick  below  the  stick  pivot  point.  Nonlinearities 
in  the  control  system  due  to  friction  have  been  neglected. 

The  stick  and  sprashpot  signal  equations  are  given  by 
equations  (A17)  and  (A18)  which  are  reproduced  below. 

OO  N  O  O  o 

T(s)  +  H30t(s)  -  2  H3n4n(s)  +  H3T1n(s)  (15) 

q(s)  +  H40ri(s)  -  H40t(s)  (16) 

Simple  harmonic  motion  input  into  these  expressions  re¬ 
duces  them  to  the  following  relationships: 

t(o))  -  [-k2Z  iH3n^n^  +  J/fH30  “  K  ^  (17> 


T)(<o)  “  [H40(H40  -  1k)/(H40  +  k2)]t(gd) 


(18) 


where 


t(co)  -  Tp(o))  +  ixt(a>) 

n(«)  -  *lp(®)  +  i\(“) 


(19) 


The  stabilizer  deflection  results  from  a  superposition 
of  input  signals  from  the  autopilot  and  the  stick.  The 
equation  describing  the  stabilizer  rotation  is  reproduced 
below  from  equation  (A26) . 

*(s)  +  H2Q7(s)  +  H217(s)  -  H21^(s)  +  H22t(s)  (20) 

Substitution  of  simple  harmonic  motion  into  the  above 
equation  yields, 


11 


7(o>)  -  [N7lp.(a>)  +  N72t(o>)  3/D7 

where 

7<«>)  “  7p(<u)  +  i7t(“) 

and 

■71  "  C (H2i  "  k2)  -  i*H20]H21 

N72  "  1 <H21  ‘  k2)  “  1kH201H22 
D7  "  <H21  "  k2>2  +  A  20 


(21) 


(22) 


(23) 


J 
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III.  Indicial  Lift  and  Generalized  Force  Functions 

Indlclal  lift  functions  describe  the  unsteady  lift 
build  up  on  rectangular  planforms  due  to  planform  motion, 
gust  on  the  planform,  and  downwash  effects  on  an  aft 
planform  due  to  either  gust  or  motion  of  the  forward 
planform.  In  the  case  of  a  swept  planform,  a  gust  or 
downwash  envelops  successive  s panwise  portions  of  the 
planform  at  different  times,  and  therefore,  the  lift 
variation  on  the  planform  is  changed  by  the  sweep  of  the 
planform  leading  edge.  The  effect  of  sweep  is  introduced 
through  the  development  of  indicial  generalized  force 
functions  for  eadh  mode  (rigid-body  and  flexible  modes) 
which  depend  upon  the  indicial  lift  f  met  ions,  planform 
sweep,  mode  shapes,  and  aerodynamics . 

A  discussion  emphasizing  some  of  the  more  significant 
aspects  of  the  analysis  for  the  indicial  lift  and  general¬ 
ized  force  functions  follows. 

The  present  analysis  is  based  on  indicial  lift  func¬ 
tions  for  the  lift  build  up  on  rectangular  planforms  due 
to  a  stationary  step  gust  and  a  step  change  in  angle  of 
attack  due  to  motion.  The  lift  build  up  on  a  rectangular 
planform  due  to  a  step  change  in  angle  of  attack  as  a  re¬ 
sult  of  a  vertical  gust  with  zero  gust  front  velocity, 
normalized  with  respect  to  the  steady-state  lift  and  in¬ 
cluding  compressibility  and  aspect  ratio  effects,  is 
referred  to  as  the  Kffssner- type  f motion.  This  function 
is  approximated  by  an  analytical  expression  made  up  of 
the  sum  of  exponential  terms  as  given  by  either  of  equa¬ 
tions  (B7)  or  (BIO)  . 

He)  -  *(0)  +  L  ak(l  -  e  °k*)  (24) 

k-1  * 


where  the  letters  W  and  T  in  equations  (B7)  and  (BIO) 
representing  gust  on  the  wing  and  gust  on  the  tail,  re* 
spectively,  are  deleted  here  for  convenience.  The  maxi¬ 
mum  index,  K^,  designates  the  number  of  exponentials  used 

in  approximating  the  indicial  lift  function  due  to  gust. 
The  purpose  of  the  constant,  if (0),  is  to  permit  the  re¬ 
duction  of  the  equation  to  the  limiting  case  where  the 
lift  build  up  due  to  gust  is  assumed  Instantaneous,  that 
is,  where  no  transient  effects  are  present.  The  reduction 


13 


of  the  equation  to  this  Instantaneous  lift  build  up  may 
be  accomplished  by  letting  0(0)  -  1.0  and  a^  -  0  for 

k  -  1,  2,  K^.  The  constants  a^  and  will  be  a 

function  of  the  Mach  number  and  aspect  ratio  (see  Figure 
B2)  . 


The  lift  build  up  on  a  rectangular  planform  due  to 
a  step  change  in  angle  of  attack  as  a  result  of  motion, 
normalized  with  respect  to  the  steady-state  lift  and 
including  compressibility  and  aspect  ratio  effects,  is 
referred  to  as  the  Warner- type  function.  This  function 
is  also  approximated  by  an  analytical  expression  consist¬ 
ing  of  a  constant  plus  the  sum  of  a  number  of  exponential 
terms.  (See  equations  (B15)  and  (B17)).  Once  again  the 
expression  is  reproduced  below  without  the  identifying 
letters  W  or  T. 

-as 

0<s)  -  0(0)  +  S  b  (1  -  e  k  )  (25) 

k-1  k 


The  constant  term  arises  from  the  instantaneous  lift  at 
time  zero,  which  for  incompressible  flow  and  infinite 
aspect  ratio,  is  equal  to  one-half.  Hie  other  coeffi¬ 
cients  and  exponents  appearing  in  the  Wagner-type  func¬ 
tion  are  evaluated  in  a  manner  similar  to  the  constants 
in  the  Kdssner-type  function.  A  typical  plot  of  a 
Wagner-type  function,  and  a  comparison  with  a  Ktfssner- 
type  function  is  shown  in  Figure  B4. 

The  normalized  unsteady  downwash  experienced  at  the 
tall  due  to  step  gust  on  the  wing  or  a  step  change  in 
angle  of  attack  due  to  wing  motion  is  referred  to  as 
an  indlclal  downwash  function  (see  Figure  B5) .  Since 
the  downwash  has  the  characteristics  of  a  gust  distur¬ 
bance,  the  lift  build  up  on  the  tail  due  to  downwash  is 
assumed  to  be  given  by  Duhamel's  integral  in  terms  of 
the  indicial  downwash  function, 

w(Mg,  M,  ABg ,  8*,  8*. 

and  the  tail  KUssner-type  function.  (See  equations  (B20) 
and  (B21)) . 

«(•)  -  !l  w(Mg,  M,  ARy,  8*.  s|,  a)^(8-«)<i«  (26) 
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where,  for  downwash  due  to  motion,  the  gust  front  Mach  num¬ 
ber  is  taken  as  Mtt  -  °°,  and  for  downwash  due  to  a  station- 

ary  gust,  Mg  -  0.  The  parameters  s£  and  s^.  in  the  indicial 

downwash  function  are  defined  in  Figure  B1  ,  and  ARy  is 

the  aspect  ratio  of  the  wing.  Hie  time  history  variation 
given  by  equation  (26)  is  indicated  in  Figure  B7. 


The  indicial  generalized  force  functions  developed 
herein  are  defined  in  terms  of  a  strip  theory  representa¬ 
tion  of  the  unsteady  aerodynamics.  The  indicial  general¬ 
ized  force  functions  describe  the  normalized  time  varia¬ 
tion  of  the  generalized  forces  on  a  swept  planform, 
asstiming  a  lift  build  up  variation  on  each  strip  given 
by  the  Ktfssner-type  function,  as  a  step  gust  envelops 
succeeding  spanwise  portions  of  the  planrorm.  In  the 
present  analysis,  which  is  for  a  stationary  gust,  it  is 
assumed  that  the  spanwise  rate  of  gust  envelopment  is  a 
function  of  wing  and  tail  leading-edge  sweep.  Similarly, 
the  rate  of  downwash  envelopment  of  the  tail  is  assumed 
to  be  a  function  of  the  tail  leading-edge  sweep. 

To  develop  the  indicial  generalized  force  functions, 
the  wing  and  tail  are  divided  into  a  number  of  panels, 
where  the  lift  and  pitching  moment  about  the  reference 
axis  on  each  panel,  determined  by  the  local  geometry  and 
aerodynamic  properties,  are  assumed  to  start  at  the 
time  the  gust  (or  downwash)  reaches  the  panel  spanwise 
midpoint  at  the  leading  edge.  Assuming  a  KUssner-type 
variation  of  the  lift  build  up,  the  total  generalized 
force  after  the  gust  or  downwash  has  progressed  to  a 
given  outboard  station  is  determined  by  the  sum  of  the 
generalized  forces  over  all  the  inboard  panels. 

Analytical  approximations  to  the  indicial  generalized 
force  functions  are  obtained  by  considering  each  function 
to  consist  of  a  number  of  damped  trigonometric  and/e* 
pure  exponential  terms.  For  gust  on  the  wing,  the  expres¬ 
sion  for  the  indicial  generalized  force  function  is 
given  by  equation  (C13)  which  is  reproduced  below. 

2  -JJ  s 

•  *w(0>  +  2,  «it<1  -  «  ^  > 


+  [«U3tl  -  •  °w* 

-  *&*  ,ln  <ctW3*)  ] 


(27) 
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In  the  case  of  gust  on  the  tall,  the  indlcial  generalized 
force  function  is  similar  to  equation  (27)  with  the 
letter  W  being  replaced  by  T.  Furthermore,  the  downwash 
effects  are  represented  by  the  same  type  of  expression, 
where  now  the  left  hand  side  of  equation  (27)  is  replaced 

by  6  (s)  for  downwash  due  to  motion,  and  by  6u(s)  for 

downwash  due  to  gust;  also,  the  symbol  W  is  replaced  by 
M  and  G,  respectively.  It  is  noted  that  the  assumption 
of  rigid  tail  in  the  present  study  renders  the  tail 
indicial  generalized  force  functions  for  all  modes  nearly 
identical  numerically  to  those  for  the  first  mode,  that 
is,  according  to  equations  (Cl5),(C17),  and  (C19), 

♦*(.)  - 

«£<•>  -  «i<*> 

6“(s)  Sr  6«(S) 

TPM 

Therefore,  the  functions  tf/£(s),  6^(s),  S^s)  are  used  for 
the  other  modes  as  well. 

For  motion  of  the  wing  and  tail,  the  lift  build  up 
occurs  simultaneously  over  the  entire  planform.  Therefore, 
sweep  effects  do  not  modify  the  wing  and  tail  Wagner-type 
functions,  and  equation  (25)  is  also  a  valid  approxima* 
tion  to  the  indicial  generalized  force  functions. 


V  for  all  m 


(28) 


J 
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IV.  Equations  of  Motion  and  Associated  Aerodynamic  Lag 

Equations 

In  this  analysis  it  is  assumed  that  the  airplane  en¬ 
counters  a  stationary  gust  symmetrically,  and  that  no 
spanwise  variation  in  the  gust  exists.  The  dynamic  response 
of  the  airplane  due  to  the  gust  is  obtained  by  a  modal 
analysis  in  which  the  airframe  is  allowed  six  degrees  of 
freedom:  two  symmetric  rigid-body  modes  (vertical  trans¬ 
lation  and  pitch) ,  and  four  symmetric  vibration  modes  of 
the  wing  and  fuselage  with  the  stabilizer  assumed  rigid. 

Note  that  the  flexible  modes  considered  here  are  coupled 
bending- tors  ion  modes,  and  that  the  word  "bending"  in  a 
designation  such  as  "first  symmetric  wing  bending  mode" 
is  intended  only  as  an  indication  of  the  dominant  feature 
of  that  mode. 

Airplane  longitudinal  motion  is  described  by  a  linear, 
multi-degree-of-freedom  system  of  forced  response  equa¬ 
tions  in  terms  of  the  aforementioned  modes,  where  the 
basic  inputs  to  the  system  are  the  generalized  forces 
due  to  motion  and  gust.  Structural  damping  is  conserva¬ 
tively  neglected  throughout  the  analysis.  In  formulating 
the  generalized  forces,  strip  theory,  appropriately  modi¬ 
fied  to  yield  correct  complete-airplane  steady-state 
stability  derivatives,  is  used.  The  lifting  elements 
(strips)  of  the  airplane  are  ten  wing  panels  and  six  tail 
panels  per  side.  Aerodynamic  loading  of  the  fuselage  is 
assumed  to  be  acting  on  the  wing  sections  obtained  by 
extending  the  wing  to  the  centerline  of  the  airplane.  The 
differential  equations  used  in  obtaining  the  airplane  re¬ 
sponse  are  derived  from  the  following  special  form  of 
Lagrange's  equations  of  motion, 

ar  (f^y)  +  +  <£<'>  (29) 


for  m  -  1,  2,  . .  .N 

where  T(t)  and  V(t)  are  the  kinetic  and  potential  energies, 
and  Q^(t)  and  Q^(t)  ate  the  generalized  forces  in  the  mfc^ 
mode  due  to  motion  and  gust,  respectively. 

In  converting  equation  (29)  to  non-dimensional  time, 
the  following  substitutions  are  made: 
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(30) 


d 

Ht  = 


(31) 


It  is  noted,  for  example,  that  the  generalized  force  due 
to  motion,  Q^(t) ,  after  conversion  should  be  expressed 

but  for  brevity  and  to  reduce  excessive  symbols  for 
new  functional  notations,  it  is  designated  simply  as 
<£(s).  After  conversion,  V(t),  £m(t),  and  Q^(t)  are  simply 
expressed  as  V(s),  £m(s)  and  Q^j(s)  .  The  kinetic  energy 
terms,  after  substitution  of  equation  (31)*  become 


d  (  ST(s) _ \  _  d  /dT(s)  \ 

( V  \  (v_)  a|m(S)  J  ds  [  a|m(s)  j 


(32) 


where  the  circle  over  the  variable  indicates  differen¬ 
tiation  with  respect  to  the  non-dimensional  time,  s.  Equa¬ 
tion  (29)  becomes 


fdSfeHSw  <33> 


The  kinetic  energy,; is  expressed  in  terms  of  the  normal 
velocity  coordinates,  £  (s),  the  stabilizer  rotation 

O  ^ 

velocity,  7(s),  and  the  generalized  mass  terms  as  follows: 


T(s)-(V/bR)2j  jui^£(.)  +  2Mry;(8)ln(s)l 


+  M^72  (8) 


/2g 


(34) 


Note  that  the  factor  "g"  is  introduced  because  the  general¬ 
ized  mass,  as  used  herein,  is  defined  in  terms  of  pounds  mass 
rather  than  slugs. 
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The  potential  energy  is  expressed  in  terms  of  the 
normal  coordinates,  £  (s),  and  the  generalized  spring 
terms  as  follows*  m 

V(s)  =  %  2  (»)  I/g  (35) 

m=l 


<.u 

Hie  generalized  force  in  the  nrn  mode  due  to  motion, 

M 

Q^(s),  is  composed  of  the  four  distinct  inputs  associated 

with  wing  motion,  tail  motionj  stabilizer  rotation,  and 
downwash  on  the  tail  due  to  wing  motion.  The  wing  and 
tail  generalized  forces  due  to  wing  motion  and  tail  motion 
are  made  up  of  non-circulatory  contributions  (apparent 
mass  effects)  described  in  terms  of  acceleration  and 
velocity,  and  circulatory  contributions  given  in  terms  of 
velocity  and  displacement.  The  circulatory  contribution 
is  included  in  the  analysis  through  the  use  of  Duhamel's 
integral  in  conjunction  with  the  Wagner-type  function 
given  by  equation  (25) .  Use  of  Laplace  transform  tech¬ 
niques  results  in  an  expression  for  the  generalized  force 
due  to  motion  which  is  made  up  of  a  constant  term,  due  to 
0(0)  in  equation  (25)  and  a  set  of  first  order  differential 
equations  due  to  its  exponential  terms.  The  downwash  on 
the  tail  due  to  wing  motion  is  assumed  to  result  from  the 
rigid-body  airplane  response  only,  the  reasoning  being 
that  downwash  due  to  flexible  wing  motion,  which  has 
relatively  small  amplitude  and  high  frequency,  decays  be¬ 
fore  it  reaches  the  tail.  Hie  downwash  is  assumed  to  be 
constant  across  the  tail  span  and  to  be  defined  in  terms 
of  the  wing  motion  at  a  s panwise  station  corresponding  to 
the  tail  mean  aerodynamic  chord.  In  the  equations  pre¬ 
sented  herein,  the  indicial  lift  functions  due  to  tail 
motion  and  stabilizer  rotation,  have  been  assumed  equal. 
Thus,  these  two  motions  are  combined  in  a  single  set  of 
equations . 

The  expression  for  the  total  generalized  force  due 
to  wing  motion,  Q^m(s),  tail  motion  and  stabilizer  rotation, 

Q2m(8)»  and  downwash  on  the  tail  due  to  wing  motion,  Q3m(®)> 

is  given  by  the  sum  of  the  three  component  generalized 
forces : 


<£<■> 


3 

"r-lQrxn(S)  *  1  -  m  -  N 


(36) 
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The  equations  of  the  generalized  forces  for  the  three 
types  of  inputs  are  summarized  below. 

For  wing  motion,  (r  =  1),  the  generalized  force  con¬ 
sists  of  a  non-circulatory  contribution,  $>^(s),  plus  a 

wc  ® 

circulatory  contribution,  $  (s),  multiplied  by  a  constant, 

-Cy,  which  is  factored  out  of  the  final  equations. 

Qlm(s)  -  -Cv[tf(S)  +o”C(s)]  (37) 


for  m  *  1,  2,  . .  .N 


2  2 

where  Cy  -  V  /(bpg) .  The  expression  for  the  non-circulatory 
contribution  to  generalized  force  is 


*f(s) 


N 

2 

n=l 


[A, 


WN 

inn 


oo 

Sn(s) 


+  B 


WN 

mn 


«»(•)] 


(38) 


where  the  coefficients  in  the  above  expression  are  defined 
by  equation  (E24)  and  (E25) .  The  circulatory  contributions 
to  the  generalized  force  are  separated  into  two  parts,  the 
instantaneous  lift  term  and  the  unsteady  lift  which  is 
approximated  by  two  term  exponential  representation  of  the 
Wagner-type  function.  From  equations  (E38)  to  (E40), 

+  <39> 

where 

*>)  -  «V(0)F^(s)  (40) 


-,*♦!*<•>  <41> 


The  unsteady  lift  terms  of  equation  (41)  are  defined  by  the 
solutions  to  the  differential  equations  given  by  equation 
(E41), 
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*mk<8>  +  ^mk<8>  “  bXFm<8>’  k  "  I-  2 


(42) 


The  input  to  both  equations  (40)  and  (42)  is  defined  in 
terms  of  the  wing  motion  by  the  following  relationship: 

f"(s)  =  S  [bJ^L(s)  +C™£(s)]  (43) 

m'  _  i  mn  n  '  ran  n 
n=i 

The  coefficients  in  these  expressions  are  defined  by 
equations  (E28)  and  (E29) . 

The  expressions  for  the  generalized  forces  for  tail 
motion  and  stabilizer  rotation  (r  =  2)  are  similar  to 
those  for  wing  motion.  These  expressions  are  reproduced 
from  Appendix  E,  Section  II.  The  generalized  force  is 

<ws>  - 


for  m  =  1,  2,  3,  ...N.  The  non- circulatory  contribution 
to  the  generalized  force  is  given  by 


*™(.) 

m 


+ »™?<«> 


(45) 


In  a  manner  similar  to  the  circulatory  contributions  to 
the  wing  generalized  force,  the  expressions  for  the  tail 
are 

^(3>  -  *>>  +  "£<«>  (46) 

*L<»>  -  M°>Fm<s>  <47) 

£<•>  -£•*<•> 

T 

where  ^^(s)  are  solutions  to  the  differential  equations 


(48) 
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*Ik<s>  +  eKkf*)  -  »  -  1,  2 


and 


N 


F»  =S.(b!5L(s)  +C«LW1  +B^7(s) 


in 


n=l 


nrn^n 


tim^n 


my 1 


+  C^(S) 


(49) 


(50) 


Finally  the  generalized  force  due  to  downwash  on  the 
tail  caused  by  wing  motion  (r  =  3)  is  defined  in  terms  of 
the  downwash  of  the  wing  three  quarter  chord  point  and  an 
assumed  lift  build  up  which  includes  the  unsteady  down- 
wash  as  described  by  Hobbs  in  Reference  1  .  The  expression 
for  the  generalized  force  is  reproduced  from  equation  (E83) . 

Q3m<s>  -  -VESW  (51) 


M 

The  constant  is  defined  (equation  (E84))  in  terms  of 

the  stabilizer  geometry  and  aerodynamics .  The  generalized 
force  variation,  A"(s) ,  is  given  as  the  sum  of  the  com- 

ponent  parts,  ^"^(s),  obtained  as  the  solutions  to  a  set 

of  three  differential  equations.  The  input  to  these  equa¬ 
tions  is  the  wing  motion,  composed  of  normal  mode  velocity 
and  displacement  terras,  modified  by  the  downwash  slope, 
de/da.  Because  of  the  assumption  of  a  rigid  tail,  the 
generalized  force  variation  is  numerically  nearly  indepen¬ 
dent  of  mode  and  is  therefore  approximated  by 


AM(s) 
mv  ’ 


2  A?,,(s) 


k-0 


mk 


for  ra  =  1 


a«(s) 

m 


A«(3) 


2  £  m  £  N 


A 


}  (52) 


J 


The  right  side  of  the  first  of  equations  (52)  is  obtained 
from  the  solution  to  the  following  set  of  differential 
equations  under  the  assumption  that  the  indicial  general¬ 
ized  force  function  may  be  represented  analytically  by  a 
constant  plus  a  maximum  of  three  terms,  two  of  which  are 
pure  exponentials  and  the  other  a  damped  sinusoid. 
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These  differential  equations  are  reproduced  from  equations 
(E88)  to  (E90)  as  follows: 


A10<8>  -  6M<0>FM<s’8M> 

^(s)  +R1a”1(s)  = 

A12(s)  +R3A12(s)  =  R4FM<8“8M)‘ 
Al3(8>  +R5A13(s)  +  Vo[A13(6) 
R8FM^s“sm^»  R  ^ 


k  ■=  1 

k  =  2 

+  R7FM^“8M^ 


(53) 

i 

| 

1  (54) 


J 


Note  that  in  the  case  where  three  exponential  terms,  rather 
that  two  exponentials  and  a  damped  sinusoid,  approximate 
the  indicial  generalized  force  function  the  last  of  equa¬ 
tions  (54)  may  be  reduced  to  this  condition  by  letting 
Rfi  **  0.  The  input  to  equations  (  53)  and  (54)  is  given  by 
equation  (E91)  as 

fm(,)  +c$n<8)1  <55> 


where  the  coefficients  in  equation  (55)  are  defined  by 
equations  (E92)  and  (E93)  .  The  R^  coefficients  in  equa¬ 
tion  (54)  are  easily  evaluated  by  comparing  equations  (54) 
with  equations  (E88)  to  (E90) .  Therefore, 


4  -  aXLall 

13  “  a12 

14  "  a12a12 
L5  "  2a13 


-a? 


13 


! 

(“13) 2  +  (^)2  r8  "  a13c?3"*l3tD13  1 


f  (56) 
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The  time  delay  is  the  time  between  wing  motion  and  the 

assumed  start  of  the  resulting  tail  lift  build  up  due  to 
downwash. 

The  generalized  force  in  the  nr  mode  due  to  gust, 

(^(s)  in  equation  (33),  consists  of  the  following  three 

parts:  gust  on  the  tail,  gust  on  the  wing,  and  downwash 
on  the  tail  due  to  gust  on  the  wing.  Transport  time  delays 
required  for  the  gust  front  to  travel  (the  gust  front  is 
actually  stationary  in  space  but  does  "travel"  relative 
to  the  airplane)  from  the  assumed  zero  reference  time  to 
the  wing  and  tail  and  for  the  wing  downwash  to  reach  the 
tail  are  also  included  in  the  definition  of  the  generalized 
forces . 

The  total  generalized  force  due  to  gust  may  now  be 
defined  in  terms  of  the  following  contributions :  down- 
wash  on  the  tail  due  to  gust  on  the  wing,  i  gust 

on  the  tail,  Q5m(s) ’  and  gust  on  the  wing,  Q^m(s) . 

c£(s)  -  2  o  (s)  (57) 

r=4 


Expression  for  evaluating  the  generalized  force 
arising  from  downwash  on  the  tail  due  to  gust  on  the 
wing  (r=4)  is  given  by  equation  (D22)  and  is  reproduced 
below. 


<WS>  -  CVDm£l<s> 


(58) 


2  2  G 

where  =  (V  /b^g) ,  and  the  coefficient  depends  on 

the  tail  geometry  and  aerodynamics  and  is  evaluated  by 
equation  (D23) .  The  generalized  force  variation* 

Q 

A^(s),  is  evaluated  for  the  first  mode  only,  because 

the  assumption  of  a  rigid  tail  makes  the  unsteady  aero¬ 
dynamic  functions  due  to  gust  and  downwash  nearly  in¬ 
dependent  of  the  mode  shape  for  the  six  normal  modes 
considered  in  this  problem.  The  expression  for  evalu¬ 
ating  the  generalized  force  variation  is  given  by 
equation  (D24)  as 

3 

A?(s)  -  2  A<jWs)  (59) 

1  k-0  iK 
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The  component  parts  of  the  generalized  force  variations  as 
given  by  equations  (D25)  to  (D27)  are  reproduced  below. 


4><-> 


*G(0)x(s-sG) 


(60) 


A^Cs)  +  Q1A^1(s)  -  Q2x(s-sg),  k 


Ai2(s)  + 

A^s  (s)  +  Q5A^3(s) 
Q8X(8-sg), 


=  Q4x(s-sg),  k  -  2 


(61) 


+  Q6/o[A13(<3>  +  Q7X(e*8G)]d6- 
k  =  3 


J 


where  sG  is  the  time  delay  between  the  airplane  zero  time 

reference  point  and  the  assumed  start  of  tail  lift  build 
up  due  to  downwash  resulting  from  gust  on  the  wing.  The 
coefficients  appearing  in  equations  (61)  are  defined  as 


*2 

“3 

<>6 

<>7 

Qe 


-a 


11 

G 

11 


“llbe/Boj 


‘12 


”a12a12 

2al3 


de/£a 


(“l3) 2  +  (CD^)  2 


+a 


13 


de/daj 

"ta13a13*  ^3^13^  ]^e/^a| 


k  -  1 

k  -  2 

•  k  -  3 

> 


> 


J 
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where  the  constant  (3e/da|  is  the  absolute  value  of  an 
assumed  negative  downwash  slope  due  to  gust  on  the  wing 
(also  wing  motion) . 

The  form  of  the  expressions  describing  the  generalized 
forces  for  gust  on  the  wing  and  tail  is  similar  to  the 
expressions  for  downwash  on  the  tail  due  to  gust  on  the 
wing.  These  relationships  are  found  in  Appendix  D, 

Sections  I  and  II,  and  they  are  reproduced  below. 

The  generalized  force  due  to  gust  on  the  tail  (r  =  5) 
is  given  by 

Q5m(s)  -  Cvd£i£(s)  (63) 


T 

where  is  evaluated  using  equation  (D7) .  The  generalized 
force  variation  is 


3 

2  ¥ 

k=0 


T 

lk 


(s) 


(64) 


where  the  component  parts  of  the  generalized  force  are 
given  by 


¥ 


io<s>  =  ^t(0)x(s-st) 


(65) 


^x(s)  +  OjY^s)  =  02x(s-st),  k  -  1 


J^2(s)  +  03^2(s)  =  04x(s-st),  k  -  2 


■\ 


°T 


¥ 


i3(s)  +  05yj3(s)  +  o6/gt^3(d)  +  o7x(d-sT)]dd 


°8x(s-st) 


k  -  3 


(66) 


and  s^  is  the  time  delay  between  the  airplane  zero  time 

reference  point  and  the  start  of  tail  lift  build  up  due 
to  gust  on  the  tail.  The  coefficients  appearing  in  equa¬ 
tion  (66)  are  defined  by 
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°1  "  all 

T  T 
°2  “  allall 

°3  “  a12 

T  T 
°4  "  a12a12 

°5  “  ^a13 

°6  *  <a13)2  +  ^13)2 
°7  -  -a13 

°8  *  a13°13  “  b13CD13 


k  -  1 


k  -  2 

) 


(67) 


k  =  3 


J 


J 


Finally  the  generalized  force  due  to  gust  on  the  wing 
(r  «  6)  is 


<ws>  ■  vK<*> 


(68) 


w 

where  D  is  evaluated  using  equation  (D  7  ) .  The  general- 


m 


ized  force  variation  and  its  component  parts  are  dependent 
on  the  mode  shapes  (including  rigid  and  flexible  airplane 
modes) ,  and  are  given  by 

3 


*£<«> 


2  TW,  (s) 
k=0  ^ 


where 


^mO^  *  ^O(-)  =  %(°>X(S-SW) 


(69) 


(70) 


^1<S>  +  -Nm2X(s-sW>>  k 

^2<s)  +Nm3*S2<s>  *Nm4X(s-sw),  k 


1 

2 


(71) 


*£3<s>  +Nm5*m3(s>  +  IW>m3<8>  +  ld5 


J 
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The  coefficients  in  equation  (71)  are 


Nml  "  “ml 


N. 


m2 


W  W 
aml“ml 


N. 


m3 


a. 


,W 


'm2 


Nm4  *  am2am2 
Nm5  ’  2* 


N 


m6 


«&)2  +  (^3) 2 


A 


N 


m7 


-a 


m3J 

W 

m3 


N. 


m8 


W  W  KW  „  W 
am3am3  tnS^mS 


J 


(72) 


The  final  form  of  the  equations  of  motion,  equation 
(33),  may  now  be  obtained  by  evaluating  the  contributions 
of  the  kinetic  and  potential  energies  using  equations  (34) 
and  (35)  and  substituting  for  the  generalized  forces  due 
to  motion  and  gust  using  equations  (36)  and  (57) . 

n!1tAm.«n<*>  +  BmJn<8>  +  «««»<•>  1 

oo  ° 

+  V7  <s)  +  V7(s)  +  W(s) 


my  ‘ 

+Jt»"k(s)  +  *L(S))  +  d“  .sA(.) 


k-1 


my 1 


m 


k-0 


‘lk' 


s  -tDSAlk<8>  +  +  o"i(8)]  (73) 


k-0 


'mxlk 


m  mk 


for  m  -  1,  2,  . . . ,  N 

The  constants  in  equation  (73)  are  defined  as  follows: 
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A 

inn 

L 

mn 

mn 


Anry 

®nry 

^nry 


AI^  +  A™  +  Mn  ^ 

mn  mn  ran 

BWN  +  gTN  +  0  ( 0)BWC  +  0fp(O)BTC 
mn  "mn  ^Wv  '  ran  Tv'tnn 

V<»«  +  'V°>C£1  +  (VV>2“Xn 

A™  +  %  1 

b™  + 

V°>C-£ 


(74) 


(75) 


In  equation  (73)  the  unsteady  aerodynamics  associ¬ 
ated  with  gust  on  the  wing,  gust  on  the  tail,  downwash 
due  to  gust  on  the  wing,  and  downwash  due  to  wing 
motion  are  summed  from  k  ■  0  to  k  *  3.  Hie  component 
parts  represented  by  the  k  m  0  contribution  are  included 
in  order  to  be  able  to  introduce  the  simplification  of  in¬ 
stantaneous  aerodynamics  in  the  equations  of  motion. 

This  simplification  is  introduced  in  the  equations  by  set¬ 
ting  the  constant  terms  in  equations (53) ,  (60),  (65)  and  (70) 
equal  to  unity,  that  is,  6^(0)  *  <5G(0)  “  ^T(0)  -  ^w(0)  -  1.0. 

Associated  with  this  change,  the  unsteady  aerodynamic 
terms  as  defined  by  equations  (54),  (61),  (66)  and  (71) 
are  made  equal  to  zero  through  appropriate  use  of  equations 
(56),  (62),  (67)  and  (72).  In  the  case  where  unsteady 
aerodynamic  terms  are  included,  the  constants  6M(0),  6q(0), 

^,(0)  and  if/y( 0)  are  set  equal  to  zero. 


Before  simple  harmonic  motion  is  introduced  into 
equation  (73),  the  expressions  for  solving  the  unsteady 
circulatory  generalized  forces  due  to  motion  are  rearranged, 
that  is,  equations  (42)  and  (49)  become, 


for  wing  motion: 

&<•> + 

*»2<s>  +  F2<2<8> 


j'VliW+'.t.Wl 

n-i 

F3'*ml<8>  +  FAW1 


(76) 
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where 


Fn  -  P 


W 


ran 


'ran 


RW 

,  WAW  /,  wftw 

b2p2/bl^l 

,WRW_WC 

blPlBran 

,  WftWrWC 

VlCmn 


and  for  tail  motion: 

Sral<s>  +  Jl*ml<8> 


N 

2  t  Itnn^_  (8)  +  «J_n£_  (s)  ] 
_  i  ran  n  ran  n 

n- 1 


<■>  + 


+  1^7(8)  +  ^7(8) 
J3[^l(s)  +  V>)1 


where 

J1  " 

J2  -  P? 


J3  - 

m  rp  rp 

bjpj/bjp 

*ran” 

vT-T-TC 

blPlBmn 

J  - 

b^C™ 

mn 

llran 

I  - 

nry 

vTflT«TC 

blPlBnry 

nry 

KTRTrTC 

blPlCra7 

(77) 
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(78) 


(79) 
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Introduction  of  simple  harmonic  motion  in  equation  (73) 
as  discussed  in  Section  I  yields,  for  the  equations  of 
motion, 


Hie  expressions  for  the  aerodynamic  lag  equations  in  their 
order  of  appearance  in  equation  (80)  are  summarized  below. 


<a<®> 

■&<<“> 

•i<«> 


[  (F^IO/CfJ+k2)  ] 

[  (F2-1k)/(f|+k2)  ]  F3tF1  +  1*14^  <u>) 


N 


[  (J^/Cjf+K2)  ]  U  [J^  + 


N 


+  [Jnry  +  ^nry1  *<®> 


[  (J2-iK)/(J2+K2)  ]  +  ix]  ^(ui) 


J 


J 
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,M 


A 10 (a.)  -  6m(0)  fM(03) 

A^Cco)  -  [(R1-iK)/(R^+K2)]R2fM(oD) 
A*2(o>)  -  [(R3-iK)/(R2+K2)]R4fM(03) 

A  13(03)  =  [  ^m/^M  ^  f  M  (^) 


Where 


>  (83) 


Nm  -  (R5R8  +  R6R7)k2  -  R2R7  +  itR5R6R7  +  R8(R6-k2)1k 


and 


°M-R6 

9  9  U. 

-  (2R6  -  R3)k  +  k 

a 

0 

N-/ 

[  2 
2 

[C„  +  IkbJJh _(03)  |  [cos 

^n-1  i 

Ai0(o>)  ■ 

-  - 

de/da  j  6G(0)fG(a3) 

A11  (03) 

-  1  (Q1-1k)/(Q^+k2)  ]  Q2fG(o3) 

A^2(03) 

-  [  (Q3-1k)/(Q3+k2)]  Q4fG(oi) 

A^(03) 

-  ING/DG]fG(o)) 

J 


where 


(84) 


Ng  -  (Q5Q8  +  Q6Q7)*  “  0?*7  +  +  %<Q6 >H 


and 


D0  -  q|  -  (2Q6  -  sf)*2  +  K4 


fG(ui)  -  cos  ksg  -  1  sin  ksg 
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(co)  -  VT(0)f T(o)> 


^11 C00)  "  U01-iK)/(o2+K2)]o2fT(a)) 
^12  "  I  (03-iK>/((^+K2)]04fT(ai) 


13  (ai)  -  [NT/DTlfT(cb) 


where 


(85) 


and 


NT  "  (°5°8  +  °6°7)k2  ‘  °6°7  +  i[05°6°7  +  °8(06 
DT  -  0 l  -  (206  -  02)k2  +  k4 


fT(u>)  -  cos  kst  -  i  Sin  kst 


sJJq^  "  ^w(°)fw(^) 

*2l<“)  -  t(Nml-lK)/(H21+K2)]  Nm2£w(oo) 
*22 <“)  *  l(Nm3-lK)/(N^+K2)]  N^fyCo)) 


k2)  Ik 


'  [Nwa/DWa.1£W(‘°> 


where 


«lta  -  (NmsNmS  +  »m6“«7>*2  *  N^6Nm7 


°Wm 


N. 


m6 


+  +  Nm8<Nm6  ‘ 

-  <^*6  -  t&)K2  +  k4 


k2)]k 


and 


J 


(86) 


J 


fw(uj)  -  C08  KSy  -  i  sin  K8y 


V.  Airplane  Structural  Loads 

The  wing  incremental  structural  loads  are  derived 
herein  using  the  mode  acceleration  method,  with  the  aero¬ 
dynamic  forces  being  approximated  by  the  rigid-body  air¬ 
load  distribution  which  would  produce  the  previously 
calculated  rigid-body  vertical  acceleration.  The  incremen¬ 
tal  shear  at  the  tail  root  is  based  on  the  force  summation 
method,  with  the  unsteady  aerodynamic  forces  due  to  motion, 
gust,  and  downwash  being  summed  over  the  tail  panels.  The 
associated  tail  torque  and  bending  moment  are  obtained  by 
multiplying  the  lift  by  appropriate  moment  arms. 

W 

The  final  equations  for  the  wing  total  shear,  S.  (s), 

IT  ^ 

bending  moment,  M^(s),  and  torque  about  an  axis  perpendicu¬ 
lar  to  the  airplane  center  line  at  the  intersection  of  the 
elastic  axis  and  the  i*"*1  spanwise  load  station,  T^(s),  are 

summarized  below,  from  equations  (F66  ),  in  terms  of  the 
aerodynamic  and  inertial  contributions . 


S?(s) 


M?(s) 


T?(s) 


where 


N  OO 

cv  2  at;  £n(s) 

vn-l  ln  n 


N  r.Tjur  OO 

n-i 

N  UT  OO 

Cy  ZlAin  ^n(s) 
n-i 


(g/ga)  *  i 


+  (g/ga)Aia  ?iw(s) 
+  <8/*.>AU  ?»(•> 

+  <8/8«>AU  ?1WC> 


pounds  (mass) 
pounds  (force) 


(87) 


(88) 

(89) 


The  coefficients,  A^,  A^,  and  A^J,  used  in  the  mode- 

acceleration  method  are.  evaluated  from  the  rigid-body 
aerodynamic  loading  distribution  associated  with  an  in¬ 
cremental  one  MgM  vertical  acceleration.  The  acceleration 

OO 

term,  ( 8) ,  is  defined  as  the  effective  airplane  rigid- 

body  acceleration  with  the  tail  aerodynamic  force  removed. 
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This  variable  is  evaluated  by  equation  (F34) . 


w  uu 

^ly(s)  “  l^(s)  +  ^O^t(s)/Mh 


(90) 


where,  from  equation  (F  12 ) , 

N  rp  OO  rp  O  rp 

Vs>  -  2  (A l  «„(»)  +  B*  5  (s)  +  C*en(s)l 

n=l 


+  A*  7  (8)  +  B*  7(8)  +  C^7(S) 

+  Vl|1-lk<S)  +  V2|0[alk(s>  +Alk<s> 


+  '»lk<s>1 


(91) 


The  equations  for  the  tail  shear,  bending  moment,  and 
torque  at  the  intersection  of  the  fuselage  side  and  torque 
tube  (stabilizer  station  14)  also  are  expres&ed  in  terms 
of  aerodynamic  and  inertial  contributions.  These  equa¬ 
tions,  reproduced  from  equation  (F67  ),  are 


N  TS 


Sj(s)  -  C  [J  (s)  +  2  Aj^  5  (s)  +  A™  7(s)1 

n»l  ' 

m[(8)  -  Cv[y0/T(s)  +  2  A™  fn(»)  +  A®  y  (s)]  (93) 

n*l  * 


rp  N  rifP  oo  mM  OO 

t{(8)  -  Cv[x0i  (s)  +  2  A"  en(s)  +  AF  7  (»)] 

n-1  • 


(92) 


(94) 


Introduction  of  simple  harmonic  motion  into  the  wing 
and  tail  load  equations  reduces  equations  (87)  to  (94) 
to  the  following: 

s”((D)  -  -Ac,!  i*  en(»)  +  (g/g,)A^?lw(»)]  (95) 


(96) 


T”(o>)  =  •I!(C!^!n(,)  +  (g/g  .)*S«1W(»)1  (97) 

^iw(cc>  =  ^l(Q)  *  [V0iT(co)]/K2M11  (98) 

N 

Vm)  -  n7=  L[  (Cn  -  kA„)  +  lKB^?n('-“) 

2 

+  [  (C^  -  kV)  4-  iKB^]7<co)  +  2  <J>^k(a>) 

k  =  l 

3 

+  +  Aik(a))  +  ^ik(oo)]  (") 

-SjO  -  cV[iT(a))  "  k2  2  Al^n(a))  "  k2aP7(o))]  (100) 

n»l  ■L' 

Ml(m>  *  Wt^  *  k2  ^,A™5n(i3)  -  k2A™7(-j)))  (101) 

n-1  • 

Tl(j>)  -  Cvtx02T^JJ)  *  K2n^1Aln?n(“>  '  k2a”7(j3)  1  (102> 
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VI.  Airplane  Response  Equations 

The  equations  for  the  normal  and  pitching  accelerations 
and  velocities  at  any  point  on  the  airplane  are  expressed 
in  units  of  real  time  although  they  are  dependent  upon  and 
presented  as  functions  of  non-dimensional  time.  These  and 
the  displacement  equations  are  given  by 


h^s)  -  (V/h^)2  "  t  (V  ?„(.) 

n»l  n 
N 

“!<«>  -  WV2^  *C<8> 

V8>  ■  <V/V  "  fh  <vi>  e„<8> 

n  “1  n 

N 

at(s)  -  (V/l^)  2  f  (v£)  £n(®> 
n*  1  n 

hi<8>  ■  J1thn<’i)tn<8) 

ai<8>  -  vw*> 


(103) 


J 


where  the  subscript  1  is  assigned  integer  values  which  are 
associated  with  particular  points  on  the  aircraft.  For 
example,  i-1  represents  the  cockpit  in  the  present  study. 
The  parameter  denotes  the  location  of  the  1th  response 
point.  For  simple  harmonic  motion,  equations  (103)  reduce 
to  the  following  relationships: 


ht(<i>)  -  -(V/bg)2  r2  X  £  (vi)«n(a>) 

n 

\(m)  -  - (V/b,) 2  K2  S  £  (*.)((«) 

n-1  'h 


(104) 
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N 


h±  <CJL>) 

-  MV/b^ 

s 

n-l 

\ 

(vi)4n(^) 

• 

N 

0^(0)) 

-  iK<V/bj^) 

2 

n-l 

(vp^Co)) 

N 

ht  (i)) 

•  n-l  V 

'l>«n 

t(0>) 

> 

N 

L(o>) 

(104) 
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VII .  Power  Spectral  Techniques 

The  basic  expressions  required  for  the  application  of 
power -spectral  techniques  in  the  present  analysis  are 
developed  and  discussed  in  References  2  through  6  .  In 
this  section,  these  expressions  are  summarized  in  the  no¬ 
tation  of  this  report,  and  the  numerical  evaluation  of  an 
integral  which  occurs  is  discussed.  Two  sets  of  the  basic 
expressions  are  given,  one  in  terms  of  the  frequency,  co,  in 
radians  per  second,  and  the  second  in  terms  of  the  reduced 
frequency,  ft',  in  radians  per  foot,  where 

ft  -  cd/V  (105) 


Hie  power  spectrum  of  the  vertical  gust  velocity,  not 
including  any  spanwise  variation,  is  given  by  equation  (1) 
of  Reference  3  . 


(ft) 


L  1  +  3ft2L2 
T  [1  +  ft2L2]2 


(106) 


where  L  is  the  scale  of  turbulence  in  feet.  Hie  relation¬ 
ship  between  the  power  spectra  as  a  function  of  the  inde¬ 
pendent  variables  ft  and  ao  is  given  by  equation  (15),  Ref¬ 
erence  2  ,  which  is  reproduced  below. 

\  (*)  -  %  (n)/v  (107) 

WG  G 


Combining  equations  (105)  to  (107), 


(*) 


2  JL  1  +  3mV/V^ 

WG  TV  [1  +  'jo2LZ/V2]2 


(108) 


Representing  the  transfer  function  of  the  response 
variable  under  consideration  by  T^(ft)  and  the  output  power 

spectrum  by  $y(ft),  the  expression  relating  these  two  quan¬ 
tities  and  the  gust  power  spectrum  is  given  by  equation  (18), 
Reference  2  ,  as 


y  g 


(109) 
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or  in  terms  of  the  frequency,  a>. 


y--°> 


(110) 


The  response  power-spectral-density  function  may  now 
be  used  to  define  the  characteristic  frequency,  Nq  ,  and 

y 

root-mean- square  value,  fi  ,  of  the  particular  variable  in 
question.  These  are  given  in  Reference  3  as 


/o^2$vC-a)d^ 
/  o*y(to)d'-o 


(111) 


6y  =  /jFy(fl>)dD  (112) 

The  constant  Nq  gives  the  expected  number  of  times  per 

y 

second  that  the  variable  y  crosses  the  value  zero  with  a 
positive  slope. 

A  more  general  parameter,  N(y),  which  is  dependent 
on  Nq  ,  is  defined  as  the  expected  number  of  times  per 

y 

second  that  the  variable  crosses  some  finite  value  of  y 
with  a  positive  slope.  This  parameter  is  also  an  approxi¬ 
mation  for  the  average  number  of  peaks  per  second  exceed¬ 
ing  a  given  value  of  y,  and  will  be  used  in  this  sense 
throughout  this  report.  An  analytical  representation  of 
N(y)  is  given  by  equation  (21)  in  Reference  6  .  This  ex¬ 
pression  contains  two  terms  to  account  for  both  storm  and 
non-storm  turbulence.  For  the  present  study,  which  is 
restricted  to  low-altitude  flight,  the  second  term  asso¬ 
ciated  with  storm  turbulence  is  neglected.  The  expression 
for  N(y)  then  becomes: 

-  y 

N(y)  -  PjNq  e  B^y/6wG)  (113) 

y 

where  is  the  proportion  of  total  flight  time  in  non¬ 
storm  turbulence,  and  B  is  a  scale  parameter  in  the  proba¬ 
bility  distribution  of  root-mean- square  gust  velocity. 
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Note  that  for  y  -  0  and  ■  1.0,  equation  (113)  reduces 
to  N(0)  -  N0  . 

y 

At  this  point,  it  is  noted  that,  since  equations  (111) 
and  (112)  require  integration  from  zero  to  infinity,  the 
values  of  the  integrals  may  become  infinite  if  the  response 
power  spectrum  does  not  decay  fast  enough.  In  order  to  in¬ 
vestigate  these  integrals  further,  equation  (110)  is  intro¬ 
duced  and  the  integrals  are  written  in  general  form  as 


(114) 


Also,  the  expression  for  the  gust  power  spectrum  (equa¬ 
tion  (108)),  is  approximated  as  follows: 


(i>)  -  S  W  »  for  '■»  small 
G 


7 r  V 

(2 


$  (jo)  -  — —  ,  for  jo  large 

w. 


TT  L'JO 


(115) 

(116) 


Equating  the  above  two  expressions  results  in  a  "break¬ 
point"  frequency. 


vl>  -  ^ ^  (in) 

A  typical  gust  power  spectrum,  as  given  by  equation  (108), 
and  its  approximation,  as  given  by  equations  (115)  to  (117), 
are  shown  in  Figure  1. 

Combining  equations  (114)  to  (117) , 

0D2k|Ty(jD)|  2do> 

oo2  (k-1)  |xy(jD)|  2  dec]  (118) 


For  some  reasonable  value  of  the  scale  of  turbulence,  the 
first  integral  of  equation  (118)  always  possesses  a  finite 
value  if  the  transfer  function  T  (oo)  is  bounded.  For  the 
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second  integral  of  equation  (118),  a  stronger  restriction 
must  be  associated  with  the  function  T  (jo);  that  is,  it 

must  have  a  certain  rate  of  decay  at  t&e  high  frequencies. 
The  assumption  is  now  made  that  at  the  high  frequencies  the 
transfer  function  depends  upon  the  frequency  in  the  follow¬ 
ing  manner : 


(119) 


Introducing  equation  (119)  into  the  second  integral  of 
equation  (118)  yields 


f  J2(k-!)  +  2z)  ^ 

“b 


a)[2k-l+2e] 


2k- 1+2  e 


“b 


(120) 


Therefore,  in  order  for  the  integral  In  equation  (114)  to 
have  a  finite  value,  the  following  restriction  must  be 
imposed: 


2k- 1+2  e  <  0 


(121) 


For  the  response  parameters  considered  in  the  present 
study,  i.e.,  displacements,  velocities,  accelerations,  and 
structural  loads,  an  upper  bound  to  the  value  of  the  ex¬ 
ponent  in  equation  (119)  is  s  -  0.  For  certain  other 
response  parameters,  such  as  the  rate  of  change  of  accelera¬ 
tion  (jerk),  e  could  be  positive,  but  these  will  not  be 
considered  herein.  For  e  -  0  then, equation  (120)  becomes 


/“  .2<k-1>  |  T  (cjd)|  2doo  -  co2k”1/(2k-l) 

rw.  y 


(122) 


Thus,  equation  (121)  reduces  to 

2k- 1  <  0 


(123) 
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which,  since  k  is  an  integer,  becomes 

k  1  0 

Introducing  the  requirement  of  equation  (124)  into 
equation  (114)  and  then  comparing  this  equation  with 
equations  (111)  and  (112) ,  it  is  seen  that  the  integration 
for  the  root-mean-square  response,  6^,  is  bounded  whereas 

that  for  the  characteristic  frequency,  Nq  ,  may  become  in- 

y 

finite.  This  last  conclusion  is  of  course  inconsistent 
with  realistic  considerations,  and  to  eliminate  this  con¬ 
dition  it  has  been  proposed  in  Reference  4  that  the 
expression  for  the  gust  power  spectrum  as  given  by  equation 
(108)  be  modified  so  that  theoretical  response  peak  calcu¬ 
lations  become  consistent  with  flight  test  data.  This 
modified  gust  power  spectrum  is 

1  +  3'jo2L2/V2  C4  +  6C2'jd2L2/V2  -  3  jd4L4/V4 

[1  +  oo2L2/V2]3  +  [C2  +  o>ZL2/V2]3 


*  wr 

«  (CD)  «  — 2- 
G  t rV 


(125) 

A  value  of  50  for  the  parameter  C  in  this  equation  has  been 
proposed  in  the  above  reference. 

With  the  aid  of  equation  (125),  an  investigation  was 
made  to  determine  an  acceptable  upper  limit  of  integration 
for  the  response  power  spectra  being  considered  in  this 
study.  It  was  concluded  that  the  use  of  an  upper  limit  of 
200  rad. /sec.  in  conjunction  with  the  gust  power  spectrum 
of  equation  (108)  would  yield  characteristic  frequencies 
which  are  in  good  agreement  with  those  obtained  by  the  use 
of  an  effectively  infinite  limit  of  integration  in  conjunc¬ 
tion  with  the  power  spectrum  of  equation  (125) .  Also,  this 
upper  limit  of  200  rad. /sec.  was  found  to  be  more  than 
sufficient  for  accurate  calculation  of  root-mean- square 
response  values. 
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Figure  1 

GUST  POWER-SPECTRAL-DENSITY  FUNCTION 

AND  SIMPLIFIED  APPROXIMATION  FOR  6.  -  1  FT. /SEC. 

— - - - - "  ■  1 


L  ■  500  Ft.. 


1000  Ft. /Sec. 


I 

\ 

CM 

$ 


|e 


l 

co 

1 

ii 

« 

s 


fg)  to  (117) 


.1 


10 


10“ 


Frequency  -  gd,  R/T 


REPORT 

DATE 


ADR  06-14-63.1 
OCTOBER  1963 


GRUMMAN  AIRCRAFT  ENGINEERING  CORPORATION 


44 


APPENDIX  A 

AUTOPILOT  AND  CONTROL  SYSTEM  EQUATIONS 


I .  Autopilot  Signal 

A  mathematical  description  of  the  A-6A  (A2F-1)  auto¬ 
pilot  and  its  related  systems  is  given  in  Reference  7 
In  particular,  Figures  2  and  3  of  Reference  7  provide 
block  diagrams  for  the  autopilot  system  in  terms  of  the 
five  following  autopilot  modes :  stability  augmentation, 
attitude  hold,  Mach  hold,  altitude  hold,  and  G  command. 

For  the  present  gust  response  study,  only  two  of  these 
autopilot  modes  are  included  in  the  analysis,  namely, 
stability  augmentation  and  attitude  hold.  Considering  only 
these  two  modes,  the  block  diagrams  in  Reference  7  are 
translated  into  the  following  equation  in  terms  of  Laplace 
transform  operators : 


TT(p)  *  [2pk7/(2p  +  l)]£[a(x“,t)]  +  kg££a(x“  t)]  (Al) 


As  written,  equation  (Al)  is  the  transform  equation  of  the 
autopilot  system  for  the  attitude  hold  mode.  It  may  be 
reduced  to  the  stability  augmentation  mode  by  simply 
equating  kg  to  zero.  The  values  of  the  constants  obtained 

from  the  block  diagram  in  Reference  7  are 


j  (A2) 

kg  “  1,  R/R 


The  autopilot  input  is  now  expressed  in  terms  of  the 
normal  modes  and  mode  shapes  at  the  location  of  the  various 
sensing  devices  as  follows: 


»<*£.  *>  *  l  £an<X^n(C) 
a<V  e>  '  l 


(A3) 
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Also,  to  be  consistent  with  other  equations  to  be  developed 
subsequently,  the  real  time,  t,  is  replaced  by  the  non- 
dimensional  time,  s,  by  use  of  the  following  relationships: 


s  -  (V/bR) t 


d/dt  =  (V/t^Hd/ds) 
d2/dt2  =  (V/bR) 2 (d2/ ds2) 


(A4) 


J 


Substituting  equation  (A3)  into  equation  (Al) , 
and  then  taking  the  inverse  Laplace  transform  of  this 
expression  yields  the  autopilot  equation  as  a  function  of 
real  time,  t.  Conversion  to  non-dimensional  time  may  now 
be  accomplished  through  use  of  equations  (A4) .  The  final 
expression  for  the  autopilot  signal  is 


o  oo 

M-(s)  +  L10n(s)  “  S'Lln«„<s>  + 

n 

where  the  coefficients  in  equation 
Lio  =  bR/2V 
Lln  -  (V/bR)k7f0n(x“) 

L2n  “  k8fan  <*£> 

L3n  *  (V2V>k8fan(<) 


L2J„<S>  +  L3A<s>]  <A5> 

(A  5)  are  given  by 

\ 

(A6) 
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APPENDIX  A 


II .  Stick  and  Sprashpot  Signals 

A  complete  description  of  the  A-6A  manual  control 
system,  together  with  numerical  values  of  the  properties 
of  the  various  components,  is  presented  in  Reference  8, 
Reference  9  summarizes  the  basic  data  presented  in  Refer¬ 
ence  8  ,  and  also  presents  equations  of  motion  for  the 

rigid-body  airplane  together  with  the  expressions  for  the 
stick  and  sprashpot  systems  as  a  function  of  the  rigid 
airplane  vertical  and  angular  accelerations.  In  the 
present  analysis,  airplane  flexibility  will  also  be  con¬ 
sidered.  Therefore,  the  stick  and  sprashpot  equations 
are  rederived  in  this  section  in  terms  of  the  airplane 
normal  modes.  The  basic  derivation  follows  the  procedure 
presented  in  Reference  10,  pages  22  to  24,  which  applies 
to  a  rigid  airplane. 

A  schematic  diagram  of  the  manual  control  system 
(stick  and  sprashpot)  is  shown  in  Figure  Al,  which  is  a 
modification  of  Figure  2.1  in  Reference  8  .  The  control 
system  without  bobweights  has  a  basic  weight  moment  of 
inertia  about  the  stick  pivot  point  denoted  by  IgQ*  The 

expression  for  the  additional  control  system  inertia  due 
to  bobweights,  1^,  is  now  obtained.  The  effective  kinetic 

energy  of  the  stick  may  be  related  to  the  kinetic  energy 
of  the  bobweight  system  through  the  following  expression. 

i(lsl/g)T(e)2  =  ^(la/g)(NaT(t))2 

+  %(If/g)(Nft(t))2  (A7) 

where  and  are  the  ratios  of  the  bobweight  angular 

motions  to  control  stick  angular  motion  for  the  aft  and 
forward  bobweights,  respectively.  The  bobweight  inertias 
about  their  own  pivot  points,  disregarding  the  contribu¬ 
tions  about  their  own  centers  of  gravity,  are  given  by 


(A8) 


Combining  equations  (A7)  and  (A8)  will  provide  the  moment 
of  inertia  about  the  stick  pivot  point. 


a 


maRa 
a  a 


mfRl 
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I.l  -  mfNfRf  +  maNaRa 


(A9) 


The  total  control  moment  inertia  referred  to  the  stick 
pivot  point  is  therefore  given  by 


I.  * 


Is0  +  Isl 


s  0 


2  2 
mfN|R| 


2  2 

Wi 


(A10) 


It  is  assumed  in  the  above  derivation  that  the  forward 
and  aft  bobweights  are  the  only  control  system  components 
sensing  airplane  acceleration.  The  further  assumptions 
made  are  that  normal  (positive  down)  and  pitch  (positive 
nose-up)  acceleration  sensing  are  not  dependent  upon  stick 
position,  and  that  forces  normal  to  the  flight  path'  are 
very  nearly  perpendicular  to  the  airplane  fuselage  ref¬ 
erence  line. 

The  moment  about  the  stick  pivot  point  due  to  airplane 
and  control  system  accelerations,  M^T(t),  may  now  be 

written  as  follows: 

MlT(t)  -  (Is/g)lV(t)  +a(xg,t)] 


-  (mfRfNf/g)[h(xf,t)  -  Rfa(xf,t)] 

+  (maRaNa//g)[l^xa,t)  +Ra“0ca,t0]  <A11) 


where  the  indicated  normal  and  angular  accelerations  about 
the  stick  pivot  point  and  the  forward  and  aft  bobweight 
pivot  points  are  given  in  terms  of  the  normal  modes  by 


“<*«•*>  - 

a(xf,t)  -  2f  (xf)£  (t) 
n  n 


h(xf,t)  -  2fh^(Xf )’^n(t) 


(A12) 


h(Xa,t)  -  2fhn(x,)en(t) 
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Additional  moments  about  the  stick  pivot  point  resulting 
from  sources  other  than  the  accelerations  are  due  to  the 
sprashpot  and  artificial  feel  spring.  Under  the  assump¬ 
tion  that  the  spring  moment  is  proportional  to  the  stick 
angular  displacement,  and  the  sprashpot  contribution  is 
proportional  to  the  sprashpot  angular  velocity  measured 
at  a  point  between  the  spring  and  the  damper,  the  addi¬ 
tional  moment  on  the  stick  is  given  by 

M2r(t)  «  KxT(t)  +  CyKt)  (A13) 

where  Kt  and  are  the  angular  spring  constant  and  angu¬ 
lar  damping  constant  about  the  stick  pivot  point  of  the 
artificial  feel  spring  and  sprashpot,  respectively.  The 
total  moment  about  the  stick  pivot  point,  without  any  ex¬ 
ternal  forces  such  as  pilot  input,  is  given  by 

MlT(t)  +  M2x(t)  =  0  (A14) 

Combining  equations  (All)  and  (A13) , 

(ls/g) [  V(t)  +  a(xg,t)]  -  (mfRfNf/g)[h(xf,t)-Rfa(xf,t)] 


+  (maRaNa//g^ [  Mxa,  +  Ra“^xa»  ^ 1 

+  KTx(t)  +  C^nCt)  *  0  (A15) 

The  expression  relating  the  stick  angular  displacement 
and  the  sprashpot  equivalent  angular  motion  is  given  by 

C^(t)  +  K^txiCt)  -  x(t>]  -  0  (A16) 

where  is  the  angular  spring  constant  of  the  sprashpot 
about  the  stick  pivot  point. 

After  introducing  the  non-dimensional  time  and  making 
the  appropriate  substitutions  into  equations  (A15)  and 
(A16),  the  final  expressions  for  the  stick  and  sprashpot 
become 
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00  6  oo  ° 

T(s)  +  h30t(s)  -nSH3nln(s)  +  H3tln(«)  (A17) 

and 

n(s)  +  H40ti(s)  =  H4Qt(s)  (A18) 


where  the  coefficients  are  given  by 

h30'<W(VV>2« 

t 

h3t,  -  -«yis><bR/v>g 


H 


3n 


-  -<i/i8> 


IsO£cn(xs) 


(A19) 


+  ™fRfNf[Rf(Nff  (x  )  +  f  (xf))  -  £h  (xf)] 

n  n  n 

+  %RaNa[Ra(Na£  (xs)  +  £  (xa))  +  fh(xfl)  1 

n  n  n 


(A20) 


H40  *  <YV ( YV> 


(A21) 
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III.  Stabilizer  Rotation 


A  block  diagram  of  the  A-6A  (A2F-1)  longitudinal 
power  control  system  in  the  manual  and  series  modes  is 
shown  in  Figure  4.2  of  Reference  11.  For  the  frequency 
range  of  interest  in  the  present  gust  response  studies, 
it  would  have  been  sufficient,  from  the  standpoint  of 
the  stabilizer  response  characteristics,  to  approximate 
the  power  control  system  by  a  first-order  differential 
equation  representing  the  linear  factor  of  the  power 
control  actuator  dynamics.  For* the  calculation  of 
inertia  loads,  however,  it  was  required  to  include  a  sec¬ 
ond  order  or  angular  acceleration  term.  For  this  purpose, 
the  power-control- system  transfer  function  was  represented 
by  the  expression 

7<P)/7S(P)  *  l/(TlP+l)(T2p+l)  (A22) 

In  this  expression,  the  dynamics  of  the  secondary  ac¬ 
tuator  shown  in  Figure  4.2  of  Reference  11  are 
neglected.  The  first  time  constant,  Tp  is  the  same  as 

that  inherent  in  the  first  factor  of  the  original  third- 
order  power-control-actuator  equation.  The  second  time 
constant,  T2,  is  an  approximation  which  reduces  the  third- 

order  power-control-actuator  equation  to  a  second-order 
one. 


The  actual  differential  equation  for  stabilizer  ro¬ 
tation  used  in  the  present  study  is  given  as 

a77(t)  +  b7*y(t)  +  7(t)  -  7s(t)  (A23) 

where 

ay  «  V2  (A24) 

by  -  Tx  +  T2  (A25) 

and  7_(t)  is  the  total  signal  input  to  this  equation. 

O 

The  input  consists  of  two  signals,  one,  due  to  the  auto¬ 
pilot,  p.(t),  and  the  other  due  to  the  stick  including  the 
stabilizer-to-control-stick  linkage  ratio,  k  r(t).  The 


51 


APPENDIX  A 


final  form  of  the  expression  describing  the  stabilizer 
rotation  is  given  in  terms  of  non-dimensional  time  as 


7(s)  +  H207(s)  +  H217(s)  -  Hnn(s)  +  H22t(s)  (A26) 

where  the  coefficients  are  defined  by 

H20  -  <VS7)  (VV)  (A27) 

H21  ’  (V«7)(VV>2  (A28) 

h22  *  (kT/a7)(VV)2  (A29) 
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UNSTEADY  AERODYNAMICS  OF  RECTANGULAR  PLANFORMS 
DUE  TO  GUST,  MOTION,  AND  DOWNWASH 
(Indie la 1  Lift  Functions) 


I.  Indicial  Lift  Functions  for  Gust  on  a  Planform  (Ktissner- 
Type  Functions) 

A.  Planform  Ktissner-Type  Function 

The  idealized  model  of  a  two-dimensional  wing- tail 
combination  flying  at  a  constant  subsonic  Mach  number,  M, 
and  encountering  a  stationary  vertical  gust  with  a  verti¬ 
cal  gust  velocity  wG  (wG  «  V)  is  represented  in  Figure  Bl. 

Tail  position  with  respect  to  the  wing  (defined  by  the 
parameters  s*  and  s^)  is  important  in  evaluating  lift  build 

up  on  the  tail  due  to  wing  downwash.  Downwash  effects, 
whether  due  to  wing  motion  or  gust  on  the  wing,  will  be 
discussed  in  Section  III  of  this  appendix. 

The  Kussner-type  function  is  defined  herein  as  the  time 
variation  of  lift  for  an  airfoil  of  arbitrary  aspect  ratio 
and  flying  at  an  arbitrary  subsonic  Mach  number,  due  to  a 
step  change  in  vertical  gust  velocity,  normalized  with 
respect  to  the  steady-state  lift  at  infinite  time.  (The 
Kussner  function  itself  is  restricted  to  apply  only  to  the 
normalized  lift  build  up  due  to  a  step  gust  on  a  two- 
dimensional  airfoil  (AR  *  °°)  in  incompressible  flow  (M  *  0)  . 
The  Kussner-type  function  is  dependent  on  a  non-dimensional 
time,  t),  defined  in  terms  of  planform  (wing  or  tail) 
reference  semichords  traveled. 

n  -  Vt/b^  (Bl) 


where  V  is  the  airplane  forward  velocity  and  b^  is  a 

reference  semichord  for  the  planform  under  consideration. 
For  an  untapered  planform,  b^  would  be  identically  equal 


to  half  the  geometric  chord,  whereas  for  a  tapered  plan- 
form  it  would  be  half  the  chord  at  a  reference  spanwise 
station  such  as  the  mean  aerodynamic  chord.  The  Kussner 
type  function  increases  from  an  initial  value  of  zero  to 
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a  steady-state  value  of  unity  at  tj  equal  to  infinity.  As  .. 

shown  in  Figure  B2,  the  effect  of  Mach  number  on  the  basic 
Kussner  function  is  to  decrease  the  rate  of  build  up, 
whereas  finite  aspect  ratio  has  the  opposite  effect. 

•• 

Kussner-type  functions  for  two-dimensional  airfoils  at 
Mach  numbers  other  than  zero  and  for  finite-aspect-ratio 
airfoils  in  incompressible  flow  have  been  derived  in  vari¬ 
ous  technical  reports .  (See  References  12  ,  13  ,  14.)  By 
interpolation  or  extrapolation  .from  these  two  groups, 
values  of  the  Kussner-type  functions  for  any  Mach  number 
and  aspect  ratio,  other  than  the  given  curves,  may  be  es¬ 
timated.  The  method  used  in  the  present  investigation  is  i 

based  on  the  following  assumed  formula: 

*<M,  AR,  q)  -  00  >  Ti)  (B2) 


where  °°,  q)  is  the  incompressible  two-dimensional 

(Kussner)  function  (See  Reference  15,  page  288)  and 
^(0,  AR,  q)  and  ^(M,  °°,  q)  are  incompressible  three- 
dimensional,  and  compressible  two-dimensional  Kussner- 
type  functions  such  as  those  evaluated  in  the  references 
stated  previously. 

Equation  (B2)  as  given  provides  a  time  relationship 
for  the  indicial  lift  function  based  upon  a  non-dimensional 
time  measured  in  planform  semichords.  Subsequent  to  evalu¬ 
ating  this  function  for  the  given  conditions,  it  is  desirable 
to  approximate  the  results  by  simple  analytical  expressions 
to  facilitate  the  use  of  indicial  lift  functions  in  the 
present  study.  A  simple  and  useful  approximation  is  an 
expression  consisting  of  exponential  terms,  where  the  num¬ 
ber  of  terms  is  dictated  by  the  desired  accuracy  of  the 
approximation.  The  assumed  expression  then  is  of  the 
following  general  form: 

-Vi 

Hr\)  -  *(0)  +  2  ak(l  -  e  )  (B3) 

k-1 

Note  that,  to  obtain  the  correct  steady-state  lift,  the 
following  condition  must  be  imposed  upon  the  coefficients 

3k‘  % 

*(0)  +2  a.  -  1.0  (B4) 

k-1  k 
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The  function  if/(r\)  provides  the  analytical  approxima¬ 
tion  of  the  indicial  lift  as  calculated  from  the  expression 
ip (M,  AR,  rj)  given  by  equation  (B2)  .  The  maximum  index, 

K^,  appearing  in  equation  (B3)  designates  the  number  of 

exponentials  used  in  approximating  the  planform  indicial 
lift  function.  The  purpose  of  the  constant,  ^(0),  is  to 
permit  the  reduction  of  the  equation  to  the  limiting  case 
for  which  the  lift  build  up  due  to  gust  is  assumed  in¬ 
stantaneous,  that  is,  no  transient  effects  are  present. 

The  reduction  of  equation  (B3)  to  instantaneous  lift 
build  up  may  be  accomplished  simply  by  letting  ^(0)  =  1.0 
and  all  aR  =  0. 


Evaluation  of  the  constants  appearing  in  equation  (B3) 
may  be  carried  out  by  a  graphical  procedure  as  follows. 

First  the  variable  1  -  ^(M,  AR,  T|)  is  plotted  on  semi-log 
paper  versus  the  non-dimensional  time  for  some  finite  value 
of  T).  Any  straight  line  drawn  on  the  s etni - 1  ogpa per  repre- 

-cyi 

sents  one  exponential  term,  as  defined  by  a,  e  ;  therefore, 


the  constants  aR 


and  ct^  are  easily  evaluated. 


A  number  of 


these  straight  lines  whose  sum  is  equal  to  1  -  ^(M,  AR,  t]) 
will  represent  the  variable  (1  -  ^(q))  with  the  desired 
accuracy. 


The  application  of  these  exponential  expressions  to  a 
given  wing-tail  combination  is  summarized  below. 

B.  Explicit  Expressions  for  the  Wing  and  Tail  Kussner- 
Type  Functions 

For  a  specific  wing-tail  combination,  the  real  time, 
t,  is  non-dimensional  with  respect  to  an  arbitrary  distance, 
bR,  through  the  relationship, 

s  *  Vt/bR  (B5) 


In  the  present  analysis,  the  planform  reference  semichord, 
b^,  for  the  wing  was  chosen  as  half  the  wing  mean  aero¬ 
dynamic  chord,  bw,  and  similarly,  b^  for  the  tail  was  chosen 
as  half  the  tail  mean  aerodynamic  chord,  b^.  Thus,  for  the 

purposes  of  evaluating  the  Kussner-type  functions,  the 
actual  swept  wing-tail  combination  was  replaced  by  a  rec¬ 
tangular  wing-tail  combination  as  shown  in  Figure  B3. 
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The  arbitrary  reference  distance,  bR,  used  in  non-dimen- 

sionalizing  the  time  variable  was  chosen  in  the  present 
analysis  as  bT^.  The  relationship  between  the  non- 

dimensional  time  parameters  defined  by  equations  (Bl)  and 
(B5)  is  obtained  by  combining  these  equations. 

ti  =  (bR/b^)s  (B6) 

For  the  gust  on  the  wing,  the  expression  for  the 

Kussner-type  function  is  given  by  an  expression  similar 
to  equation  (B3)  where  the  time  r\  is  replaced  by  s  and 
the  subscript  and  superscript,  W,  designates  the  wing. 

W 

*  %  u  ~CKS 

fy(s)  -  ^y(0)  +  2  a£(l-e  K  )  (B7) 

k — 1 


This  result  is  obtained  by  substituting  equation  (B6) 
into  equation  (B3),  and,  as  discussed  above,  setting 


(B8) 


w 

Therefore,  the  new  exponent  is  defined  by  the  follow¬ 
ing  expression. 


\  -W 


bw  -W 


°k 


(B9) 


*  Note:  In  order  to  keep  the  number  of  symbols  used  to  a 
minimum,  the  left  hand  side  of  equation  (B3)  will 
be  represented  by  ^w(s),  after  substitution  of 

equation  (B6)  into  equation  (B3) ,  rather  than 

bR 

tf'uCjT"  s) .  Note  that  this  procedure  in  represent - 

n 

ing  the  functional  relationships  will  be  used 
throughout  the  analysis  whenever  it  is  necessary 
to  switch  from  one  independent  variable  to 
another. 
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As  noted  earlier,  the  values  for  are  evaluated  based 

upon  rj,  and  are  dependent  on  the  wing  aspect  ratio  and 
airplane  Mach  number. 

For  the  gust  on  the  tail,  the  Kussner-type  function 

is  similar  to  equation  (B7)  with  the  subscript  or  super¬ 
script,  W,  replaced  by  T. 


^T(s) 


*T<0> 


k  ;  _  -cJs 

+  2  a£(l-e  ^  > 
k  =1  k 


(BIO) 


In  this  case  the  semichord  b 


is  defined  by 


=  bT  (Bll) 

Once  again,  bR  =  bw,  since  the  same  nondimensionalizing 

reference  distance  must  be  used  throughout  the  analysis. 
The  new  exponent  for  the  tail  is  identified  by  the  super¬ 
script,  T,  and  is  given  by 

„T  _  hR  _T  _  hf  _T 

**  S  ^  {B12) 

As  before,  the  values  of  in  equation  (B12)  are  defined 

in  terms  of  the  non-dimensional  time,  q,  and  depend  on 
the  tail  aspect  ratio  and  airplane  Mach  number. 
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II.  Indlclal  Lift  Functions  Due  to  Planform  Motion 

(Wagner-Type  Functions) 

A.  Planform  Wagner-Type  Function 

The  following  discussion  of  the  Wagner-type  function 
is  similar  to  that  in  the  previous  section  on  the  Kiissner- 
type  function.  The  Wagner- type  function  is  defined  herein 
as  the  time  variation  of  lift  for  an  airfoil  of  arbitrary 
aspect  ratio  and  flying  at  an  arbitrary  subsonic  Mach  num¬ 
ber  due  to  motion  resulting  in  a  step  change  in  angle  of 
attack,  normalized  with  respect  to  the  steady-state  lift 
at  infinite  time.  The  term  Wagner  function  as  used  here 
is  restricted  to  the  special  case  of  the  Wagner-type  func¬ 
tion  in  which  a  two  dimensional  airfoil  is  flying  in  in¬ 
compressible  flow. 

In  contrast  to  the  Kiissner-type  function,  the  Wagner- 
type  function  exhibits  an  instantaneous  partial  lift  and 
thereafter  a  build  up  to  the  steady-state  value  of  unity. 
As  shown  in  Figure  B4,  the  Wagner-type  function  is  greater 
than  the  Kussner-type  function  for  all  values  of  the  non- 
dimensional  time,  T],  but  the  functions  approach  each 
other  asymptotically  as  q  approaches  infinity. 

The  Wagner-type  function  will  be  evaluated  by  an  ex¬ 
pression  similar  to  equation  (B2) . 


0(M,  AR,  q)  -  0(M’  °°»  ^  (B13) 


where  0(0,  <»,  q)  is  the  incompressible  two-dimensional 
(Wagner)  function  and  0(0,  AR,  q)  and  0(M,  °°,  q)  are 
special  cases  of  the  Wagner-type  function  such  as  those 
found  in  References  12,  13,  and  14. 

The  assumed  approximation  for  the  Wagner-type  func¬ 
tion  is  similar  to  equation  (B3). 

-K0 

0 (q)  -  0(0)  +  2  b,(l  -  e  )  (B14) 

k-1  K 


where  now  the  constant  0(0)  possesses  a  positive 
value  less  than  one,  and  in  the  case  of  two-dimensional 
incompressible  flow  its  value  is  0.5. 


APPENDIX  B 


B.  Explicit  Expressions  for  the  Wing  and  Tall 
Wagner-Type  Functions 

For  wing  motion,  the  expression  representing  the  wing 
Wagner -type  function  for  any  Mach  number  is  given  by  a 
relationship  similar  to  equation  (B7) . 

K“  H  A* 

t^(s)  -  4^(0)  +  £  b“(l  -  e  ^  )  (B15) 


W 

where  the  exponent,  is  defined  in  a  manner  similar 
to  equation  (B9) . 


«w  \ 

*  S 


(B16) 


For  tail  motion,  the  expressions  are  similar  to 
equations  (fcl5)  and  (B16)  with  the  letter  W  replaced  by  T. 


K 


0T(s) 


0 


0  (0)  +  2  b£(l 
1  k«l  K 


e  ) 


A 


>  qT 


(B17) 

(B18) 


The  exponent  Based  upon  the  non-dimensional 

time,  q,  and  is  a  function  of  the  planform  aspect  ratio 
and  Mach  number.  Since  the  wing  and  tail  have  different 
aspect  ratios,  the  values  of  ^  differ  for  the  two 

surfaces . 
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III .  Indicial  Lift  Functions  Due  To  Downwash  on  the  Tail 
Caused  By  Wing  Motion  and  Gust  on  the  Wing  (Hobbs - 
Type  Functions) 

A.  Indicial  Downwash  Functions 


Hobbs  in  Reference  1  has  developed  expressions  for 
calculating  the  downwash  variation  on  the  tail  in  incom¬ 
pressible  flow  due  to  a  unit  step  vertical  gust  moving  at 
various  gust  front  velocities  ranging  from  zero  (for  a 
stationary  gust  field  acting  on  the  wing)  to  infinity 
(which  is  equivalent  to  wing  motion) .  For  the  purpose  of 
easy  reference  to  these  functions,  the  descriptive  state¬ 
ment  indicial  downwash  functions  will  be  used  to  describe 
the  time  variation  of  downwash  on  the  tail  normalized 
with  respect  to  the  steady-state  value.  Furthermore,  the 
downwash  effects  may  be  divided  into  downwash  caused  by 
wing  motion  and  downwash  caused  by  gust  on  the  wing,  and 
hence,  the  indicial  downwash  functions  will  be  referred  to 
as  indicial  motion  downwash  functions  and  indicial  gust 
downwash  functions,  respectively. 

The, indicial  downwash  functions  are  of  the  general 
form  presented  £n  Figure  B5*  where  the  variable 
w(Mg,  M,  /Rw,  sT,  s£,  s)  represents  the  indicial 

downwash  as  a  function  of  the  gust  front  Mach  number 
(either  zero  or  infinity  for  most  applications),  airplane 
flight  Mach  number,  wing  aspect  ratio,  horizontal  and 
vertical  positions  of  the  tail  with  respect  to 
the  wing,  and  non-dimensional  time.  The  locations  of  the 
negative  peak  values  in  the  indicial  downwash  functions 
depend  upon  the  parameter,  s£,  which  defines  the  relative 
position  of  the  stabilizer  and  the  wing  in  the  longitudi¬ 
nal  direction  (see  Figure  B6) .  For  the  downwash  due  to 
motion,  this  peak  location  will  occur  just  prior  to  a 
time  given  by  (s*  -  1)  whereas  for  the  downwash  due  to 

gust,  it  occurs  just  prior  to  a  time  defined  by  (s£  +  1). 


*Note:  For  the  present  study,  Hobbs’  curves  are  normal¬ 
ized  with  respect  to  the  steady-state  downwash 
slope,  de/da,  and  the  actual  value  of  de/da  is 
included  in  the  final  coefficients  of  the  general 
ized  forces  which  are  the  inputs  to  the  equations 
of  motion. 
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The  magnitude  of  the  peak  value  depends  on  the  height  of 
the  horizontal  stabilizer  above  the  wing,  s^,,  and  its 

maximum  value  is  obtained  when  both  surfaces  lie  in  the 
same  plane  (s^,  =  0)  . 

For  the  present  study,  unsteady  downwash  functions 
similar  to  those  shown  in  Figure  B5  are  evaluated  from 
either  the  expressions  or  the  final  curves  which  Hobbs 
has  obtained.  The  effects  of  airplane  Mach  number  on  the 
indicial  downwash  functions  and  downwash  variation  across 
the  tail  span  have  been  neglected  in  the  present  analysis. 

B.  Explicit  Expressions  for  the  Hobbs-Type  Functions 

The  unsteady  downwash  experienced  at  the  tail  has 
the  characteristics  of  gust  regardless  of  whether  it 
arises  from  wing  motion  or  gust  on  the  wing.  Therefore, 
lift  build  up  on  the  tail  due  to  downwash  is  evaluated 


will  be  used  to  represent  the  normalized  lift  build  up  on 
the  tail  associated  with  the  unsteady  downwash  due  to  a 
step  gust  on  the  wing  or  a  step  change  in  wing  angle  of 
attack.  This  variable  will  be  referred  to  as  the  Hobbs- 
type  function,  and  the  words  motion  or  gust  will  be 
added  as  appropriate  to  designate  downwash  on  the  tail 
due  to  wing  motion,  or  downwash  on  the  tail  due  to  gust 
on  the  wing,  respectively.  In  terms  of  Duhamel's  integral. 


6(Mg,  M,  AR^,  ARt,  s*  s*,  s)  - 

/*  w(Mg,  0,  /Rw,  S.p,  S,p,  <5)^(s-d)d<j 


(B19) 


Note  that  the  expression  for  ^(s),  as  given  by  equation  (BIO), 

is  a  function  of  airplane  Mach  number,  M,  and  tail  aspect 
ratio,  ARt>  Therefore,  in  spite  of  the  fact  that  the 
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indicial  downwash  function  as  shown  in  equation  (B19)  is 
independent  of  airplane  Mach  number  and  tail  aspect  ratio, 
the  Hobbs-type  function  is  dependent  on  these  parameters. 

The  effect  of  superimposing  the  unsteady  lift  build  up  on 
the  indicial  downwash  is  shown  by  the  curves  in  Figure  B7. 

When  these  curves  are  compared  to  the  indicial  downwash 
functions  of  Figure  B5,  it  is  seen  that  the  unsteady  lift 
build  up  has  the  effect  of  shifting  the  entire  curves  to 
the  right  (including  their  intersection  with  the  horizon¬ 
tal  axis) ,  and  also  reducing  the  ordinates  for  the  initial 
values  of  s. 

The  expression  for  the  normalized  lift  build  up  on  the 
tail  due  to  a  step  gust  on  the  wing  (Hobbs-tvpe  gust  func¬ 
tion)  will  be  obtained  from  equation  (Bl9)  with  Mg  *  0. 

If  the  functional  variable,  6(0,  M,  ART>  si£,  s^,,  s) , 

is  replaced  by  6g(s)  where  the  subscript,  G,  designates 

gust  on  the  wing,  the  result  is 

6G(s)  "  /q  W(0>  °*  *W»  ST>  sT»  0>^iCs-c5)dc5  (B20) 


Note  that  downwash  due  to  gust  on  the  wing  is  assumed  to 
commence  when  the  gust  reaches  a  reference  point  on  the 
wing  defined  by  the  intersection  of  the  wing  leading 
edge  and  the  spanwise  location  of  the  wing  mean  aerody¬ 
namic  chord  (see  Figure  B6) . 

The  expression  for  the  normalized  lift  build  up  on 
the  tail  due  to  a  step  change  in  wing  angle  of  attack 
resulting  from  wing  motion  (Hobbs -type  motion  function) 
will  be  obtained  from  equation  (B 19;  with  Mg  -  °°. 

Again  the  functional  variable  6(oo,  m,  ARw,  ARt,  s*,  s£,  s), 
is  replaced  by  ^(s)  where  now  the  subscript,  M,  desig¬ 
nates  wing  motion.  The  result  is 

6M(s)  "  /o  w(°°»  °»  AW’  ST*  ST»  <5>^(s-<Od<5  (B21) 

where  the  downwash  due  to  wing  motion  is  assumed  to  origi¬ 
nate  from  the  trailing  edge  of  the  wing  mean  aerodynamic 
chord  (see  Figure  B6) . 
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UNSTEADY  AERODYNAMICS  OF  SWEPT  PLANFORMS 

DUE  TO  GUST  AND  DOWNWASH 
(Indicial  Generalized  Force  Functions) 


I.  Calculation  of  the  Wing  and  Tail  Indie la 1  Generalized 
Force  Functions 


A.  Gust  on  an  Arbitrary  Planform 


Consider  a  planform  traveling  at  a  constant  horizon¬ 
tal  velocity,  V,  and  encountering  a  stationary  vertical 
gust,wG  .  For  the  representation  shown  in  Figure  Cl 

(assuming  symmetrical  conditions  half  the  planform  Is  shown) 
the  planform  apex  enters  the  gust  first,  and  the  subsequent 
spanwise  rate  of  immersion  into  the  gust  depends  upon  the 
forward  velocity,  V,  and  the  leading  edge  sweep  angle,  A^g. 

The  planform  is  divided  into  a  number  of  rectangular  panels, 
and  it  is  assumed  that  the  lift  build  up  on  each  of  the  panels 
is  defined  by  the  Kussner-type  function  evaluated  for  the 
aspect  ratio  of  the  entire  planform  in  question.  The  in¬ 
cremental  lift  build  up  time  history  on  the  ith  panel  due 
to  a  step  change  in  gust  velocity  is  given  by  the  follow¬ 
ing  basic  expression  in  terms  of  non-dimensional  time  in 
planform  reference  semichords  traveled. 

yn-n)  -  *  <g  V2)[Ay(y1)2b(y1)c/a(y1)]^G 


for  (f(y^)  £  q  £  00 


(Cl) 


^(yvr\)  -  0, 


for  0  £  n  <  6(yj) 


where  the  transport  time,  ?  (y^),  is  related  to  the  plan- 
form  sweep  through  the  following  expression: 


Z(yt) 


tan  A 


LE 


bR 


®(y1> 


(C2) 
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and  the  gravitational  constant,  g,  is  introduced  because 
the  air  density,  p,  is  defined  in  this  report  in  terms  of 
pounds  mass  rather  than  slugs.  For  a  unit  gust,  w  ■  1  foot 

per  second,  equation  (Cl)  becomes 


^(yi»n)  -  |  V[Ay(yi)b(yi)cia(yi)]^(q-6(yi>)  (C3) 

where  the  same  conditions  still  hold  as  in  equation  (Cl) . 

At  the  time  the  gust  reaches  the  panel  leading  edge, 
lift  build  up  commences  according  to  the  Kussner-type 
function,  and  associated  with  this  lift  build  up  there 
exists  a  pitching  moment  about  the  reference  axis .  Until 
the  airfoil  is  completely  enveloped  by  the  gust,  the  local 
center  of  pressure  is  located  at  a  point  forward  of  the 
steady  state  value.  For  the  present  analysis,  however, 
it  is  assumed  to  be  situated  at  the  steady-state  location 
immediately  upon  gust  contact.  The  pitching  moment  about 
the  reference  axis  at  any  time  is  therefore  obtained  by 
multiplying  the  incremental  lift  given  by  equation  (C3) 
by  the  moment  arm  -b(y^)[c(yi)  -  a(yi)]. 

-  -b(yi)[c(yi)  -  a(yi)]^(yi,Ti)  (C4) 

th 

Finally,  the  generalized  force  on  the  i  panel  for 

the  m**1  mode  due  to  a  unit  gust  is  evaluated  from  the 
expression 


%(y1>Ti)  -  -£h  <yi>Vyi,Tl)  +  £a  (yi)  Vyi>r|)  (c5) 

m  m 


where  f^  (y^)  and  f  (y^)  are  the  vertical  and  angular 
m  nm  ^ 

displacements,  respectively,  in  the  nr  mode  at  the 


reference  axis  of  the  iC  panel.  The  minus  sign  in  the 
first  term  is  due  to  the  assumed  opposite  positive  direc¬ 
tions  of  the  lift  and  vertical  displacement  vectors.  The 
index,  r,  in  the  generalized  force  symbol  is  used  to  de¬ 
note  the  origin  of  each  of  the  various  generalized  force 
components  as  follows: 


r  =  1  Motion  of  the  wing 

r  =  2  Motion  of  the  tail 
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r  ■  3  Downwash  on  the  tail  due  to  motion  of 

the  wing. 

r  *  4  Downwash  on  the  tail  due  to  gust  on 

the  wing 

r  =  5  Gust  on  the  tail 

r  *  6  Gust  on  the  wing 

Note  that  equation  (C5)  as  written  applies  only  to  r  *  5 
and  r  *  6.  For  example,  for  the  evaluation  of  the  down- 
wash  contributions  to  the  generalized  forces,  the  Kussner 
type  function,  ^(tj),  in  equation  (C3)  should  be  replaced 
by  the  Hobbs- type  function,  <5(tj)  .  Note  also  that  when  r 
is  either  2,  3,  4  or  5  the  airfoil  parameters  of  the 
generalized  force  equation  refer  to  the  horizontal  tail, 
and  when  r  is  either  1  or  6,  the  section  parameters  are 
those  of  the  wing.  Combining  equations  (C3)  to  (C5)  and 
rearranging  terms,  the  generalized  force  becomes 


qrm(yi’T1)  -  -  |  V  Pm(yi}  “)^(n-6(yi))  (C6) 

where  the  constant  PmCy^00)  is  defined  by 

pm(yi,0°)  =  [Ay(yi)Myi)c.g  (y^  fh  <yi> 

a  m 


+b(yi)[c(yi)  -  a(y1)]fa  (yt> 


(C7) 


The  total  generalized  force  is  now  obtained  by 
summing  over  the  panels. 

Q™(T,)  ’  1liqrm(>'i-T|)  <C8) 

The  final  form  of  the  total  planforra  generalized  force  for 
a  unit  vertical  gust,  after  equations  (C6)  and  (C8)  are 
combined,  may  be  given  in  terms  of  a  parameter  proportional 
to  the  steady-state  generalized  force  and  the  indicial 
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generalized  force  function. 


Qrm(n)  -  -  |  v  Pm(”)#'m(<l) 


(C9) 


where  the  indicial  generalized  force  function  for  the  m*" 
mode  is  defined  in  terms  of  the  basic  data  and  the 
Kiissner-type  function  as  follows: 


i=l  ra  1  1 


P  (®) 
mv  ' 


(CIO) 


I 

J' 

J 

I 

I 

I 

I 

I 

I 

I 

1 

1 


and  the  constant  P  (°°)  is  defined  in  terms  of  the  basic 
data.  m 


Pm(">  ’  <cll> 


The  specific  steps  necessary  for  evaluating  and 
approximating  the  wing  and  tail  indicial  generalized 
force  functions  are  indicated  in  the  following  sections. 


B.  Gust  on  the  Wing 


In  applying  the  above  expressions  to  a  specific  wing- 
tail  combination,  the  general  non-dimensional  time,  tj, 
which  applies  to  an  arbitrary  planform  considered  alone, 
is  replaced  by  the  non-dimensional  time,  s,  which  is 
defined  in  terms  of  an  airplane  reference  distance,  b_. 
Following  the  procedure  discussed  in  Appendix  B,  K 
and  making  use  of  equations  (B6)  and  (B8),  the  wing 
indicial  generalized  force  function  referenced  to  zero 
time  at  the  wing  apex  is  obtained  from  equation  (CIO) . 


s  tpm(y?>oo)^w(s“6(yi))] 


_w 


(C12) 


where  the  wing  Kiissner-type  function  ^y(s)  is  the  exponen¬ 
tial  approximation  given  by  equation  (B7),  and  the 
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letter  W  used  as  a  subscript  or  superscript  denotes  that 
the  planform  parameters  are  those  of  the  wing.  The  shape 
of  the  indicial  generalized  force  function  as  calculated 
from  equation  (C12)  depends  primarily  upon  the  particular 

mode  for  which  ( s)  is  calculated.  In  order  to  obtain  a 

m 

physical  picture  of  these  functions,  Figure  C2  has  been 
prepared  showing  some  possible  shapes. 


Analytical  approximations  to  the  type  of  functions 
represented  by  Figure  C2  are  assumed  to  be  made  up  of  a 
series  of  exponential  and  damped  trigonometric  terms  as 
follows- 


if^0) 


■ct. 


W  4 
mk" 


k=l 


+  [a"3(l-e"a"3Scos  <&.)  -  l&e'^sin  «£,.] 


U  «■*  k 

m3 


m3 


m3' 


(C13) 


The  constant,  ^W(0) ,  is  introduced  to  facilitate  considera- 

tion  of  the  quasi-steady  case,  for  which  it  is  set  equal  to 
one  for  an  instantaneous  lift  build  up  and  a^  and 

b^  are  set  equal  to  zero.  In  general,  the  constants 

appearing  in  equation  (C13)  are  not  all  zero  and  their 
evaluation  will  be  carried  out  according  to  the  technique 
known  as  Prony's  method  which  is  discussed  in  a  later 
section  of  this  appendix. 

C.  Gust  on  the  Tail 

The  tail  indicial  generalized  force  function  referred 
to  zero  time  at  the  tail  apex  is  calculated  from  an 
expression  similar  to  equation  (C12),  where  the  letter  W 
is  replaced  by  T  to  represent  the  tail. 


(C14) 
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This  expression,  like  equation  (C12),  is  derived  from 
equation  (CIO)  by  introducing  a  change  in  the  non- 
dimensional  time  with  the  aid  of  equations  (B6)  and  (Bll) . 
The  tail  Kussner-type  function,  ^(s),  is  given  by 

equation  (BIO) ,  and  the  other  quantities  in  equation  (C14) 
are  similar  to  equations  (C7)  and  (Cll) . 


If  a  flexible  tail  is  considered,  the  indicial 
generalized  force  functions  for  the  various  modes  will 
vary  considerably  as  shown  for  the  wing  generalized  forces 
in  Figure  C2.  Under  the  assumption  of  rigid  stabilizer, 
however,  tail  indicial  generalized  force  functions  become 
numerically  very  close,  since  the  relative  tail  deflections 
in  these  modes  are  themselves  very  similar.  Therefore,  it 
is  necessary  to  calculate  the  tail  indicial  generalized 
force  function  for  the  first  mode  only,  and  it  can  then  be 
assumed  that  those  for  the  remaining  modes  are  equal  to 
the  curve  calculated  for  the  first  mode.  Under  this  assump¬ 
tion  the  analytical  approximation  to  the  tail  indicial 
generalized  force  function  is  given  by 

£<•>  -  M°>  +  J.1a^k(1*e 


+  ^(l-e  m3  008  “W3S> 


A  "am3s  T  t 
b_,e  sin  aTjS] 


Jm3' 


for  m  -  1 


f(C15) 


♦I<s>  ■  ^i(s)  for  2  £  m  £  N 


Again,  ^(0)  is  introduced  to  permit  reduction  of  the  equa¬ 
tions  to  the  case  of  instantaneous  lift  build  up.  Then 
^T(0)  is  set  equal  to  one  and  the  other  constant  coeffi¬ 
cients  are  zero.  In  general,  the  initial  values  of  the 
Indicial  generalized  force  functions  are  zero,  and  the 
other  constants  will  be  calculated  by  Prony’s  method. 
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D.  Downwash  on  the  Tail  Due  to  Wing  Motion  and  Gust 
on  the  Wing 

The  tail  indicial  generalized  force  function  due  to 
downwash  caused  by  wing  motion  will  be  givejj  by  an  expres¬ 
sion  similar  to  equation  (C14),  where  the  Kussner-type 
function,  ^T(s)>  is  replaced  by  the  Hobbs-type  motion 

function  given  by  equation  (B21) . 


6”  <•)  - 


[Pm(yi'  ">  6M(s'S(yi)) 


(C16) 


The  analytical  approximation  is  given  by 


*2  <»>  -  6m<°)  +Jx  *2k  s> 


+  [am3^1-e  ^  008  ^s> 


m3 


KM  f  . 
-  bm3e 


sin  o^s] 


for  m«l 


( s )  -  «"  (s),  2  <  m  <  N 


(C17) 


In  the  case  of  downwash  due  to  gust,  the  necessary  ex¬ 
pressions  are 


«£<<o 


i*l  P“(yi’  ">  6G<S‘3(yi)) 


P>> 


(C18) 
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where  6^( s)  is  given  by  equation  (B20),  and 


6£(s)  -  6g(0)  +  (l-e'^8) 


G  G 

,  ,  G  /,  "“Vs8  G  .G  “aln38  .  G  , 

+  [am3(l-e  coe  ffi^s)  -  b^e  sin  o^s] 


for  m=l 

6G(s)  =  6?(s)  for  2  £  m  <  N 
m  1 


(C19) 


r  Figure  C3  shows  a  typical  calculated  and  approximated 
6^(s)  function.  Because  the  calculated  indicial  generalized 

force  function  has  comparatively  small  amplitudes  initially, 
it  is  extremely  difficult  to  obtain  an  approximation  based 
on  Prony's  method  which  is  accurate  throughout  the  complete 
time  history.  However,  the  introduction  of  an  incremental 
time  delay  in  the  approximation  (As^  as  indicated  in  Figure  C3) 

improves  the  representation  of  the  Hobbs-type  functions  con¬ 
siderably. 
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II.  Application  of  Prony's  Method  of  Analytical 

Approximation  of  the  Indlcial  Generalized  Force 
Functions 


As  stated  previously  the  approximations  to  the  indlcial 
generalized  force  functions  will  be  of  the  form  given  by 
equations  (C13),  (C15),  (C17),  and  (C19).  The  method  used 
in  evaluating  the  constants  in  these  expressions  is  a 
special  case  of  Prony's  method  which  is  discussed  in 
Section  9.4  of  Reference  16  .  The  following  presentation 
is  more  detailed  than  that  given  in  Reference  16  in 
order  to  summarize  the  specific  sets  of  equations  that 
must  be  programmed  on  the  digital  computer.  Some  typical 
shapes  of  the  curves  that  must  be  approximated  are  indi¬ 
cated  in  Figures  C2  and  C3. 

Prony's  method  approximates  a  curve  which  approaches 
zero  as  tne  independent  variable  approaches  Infinity; 
therefore,  the  following  new  function  is  defined. 


£<•>  - 


for  gust  effects 


or 


(8) 


for  downwash  effects 


(C20) 


where  the  boundary  conditions  associated  with  the  above 
expression  are 


f*(0)  -1.0 

£<->  -  o- 


(C21) 


A  three  term  approximation  to  either  one  of  equations 
(C20)  which  is  equivalent  to  the  form  of  equations  (C13), 
(C15),  (C17),  and  (C19)  is  as  follows. 

aA  8  aA  8  aA  8 

f*(8)  -  F(s)  -  Axe  1  +  A2e  2  +  A3e  3  (C22) 


where  the  subscript  m  is  deleted  for  convenience, 
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and  the  coefficients  and  exponents  in  the  above  expres¬ 
sion  may  be  either  all  real  or  contain  a  complex  conjugate 
pair.  An  equivalent  form  of  equation  (C22)  is 


F(s)  -  AXA®  +  A2A|  +  A3a| 


where 


(C23) 


(C24) 


To  perform  the  actual  numerical  approximation,  values 
of  the  curve  f(s)  must  be  specified  for  J  equally  spaced 
values  of  s  where  J  >  6,  since  there  are  six  constants  to 
be  evaluated  in  equation  (C23) .  By  introducing  the 
transformation, 

Sj  -  jAs,  j  -  0,  1,  2,  ...»  J-l  (C25) 


or 


(C26) 


equations  (C23)  and  (C24)  are  now  rewritten  as 


F(j)  s  Fj  -  Aj\{  +  A2X^  +  A3\J 


(C27) 


where 


a.  As 

.  _  .As  o  \ 

\  \  e 


(C28) 


Equation  (C27)  is  now  expanded  to  form  a  set  of  equations 
for  the  solution  of  the  coefficients,  A^,  assuming  for 

the  moment  that  the  X^  are  known. 
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Fq  m  +  A2  +  A3 

F1  "  A1X1  +  A2X2  +  A3X3 

•  •  •  • 

Fj  ■  ■**  ^2X2  A3X3 

•  •  •  • 

•  •  •  • 

•  7-1  0  T—  1  •  T_1 

FJ-1  ’  A1X1  +  A2X2  +  A3X3 


Y  (C29) 


This  set  of  J  linear  equations  in  the  unknowns,  Ap  A2, 
and  A^  may  be  solved  exactly  if  J  -  3,  or  approximately 
by  the  method  of  least  squares  if  J  >  3. 

Since  the  X^  are  actually  not  known,  it  is  now 
assumed  that  Xp  X2,  and  X3  are  the  three  roots  of  the 
equation, 


(X-Xx)  <X-X2)  (X-X3)  -  0 


(C30) 


where  the  coefficients  Ep  £2*  And  E3  in  the  above  ex- 

?ression  are  real  and  are  determined  accerding  to  the 
ollowing  procedure.  The  first  of  equations  (C29)  is 
multiplied  by  e3,  the  second  by  tha  third  by  ep  and 

the  fourth  by  -I. 


e3F0  "  e3^Al  +  A3) 

£2^^  ■  e2^AlXl  A2^2  +  A3X3) 

e^F,  -  e,(A,X?  +  A,X?  +  A,X?) 


(C31) 


I 


ox 


I 

I 

I 

I 

J 

J 

J 

J 

J 

I 

I 

J 

I 

I 

I 

I 

I 

1 


Since  the  X^  satisfy  equation  (C30),  the  result  of  add¬ 
ing  the  expressions  in  equation  (C31)  will  be  of  the  form. 


F3"elF2“e2Fl“e3F0  "  0 


(C32) 


A  set  of  J-3  similar  equations  Is  obtained  by  starting 
successively  with  the  second,  third,  ...,  ana  (J-3)cl;i 
expression  of  equation  (C29) .  Thus,  equations  (C29)  and 
(C30)  imply  the  J-3  linear  equations. 


F3  ■  eiF2  t  £2Fi  e3fq 


“f*  £2^2  £3Fi 


FJ  *  elFJ-l  +  e2PJ-2  +  £3Fj-3 


FJ-1  "  elFJ-2  e2FJ-3  ^  e2Fj-4 


Since  the  ordinates  Fj  are  known,  equations  (C33)  can  be 
solved  directly  for  if  J  -  6,  or  approximately  by  the 
method  of  least  squares  if  J  >  6. 

Subsequent  to  evaluating  the  e^’s,  the  X^  are  found 

as  roots  of  equation  (C30) .  The  roots  may  be  either  all 
real  or  contain  a  complex  conjugate  pair.  With  the  X^ 

known,  the  A^  are  now  determined  by  applying  the  least 

squares  technique  to  the  set  of  linear  equations  (C29). 

As  stated  previously  the  type  of  curve  to  be  approxi¬ 
mated  is  such  that  F(s)-**0  as  s~*>«°;  therefore,  the  real 
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parts  of  the  exponents,  o^,  should  be  negative.  If  the 
are  all  real, 

\  -  14*3  (C34) 

then  the  approximation  will  be  of  the  form, 

3  -a  .  s 

V>  -  kr.1  amke  <c35> 


where  the  modal  Index  m  Is  now  reintroduced  to  Indicate 
the  dependence  of  the  Indlclal  generalized  force  function 
on  the  mode.  The  coefficients  a^  are  given  by 


amk  *  (C36) 

where  the  A*,  are  obtained  from  equations  (C29)  for  the 
fch  * 

m  mode.  The  exponents  are  given  by 

■  jj  /nX^,  OO^Cl.O  (C37) 


where  the  X^  are  given  by  equation  (C34). 

The  exponential  approximation  of  the  Indlclal 
generalized  force  function,  as  given  by  equation  (C35) 
will  approximate  a  curve  of  the  type  (a)  In  Figure  C2. 

In  the  case  of  curves  (b)  and  (c).  It  appears  that  the 
approximation  will  Include  damped  trigonometric  terms. 

For  this  condition,  the  solution  of  equation  (C30)  will 

consist  of  complex  roots,  or  more  specifically,  two 

roots  must  be  a  complex  conjugate  pair  (since  the  are  real) 
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and  one  of  the  roots  must  be  real.  That  is 

X1  “  Xlp 

and 

X3  “  X3p  +  iX3t 
X2  "  X3p  "  iX3t 


(C38) 

(C39) 


These  roots  are  now  converted  Into  the  exponents  in 
equation  (C22)  with  the  aid  of  equation  (C28) . 


KS  M  xlp 


Y 


“A2  "  TS  A‘<x3p  _ix3 «,) 


'  ^t)-i  tan'1 


(C40) 


Es  to<x3p+lx3p 
Sjt%A>(X.3p  +  »jt>  +  i  tan*1 


/ 


The  corresponding  coefficients  in  equation  (C22),  which 
are  obtained  as  solutions  to  equation  (C29),  are  of  the 
form, 


(C41) 


A3  "  A3p  +  1  A3t 
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Equation  (C22)  therefore  becomes 

“A  •  %o'i0A3l)S 

F(s)  -  Alpe  Alp  +  (A3  -  iA3t)e  3P 


(a.  -iaA  )s 
+  (A3p+  ^31 )e  3P  3t 


(C42) 


where 


a.  ■  xr  t.n 
Alp  As  lp 


aA3p  "  hi  /n*X3p  +  X3i} 
1  -1  *3i 

X  ^ 


y  (C43) 


J 


Finally,  rearranging  terms  in  equation  (C42)  yields  the 
final  form  of  the  approximation  used  for  indicial  gener¬ 
alized  force  functions  of  the  type  illustrated  by 
curves  (b)  and  (c)  of  Figure  C2. 

Fm(s)  “  amle  a“l8 


+  e 


-Vi8 


tam3  008(0)^8)  +  b^sinfai^s) )  (C 44) 


where 


“ml  ’  Alp 
am3  ‘ 


/ 


(C45) 
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a“X  ‘  ~\p 

am3  '  “°A3p 


(C45) 


) 


An  alternate  approach  to  solving  for  combined  ex¬ 
ponential  damped- trigonometric  approximations  Is  as 
follows.  The  X^  are  obtained  as  before,  and,  after 

having  established  that  X2  and  X^  are  complex  conjugates, 

use  Is  made  of  the  representation. 


la.  )As 
A3i 


(a*  +  ia*  )As 
e  A3p  A3i 


\  (C46) 


i 

1 

J 


Equation  (C42)  Is  now  rewritten  In  the  form, 

jaA  As  jot A  As 

Alp  A 


F,  -  A,^e  +  2A0„e  cos  jaA  As 


j  _  Alp*  T  “3p' 


3i 


ja.  As 
*3P 


sin  ja.  As 
A3i 


(C47) 


This  result  is  now  used  to  obtain  a  system  of  equations 
similar  to  equation  (C29). 
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F0  -  Bx  +  B2 


F1  “  Blvl  +  B2v2  +  B3v3 


Fj  -  +  B2vJ  +  B3v3 


•  •  T-1  •  T-l  •  T«1 

FJ-1  “  Blvl  +  B2V2  +  B3v3 


;  (C48) 


where 


B-, 


Bs 


B, 


IP 

2A3p 

-2A 


and 


3i 


ja.  As 


1  Alp 

;1  “  e  ^ 

JaA,.AS 

As 


■  V*  A  «-»« 

e  cos  ja. 

a3i 

sin  1aA  As 
A3i 


j<*A,  As 


l  (C49) 


(C50) 


The  solution  of  equations  (C48)  for  the  unknowns,  B^, 

where  the  Vj|  are  given  by  equation  (C50)  and  a*  ,  aA  , 
K  Alp  A3p 

*  F0"  Blvl  +  B2v2  +  B3v3  where  V1  “  v2  "  1  and  v3  ■  0 
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and  aA  are  given  by  equation  (C43) ,  constitutes  the 
a3i 

second  step  in  this  alternate  approximation  procedure. 
The  final  form  of  the  approximation  for  each  mode. 


F  (s)  -  a  .e 
mN  '  ml 


+  e 


-v>8 


t%3  008(^8)  +  sin  (<Bm3s)]  (C51) 


is  now  developed  by  means  of  the  relationships. 


aml  ~ 

B1 

am3  " 

B2 

bm3  “ 

B3 

<Vl  “ 

"aA 

Alp 

am3  " 

"°A 

A3p 

co  o  »  a. 
m3  A3l 

f  (C52) 


j 


At  this  point  it  is  noted  that  the  parameter  a.  in 

3i 

equation  (C47)  may  be  replaced  by  the  somewhat  simpler 
expression, 

-!  I X3t  I  (C53) 


3i 


.L 

As 


tan 


'*Tr 


with  no  loss  in  generality  for  the  present  application. 
From  equation  (C28),  it  is  seen  that  for  a  complex  root 
of  equation  (C39), 
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(aA  +  ia.  )As 
e  A3p  A3i 


a. 


k3p 


As 


[cos(aA  As) 
a3i 


+  i  sin(aA  As)] 
A3t 


(C54) 


For  the  adequate  numerical  representation  of  a  curve 
having  a  frequency  component,  aA  ,  the  points  F,  should 

A  q.  j 


be  chosen  with  a  spacing  As  of  less  than  approximately 
1.5/aA  .  Under  these  circumstances, 

A3t 


aA  As 
A3l 


< 


TT 

7 


(C55) 


and  therefore,  the  problem  of  choosing  between  multiple 
arguments  in  equation  (C43)  will  not  occur. 

Consider  now  the  sign  of  aA  .  If  equation  (C53) 

a3i 

rather  than  equation  (C43)  is  used  in  the  evaluation  of 
the  by  equations  (C50),  the  sign  of  may  be  reversed. 

However,  the  use  of  -(v^)  rather  than  the  actual  in 

equations  (C48)  will  merely  result  in  a  corresponding 
change  in  sign  of  the  unknown  B^.  The  net  result,  there¬ 
fore,  is  that  the  product  remains  unaffected,  and 

hence,  the  Fj  given  by  equation  (C47),  as  well  as  the 

corresponding  final  form  of  the  approximation  given  by 
equation  (C53) ,  will  remain  correct  in  spite  of  a  change 
in  sign  of  aA 

A3i  . 

Finally,  a  general  form  of  the  approximation  to  the 
gust  functions  is  established  by  combining  equations  (C20), 
(C35) ,  and  (C44)  as  follows. 
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^m<s> 


”amls 


)  +  a, 


m2 


-ct  0s 
(1-e  m2  ) 


4-  a 


m3 


.  i  ”am3S  , 
[1-e  cos(a> 


m3 


8)]-b, 


”am3s 


m3' 


si-nC^s)  (C56) 


The  corresponding  approximations  for  the  downwash  indicial 
generalized  force  functions  will  also  be  of  this  form,  with 
the  variable  ^m(s)  replaced  by  6m(s) .  For  an  approximation 


m 


to  consist  of  three  exponential  terms. 


m3 


CD 


m3 


(C57) 


In  the  case  of  one  exponential  and  one  damped  trigonometric 
term,  the  second  exponential  term  is  deleted,  that  is 


am2 


-  0 


(C58) 


In  conclusion  it  is  noted  that,  while  Prony's  method 
will  determine  a^'s  whose  sum  is  approximately  1,  the 

a^'s  should  be  modified  so  that  their  sum  is  exactly  1. 

These  corrections  will  yield  indicial  generalized  force 
functions  which  are  consistent  with  steady-state  values 
of  unity. 
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^3 


GENERALIZED  FORCE  DUE  TO  GUST  AND 
ASSOCIATED  DOWNWASH 


I .  Gust  on  the  Wing 

Planform  (wing  or  tail.)  generalized  forces  due  to  gust, 
which  are  part  of  the  inputs  to  the  equations  of  motion,  are 
evaluated  through  Duhamel's  integral  in  conjunction  with 
indicia 1  generalized  force  functions  developed  in  Appendix  C 
and  the  forcing  function>  x(s),  which  represents  an  arbi¬ 
trary  gust.  The  development  of  the  generalized  force  begins 
by  considering  the  incremental  lift  and  moment  existing  on 
the  i^1  panel,  and  these  quantities  are  then  converted  to  a 
generalized  force  by  making  use  of  the  mode  shapes.  The 
incremental  lift  on  the  i“*  panel  is  given  in  terms  of  the 
non-dimensional  time  s  as  follows'. 

yyrs)  =  \  |  V2  tAy(yt)  2b(yi)  cA  (yi)]f^/(s-6(yi>)  (Dl) 

In  general,  the  forcing  function,  f^(s-(5  (y^)  ) ,  on  each 

panel  is  given  by  Duhamel's  integral  in  terms  of  the  Kussner- 
type  function,  ^(s),  and  the  change  in  angle  of  attack  due 
to  gust,  x(s)/V, 


f^(s)  =  fo  tx(<0/V]  if  (s-6) do  (D2) 


Note  that  when  x(s)  *  1.0  (unit  step  gust  input)  equation 
(Dl)  reduces  to  equation  (C3) . 

|_l_ 

The  generalized  force  for  the  mc  mode  and  1  panel 
due  to  an  incremental  lift  given  by  equation  (Dl)  is  similar 
to  equation  (C6) . 

^rm(yi»s)  *  ’CV  lp^2pm(yi}0o)  1  (D3> 

where  the  constant  Cy  is  defined  as  follows: 

Cy  -  (V/V2  (VS)  (D4) 
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and  PTn(yi,C0)  is  defined  by  equation  (C7).  The  subscript 

r  is  assigned  the  value  of  5  to  denote  gust  on  the  tail 
and  6  to  denote  gust  on  the  wing. 

The  total  generalized  force  may  now  be  obtained  by 
summing  over  the  panels  and  rearranging  equation  (D3)  as 
follows'. 


Qrm<s>  * 


o  r  /j  *  j  °°\ 


m 


(D5) 


where  the  constant  Pm(°°)  is  given  by  equation  (Cll)  . 

After  substituting  equation  (D2)  into  equation  (D5) ,  the 
generalized  force  may  be  written  as 


()  (s)  =  C„  D  T  (s) 

rm  V  m  m 


(D6) 


where  the  constant,  Dm,  is  given  by 

Dm  =  -(p/V)b£  Pm(co)  (D7) 

The  generalized  force  variation,  ¥m(s),  is  given  by 

Duhamel's  integral  in  terms  of  the  gust  input,x(s),  and  the 
indicial  generalized  force  function,  ^m(s) . 

*m(s)  -  / o  X<c)  (s"6>  d6  (°8) 


where  the  indicial  generalized  force  function  is  given  by 
the  following  basic  expression,  which  is  the  same  as  equation 
(C12)  . 


I  P-Cy^00) 

Vs)  Pm-R-  -^s-g(yi)>]  <D9> 

Note  that  thus  far  in  the  development  the  letter  W  used  as  a 
subscript  or  superscript  to  designate  gust  on  the  wing  has 
been  omitted  since  the  expressions  are  general  and  equally 
applicable  to  both  the  wing  and  the  tail.  It  is  understood 
that  the  expressions  may  be  applied  in  other  parts  of  the 
report  by  the  use  of  the  appropriate  subscript. 
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The  solution  of  equation  (D8) ,  and  hence  also  equation 
(D6),  may  be  obtained  by  combining  it  with  the  analytical 
approximation  for  the  indicial  generalized  force  function 
given  by  equation  (C56).  The  procedure  that  is  followed 
here  is  to  obtain  the  Laplace  transform  of  equations  (D8)  and 
(C56)  as  follows: 


Tm(p)  =  X (P)  P  ^m(P> 


Ip  (p)  =  + 

^mv^'  p 


2 

2 

k=l 


amkamk 

p(p+amk) 


(DIO) 


+ 


[am3  am3  '  bm3mm31  p  +  am3  1  am3+  "41 
p  [p"  +  2am3p  +  (“m3  +  »m3> 1 


(Dll) 


Equation  (DIO)  is  now  redefined  in  terms  of  the  component 
parts  of  the  Laplace  transform  given  in  equation  (Dll) . 

?m(p)  -,f.0  fmk<p>  <D12> 

where  the  component  parts  of  the  Laplace  transform  are  given 

by 


fm0(p)  "  X(P) 


(D13) 


\k<P>  -  t  X(P).  for  k-1,2 


P+Ptak' 


(D14) 


*mk(P>  - 


^  amkamk  ^tnk^mk^  P 


mk  1  mk  mk‘ 


P2  +  2amkp  +  +  “bk> 


X(p) »  (D15) 


for  k-3 
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th 

Finally,  the  generalized  force,  Qrm(s)>  f°r  the  nr 

mode  due  to  the  gust  on  the' wing  is  obtained  from  equation 
(D6)  and  (D12)  through  (D15) .  Also,  the  r  is  assigned  the 
value  of  6  and  the  letter  W  is  used  as  superscript  to  indi¬ 
cate  gust  on  the  wing. 


Q6m<s>  -  cw  d"  *”(0 


V  m  m' 


The  generalized  force  variation,  ¥^(s),  is  evaluated  by 

taking  the  inverse  Laplace  transform  of  equations  (D12) 
through  (D15) . 


iW(s)  =2  ¥W 


m 


k=0 


mk 


(s) 


(D16) 


(D17) 


The  component  parts  of  the  generalized  force  variation, 

^mk ( s ) 1  are  8^ven  hy  the  solutions  to  the  following  set  of 

differential  equations,  which  are  obtained  by  taking  the 
inverse  of  equations  (D13)  to  (D15) . 


*>)  =  ^w(0)  X(s-sw) 


(D18) 


*S*(s)  +  i  *l(s>  "  X  (s-s„),  for  k-1,2 


mk  mk 


mk  mk 


(D19) 


l  X(S-SW)  +  a«k  [(a«k)2  +  (i)2]  X(s-sw) 


for  k=3 


(D20) 


J 


1 


It  is  noted  that  the  transport  time  delay  introduced 
here  accounts  for  the  time  the  gust  requires  to  travel  from 
an  arbitrary  zero  reference  point  located  ahead  of  the  wing 
to  the  wing  apex  as  defined  by  the  intersection  of  the 
airplane  centerline  and  the  wing  leading  edge. 


1 
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II.  Gust  on  the  Tail 


The  development  of  the  generalized  forces  for  gust  on 
the  tail  will  be  similar  to  that  for  gust  on  the  wing.  In 
fact,  the  final  expressions  for  the  tail  generalized  force 
will  be  the  same  as  equations  (D16)  to  (D20)  except  that  the 
subscript  r  now  has  the  value  of  5  and  the  letter  W  is  re¬ 
placed  by  T  to  represent  gust  on  the  tail. 

Since  the  indicial  generalized  force  function  for  gust 
on  the  tail  is  assumed  independent  of  the  mode  as  indicated 
by  equation  (C15),  the  generalized  force  for  gust  on  the  tail 
reduces  to 

<WS>  “  CV  Dm  <D21> 

T 

where  the  tail  generalized  force  variation,  ^(s),  is 

obtained  by  solving  equations  similar  to  (D17)  through  (D20) 
for  the  first  mode  only  (m=l) . 
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III.  Downwash  on  the  Tail  Due  to  Gust  on  the  Wing 


The  expressions  for  the  generalized  forces  on  the  tail 
due  to  downwash  as  a  result  of  gust  on  the  wing  are  similar 

in  form  to  equation  (D16),  with  the  variable  'f^(s)  being  re- 
placed  by  A^(s),  the  downwash  slope  -  de/da  being  intro¬ 
duced  on  the  right  side  to  account  for  the  strength  and 
direction  of  downwash,  and  the  index  r  in  Qrm(s)  being 

assigned  the  value  of  4. 


<WS>  -  CV  Dm 


(D22) 


(D23) 


and  that,  as  in  the  previous  section,  the  generalized  force 
variation  due  to  downwash  on  the  tail  caused  by  gust  on  the 

Q 

wing,  Am(s),  is  independent  of  the  mode  (see  equation  (Cl9.)) 

Q 

The  expression  for  calculating  Am(s)  (for  m=l)  will 
be  similar  to  equations  (D17)  to  (D20) . 


&G(s) 


3 

;  2 

k=0 


A«k(s) 


(D24) 


where 

A10<s>  -  -  IM|6G<°>  X(s-sg)  (D25) 

AL<S>  +  *5k  A?k<s>  -  -  |!l|  4  <4  x(s-sg> 


for  k*l,2 


(D26) 
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\lk(s)  +  2aJk  A^k(s)  +  UaJk)2  +  (a^k)2]  A^k(s)  = 


!  d  e 

3a 


[a 


IP 


alk  "blk  a^lk]  X(s-sG)  +  a 


G 

lk 


(  (alk)  2  +  ]x(s-^j) 


for  k-3 


(D27) 


The  time  delay  sG  is  introduced  to  account  for  the 

elapsed  time  between  the  passing  of  the  gust  over  the  zero 
reference  point  and  the  approximated  start  of  tail  lift 
build  up  resulting  from  the  wing  downwash  due  to  that  gust 
More  specifically,  sG>  is  the  sum  of  the  quantity  As"  fho^n 

in  Figure  C3  and  the  distance  in  semichords  between  the  zero 
reference  point  and  the  longitudinal  location  of  the  point 

aerodynamic  ^hordf  "ing  leadlng  edge  and  itS  mean 


GRUMMAN  AIRCRAFT  ENGINEERING  CORPORATION 


1  I  \  I 


DATE 


OCTOBER  1963 


APPENDIX  E 


GENERALIZED  FORCE  DUE  TO  MOTION 
AND  ASSOCIATED  DOWNWASH 


I .  Wing  Motion 


For  a  two  dimensional  airfoil  section  in  incompressible 
flow,  the  incremental  lift  and  pitching  moment  about  the 
reference  axis  due.  to  oscillatory  airfoil  motion  is  given 
by  Ref erencel5 ,  page  272,  equations  (.5-311)  and  (5-312), 
respectively.  These  expressions  are  reproduced  below  in  a 
slightly  modified  form;  most  of  the  symbols  used  are  de¬ 
fined  in  Figure  Cl. 


i0(t)  =7r|b2  Ch(t)  +  V4(t)  -baa(t)] 
+  27r|vbf0(t) 


(El) 


t0(t)  =  7T  |b2  [ bah(t)  -  VbG-a)A(t)  -  b2 (a2+l/8)a(t)  ] 


+  2tt  |  Vb2(a+%)f0(t) 


(E2) 


where 

f0(t)  =  "  C(x)w0(t) 


(E3) 


C(k)  is  Theodor sen's  function,  and  the  downwash  at  the 
three-quarter  chord,  w^(t),  is  given  by 


w^(t)  =  -  [h(t)+Va(t)+b(^-a)d(t)]  (E4) 

Note  that  the  density, p,  is  defined  in  units  of  pounds-mass 
rather  than  slugs,  and  is  therefore  divided  by  the  gravi¬ 
tational  constant,  g.  The  subscript,  0,  is  used  to  denote 
lift  and  moment  due  to  motion. 
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Since  the  circulatory  contribution  to  the  lift  and 
pitching  moment  given  by  equations  (El)  and  (E2)  is  defined 
in  terms  of  Theoaorsen's  function, C(k) ,  these  equations  are 
applicable  to  arbitrary  motion,  as  is  required  for  an  arbi¬ 
trary  gust,  equation  (E3)  is  redefined  in  terms  of  Duhamel's 
integral  and  the  Wagner-type  function.  In  order  to  carry 
out  this  change,  use  is  made  of  equation  (5-370),  page  285, 
and  equations  (8)  through  (10),  page  815,  of  Reference  15. 
Equation  (E3)  is  now  presented  in  terms  of  Duhamel's  integral 
as  follows: 

^(t)  =  ~f 0w<j/T^  (t_T)dT  (E5) 

Note  that  in  the  subsequent  development,  where  the  independent 
variable  is  changed  from  t  to  s,  the  function  </>(t)  is  re¬ 
placed  by  the  more  familiar  form,  0( s) . 

The  above  equations  for  section  lift  and  pitching  moment 
are  now  utilized  to  describe  the  unsteady  aerodynamic  forces 
due  to  motion  of  the  itJl  panel  of  a  three  dimensional  wing 
in  compressible  flow.  Under  the  following  assumptions, 
the  equations  are  re-written  below  in  terms  of  the  section 
properties  and  airfoil  motion  at  the  spanwise  position  cor¬ 
responding  to  the  center  of  the  ith  panel,  y^. 

(1)  .  The  circulatory  lift  and  moment  on  the  iL  panel 
is  determined  by  the  downwash  at  the  three-quarter  chord 
point  of  that  panel  (equation  (E4)). 

(2)  .  The  magnitude  of  the  circulatory  lift  and  moment 

on  each  panel  is  based  on  the  local  lift  curve  slope  (2m  is 

replaced  by  Cf  (y. )),  and  local  aerodynamic  center (the 

Ot  1 

factor  (a+^)  in  the  last  term  of  equation  (E2)  is  re¬ 
placed  by  (a(yi)-c(yi)) . 

(3)  .  The  rate  of  circulatory  lift  build  up  on  a  three- 
dimensional  airfoil  is  constant  across  the  span,  and  is 
taken  as  the  build  up  on  a  rectangular  airfoil  having 
the  same  aspect  ratio  as  the  actual  airfoil  being  con¬ 
sidered.  The  expression  for  the  circulatory  lift  due 
to  motion  is  that  given  by  the  Wagner- type  function, 
whose  analytical  approximation  is  of  the  form  given  by 
equations  (B15)  and  (B17)  for  the  wing  and  tail,  re¬ 
spectively. 
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(4) .  Non-circulatory  lift  per  unit  span  on  a  given 
panel  is  assumed  to  be  identical  to  the  lift  on  a  two- 
dimensional  section  in  incompressible  flow. 

The  modified  lift  and  moment  equations  are 


Vyi’C)  "  2l  C(yi) 


My.^  [h(yi,t)  +Vc5t(yi,t) 


-b(yi)a(yi)  aCy^t)]  +  2Vf^)(yi,t) 


(E6) 


-  £g  C(yi)b(yi)  j  b(yi)  [a(yi)h(yi> t) 


+  V  taCy^-^  3a(yi}  t)-b(yi)  [a2(yi>  +1/8]  a(yi^] 


-2V  [c(yi)-a(yi>]  f 0(yi>  t) j 

where 

C(yi)  -  Ay(yi)b(yi)c^  (yi) 
f0(yi>t)  =  -/q  0  (t-r)dT 


(E7) 

(E8) 

(E9) 


and  the  downwash  at  the  three-quarter-chord  point  becomes 
w0(yi,t)  =  -  [li(yi,t)+Va(yi,  t)-b(yi)  [a(yi)-%]  a(yi}t>]  ^ElQ^ 

The  development  of  the  present  analysis  is  based  upon 
the  non-dimensional  time,  s,  rather  than  the  real  time,  t. 

The  conversion  of  the  time  variable  is  given  by  the  following 
relationships, 
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s  - 

& 

d 

Ht  = 

/V  s  d 

3? 

d2 

(v  )2  d2 
V  z 

dt^ 

(Ell) 


J 


.  th 


Next,  the  vertical  and  angular  deflections  of  the  i  panel 
are  expressed  as  the  product  of  time  dependent  coordinates 
(normal  modes)  and  space  dependent  coordinates  (mode  shapes) . 


h(yi3s)  =  2  fh  (yi)  4  (s) 
n  n 

! 

<x(yi,s)  -  s  f  (Yl)  en(s) 

n  n  , 


(E12) 


Substitution  of  equations  (Ell)  and  (E12)  into  equations 
(E6)  to  (ElO)  yields  the  expressions  for  the  lift  and  pitch¬ 
ing  moment  on  the  it“  panel  in  terms  of  the  non-dimensional 
time,  s,  and  the  normal  coordinates,  4  n(s) .  The  expression 

for  the  lift  on  the  ifc^  panel  is 


Vyi’s)  “V(2  <4><yi>  ^n(s)+^n(yi)  ^n(s) 

+/q  fi0^yi»6)  (s-c5)dd 
where  the  function  f^( yps)  is 


*  l  +  C<LX> 


(E13) 


(E14) 
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The  coefficients  are  defined  in  terms  of  the  basic  data  as 
follows; 


=  %C(yi)b(yi)  [fh  (yi)-b(yi)a(yi)fQ^(yi)] 

-  %bRC(yi)b(y1)  f^y^ 

“  bRC(yi)  ffh  (yi)-b(yi)  ta(yi>  '^1  fa  (yi)l 
n  n 


C.ta<yi>  ■  bR  C<yi>  fan(yi) 


(E15) 


) 


A.V 

In  a  similar  manner  the  pitching  moment  on  the  i  panel  is 


t0(yt,s)-  Cyp 


l  t  At„ ^  s>  +B?n  tyi) C S  > 1 


+  ^Oft0^yi,fi^  0  (s“6)dc3 


(E16) 


where  the  function  ft0(y£>s)  is 


ft0(yi>s)  -  S  lB^n(y.)  en(s)  +  C=n(yi)  5n(s)]  (E17) 

n 


and  the  coefficients  are  defined  as  follows; 


A^n(yi>  =  h  c(yi)b^yi)Ca<yi>fh  (yi)-b(yi)  [a2(yi)+V8]f  (s)] 

n  ^n1 
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^n(yi)  “  ^bR  c(yi)b2(yi)  My*)-^  fa  (y^ 

n 


(E18) 


Btn<yi>  =  "  b<yi>  [c<yi>  -  a<yi>]  Bin(yi> 
Ctn(yi)  =  "  b<yi>  fc<yi>  ‘  a(yi)  3Cin(yi) 


The  constant  appearing  in  equations  (E13)  and  (E16)  is 
given  by 


CV  =  (V/bR)2(l/g) 


(E19) 


The  generalized  force  in  the  m 
panel  is  now  defined  as 


th 


mode  due  to  the  i 


th 


~  [Vyi,s)fhm(yi)"Vyi,s>fam(yi>]  (E20) 

The  total  generalized  force  may  be  obtained  by  summing  over 
the  panels. 


Qrm<s>  ‘•rm<>'i>s>  <E21) 

Note  that  the  index  r  is  assigned  the  values  of  1  and  2  to 
denote  that  the  generalized  force  occurs  due  to  wing  motion 
and  tail  motion,  respectively. 

The  equations  developed  thus  far  in  this  section  may 
be  used,  if  desired,  to  compute  the  total  generalized  forces 
due  to  motion.  Specifically,  equations  (E13) ,  (E16) ,  (E20) , 
and  (E21)  may  be  combined,  and  then  used  to  evaluate  the 
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contributions  of  both  the  wing  and  the  tail.  A  more  con¬ 
venient  form  of  these  equations,  however,  may  be  derived  by 
interchanging  the  order  of  summation  used  thus  far,  that  is, 
by  summing  over  the  panels  first  and  then  summing  over  the 
modes,  rather  than  vice  versa.  The  wing  (r=l)  total  gen¬ 
eralized  force  may  now  be  expressed  in  terms  of  a  non- 

circulatory  component,  <£^(s),  and  a  circulatory  component, 

*WC(s),  as  follows: 
m 

Qim<«>  -  -cv[<C<s>  +  *mC<s>]  <E22) 


where  the  non-circulatory  contribution  to  the  generalized 
force  is  expressed  in  terms  of  the  first  and  second  deriva¬ 
tives  of  the  normal  modes. 


<tJJN(s>  =  2  [a)5J  en(s)  +  i  (S)] 

m  _  mn  n'  ran  n 

n 

and  the  coefficients  A^  and  of  the  non-circulatory 

mn  mn 

contribution  to  the  generalized  force  are  defined  below. 


(E23) 


P  f tAM  V^>-A?n (yi)fam(yi>  1 


(E24) 


rWN  . 
Bmn 


(E25) 


where  the  coefficient  appearing  within  the  brackets  in 
equations  (E24)  and  (E25)  are  evaluated  from  the  expressions 
given  by  equations  (E15)  and  (E18)  based  upon  the  wing 
properties.  The  contribution  of  the  circulatory  part  to 
the  generalized  force  is  given  by 

*mC(s)  ’  fi  Fm(S)  <s-0)d6  (E26) 

o 

where  the  forcing  function  F^(s)  is  defined  in  terms  of  £  (s) 

and  its)  as  follows:  m  n 

n 
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o 


FW(s)  =  Z  [BWC  i  (s)  +  CWC  i  (s)] 
m v  '  1  mn  v  7  mn  v  J 

n 


and  the  coefficients  in  equation  (E27)  are  given  by 

‘  ■  p  f  <❖  "M>  £«  &?>] 


-  P  ?[CM£h  <J?)  <,<»?>*« 

i  m  tn 


(E27) 


(E28) 

(E29) 


where  once  again  the  coefficients  appearing  on  the  right 
hand  side  of  equations  (E28)  and  (E29)  are  evaluated  using 
equations  (E15)  and  (El8) .  At  this  point,  use  is  made  of 
the  exponential  approximation  to  the  Wagner-type  function, 
given  by  equation  (B15) ,  which  is  reproduced  below. 

K  —  g 

Vs>  =  0W(O)  +  51  bk(1‘e^  )  (E30) 

k—l 


where 

KX  =  K^  (E31) 

The  Laplace  transform  of  equations  (E26)  and  (E30)  are 


C<p>  -  *>>  p  Vp) 

ii  bX 

P  k-1  p(p4^) 


Now,  redefine  the  Laplace  transform  given  by  equation  (E32) 
as  follows: 


(E32) 

(E33) 
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^C(P)  -*Eo<P>  +  *£<P) 


where 


<(p>  -kix  «£k(p) 


(E34) 


(E35) 


and  the  component  parts  of  the  Laplace  transform  in  equations 
(E34)  and  (E35)  are  defined  by  combining  equations  (E32) 
and  (E33) . 


*£o<p>  -  v°>  *>> 


Kw' 


bWftW 

-fV  t  v  PkPk  -W,  v 
$mk^  =  |ftW 


k*=l,  2 


(E36) 

(E37) 


The  circulatory  contribution  to  the  generalized  force  is 
now  obtained  by  taking  the  inverse  Laplace  transform  of 
equations  (E34;  through  (E37) . 


<£C<s>  - 


»>>  + 


*"(s) 
m '  ' 


(E38) 


where 


*s) 

m  ' 


2 

-  2 
k-1 


I 

(E39)  , 


1 

I 

I 

I 
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and  the  component  generalized  forces  on  the  right  side  of 
equations  (E38)  and  (E39)  are  obtained  from  the  following 
set  of  differential  equations, 


*>>  -  VO)  Fm<s>  <E40> 

^  k-1.2  (E4D 
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II.  Tail  Motion  and  Stabilizer  Rotation 


The  generalized  forces  on  an  all-movable  horizontal 
tail  depend  upon  tail  motion  and  stabilizer  rotation.  The 
expressions  for  the  generalized  force  contribution  due  to 
tail  motion  will  be  exactly  the  same  as  those  for  the  wing 
which  were  developed  in  section  I.  Development  of  the 
generalized  force  due  to  stabilizer  rotation  may  be  carried 
through  by  following  the  procedure  in  section  I  with  certain 
modifications . 

Basic  expressions  required  for  evaluating  the  generalized 
force  due  to  stabilizer  rotation  are  the  panel  lift  and  pitch¬ 
ing  moment,  which  in  the  case  of  wing  motion,  are  given  by 
equations  (E6)  and  (E7)  together  with  the  auxiliary  equations 
(E8)  to  (ElO) .  These  equations  may  be  changed  to  apply  to 
stabilizer  rotation  by  replacing  the  vertical  and  angular 
displacements  and  their  associated  derivatives  (h(y^,t), 

a(y±, t) , . . .)with  the  variable  representing  stabilizer  ro¬ 
tation,  7(t).  If  the  horizontal  tail  reference  axis,  which  is 
the  torque  tube,  is  chosen  as  the  stabilizer  pivot  axis,  these 
changes  are  equivalent  to  setting 


h(yi,t)  s  0 

aCy^t)  S  7(t) 


(E42) 


Introduction  of  equations  (E42)  into  equations  (E6)  and  (E7) 
will  provide  the  panel  lift  and  pitching  moment  due  to 
stabilizer  rotation. 


E  V-y(t)  -b (yj)  a (y£)V(t) ] 


+2Vf0(yJ,t) 


(E43) 


Ill 
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t0(yT,t)  =  §gC(yi>b(yi)  |  b(yj)  [v[a(yj)-%]7(t> 


-b(yT)  C  a2 (yT)  4-  1/8  ]y(t)] 


-2V[c(yJ)-a(yJ)] 


f  e(yi»t)] 


(E44) 


where 


C(yJ)  -  Ay(yT)b(yJ)c^  (yT)  (E45) 

ct 

fgCyJ.t)  =  -/owe(yi>T)^T(t-T:)d'r  (£46) 


and  the  downwash  due  to  stabilizer  rotation  becomes 

w0(yT,t)  =  -  [V  7(t)-b(yT)  [  a  (yT) ]  7(t)l  (E47) 


Conversion  of  the  panel  lift  and  pitching  moment  to 
non-dimensional  time  may  be  accomplished  by  introducing 
equations  (Ell)  into  equations  (E43)  to  (E47) .  The  panel 
lift  and  pitching  moment  equations,  then  become,  after  re¬ 
arranging, 


M> s> 


V 


C(y?) 


b(yJ) 


^L-[bR7(s)-b(ypa(yp7(s)] 


+  /ofi0^yi’  C)0T(s-6)d6 


(E48) 
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1  b2 (yT) 

t9(yj.s)  -  Cyp  j  C(yE)  — j-i-  -%J  l'(s) 

-b<y£)ta2(y*)+l/81  y(s)] 

+/ofte^i>  «)4(s-|3>d6|  (E49) 

where 

f^0(y£,s)  =  bRC(yT)[bR7(s)-b(y^)[a(y^)-%]7(s)  ]  (E50) 

and 

£te(yi,s>  ’  -bRG (yj) b(yi) t c (y*) -a (y£> J [ bR-y(s) 

-b(y[)ta(yj)-%]7(s)]  (E51) 

The  generalized  force  due  to  tail  motion  and  stabilizer 
rotation  (r=2)  in  the  mtb  mode  due  to  the  itb  panel  is  now 
defined  as 

"l2m(yI>  s>  "  -  [£0<yi>s>fh  <yM<y^>£a  (yl)! 

ra  m 

-  Ii9(yJ,s)fh  (yj)-te(yj,s)£  <yj) 3  (E52) 

m  m 

where  the  first  bracket  on  the  right  side  of  equation  (E52) 
is  similar  to  equation  (E20),  based  on  the  tail  properties, 
and  the  second  bracket  is  the  contribution  of  the  stabilizer 

rotation.  Hie  variables,  i^(yT, s)  and  t^(yT, s),  are  given  in 

equations  (E13)  and  (E16) . 
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The  total  generalized  force  due  to  tail  motion  and 
stabilizer  rotation  is  obtained  by  summing  over  the  panels. 

02m<s>  -  f  q2m(yj,  s)  (E53) 

Following  the  procedure  of  section  I  the  tail  total  generalized 
force  may  now  be  expressed  in  terms  of  a  non-circulatory  com- 

TVf  mp 

ponent,  (s),  and  a  circulatory  component,  0^(s),  in  a 
manner  similar  to  equation  (E22). 

<WS>  -  +  *“<»>]  (E54) 


The  non-circulatory  contribution  to  the  generalized  force  is 
given  by 

$™(s)  =  2CA™  \°(s)  +  B™  |  (s)  ] 
m  mn  n v  '  mn  n 

n 


+  A™  “<s>  +  ^<s> 


(E55) 


The  circulatory  contribution  to  the  generalized  force  is 


-  «£o<s> +  «£<»> 


m 


(E56) 


where 


* l  <s> 


(E57) 


The  component  parts  of  the  generalized  forces  in  equations 
(E56)  and  (E57)  are 


*>>  -  V0)  Fm(s) 


(E58) 


O  m  m  m  m  m  m 

+  fKk<s>  -  bkPK<s>- k  - 1> 2 


(E59) 
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where  the  input  to  the  above  equation  is 


rn  rrp  0  TY1 

Fi(s)  =  2  [BiL  \  (s)  +  C1*'?  (s)  ] 
mv  '  „  L  mn  bnv  '  mnsnv  '  J 

n 


+  ;<s)  +  C™  7<s)  (E60) 

OO  0 

The  coefficients  of  £n(s),  £n(s),  and  £n(s)  in  equations 

(E55)  and  (E60)  are  given  by  expressions  similar  to  equations 
(E24) ,  (E25),  (E28),  and  (E29)  evaluated  for  the  tail. 


A™  ■  p  -  A?„(y?>fa„<yl)] 


m 


m 


B. 


TN 


,N 


N 


W  ■  p  ?[®lin(yp£h  <*i>  -  »Sn&l>£a 
I  m  m 


TC 


C  -  p  f  »M£h  <*i>  - 

l  m  m 


.TC 


*  p  fcin<y|)£h  <yi» 


(E61) 


(E62) 


(E63) 


(E64) 


Finally,  the  coefficients  of  7(s),  7(s),  and  7(s)  are  given 

by 


A™  -  p  2 

m7  5  i 


C(yJ)b2(y£)[-a(y£)fh  (yj) 

m 


+  b(y*)[a2<y£)  +  1/8]  f^  (yj)]  J  (E65) 

B™  -  P  ^  ?[c(yj)b(yj)[fh  (yf)-b(yj)[a(yj)-%]fa  (y£)]  1  (E66) 
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BTC  = 

tiry 


-pbR  2 


r 


C(yT)b(yT)[a(yJ)-k] 


f 


+  b(y^)[c(y'[)-a(y^)]f  (y *)' 
*  m 


(E67) 


2 

l 


c(yi^lfh  (yl)+b(yl)ic<yl)-a<yl) 

m 


(E68) 


Note  that  equations  (E65)  to  (E68)  may  be  obtained  directly 
from  equations  (E61)  to  (E64)  by  changing  the  mode  index  n 
to  7  on  the  left  hand  side  of  these  equations  and  substi¬ 
tuting  equations  (E15)  and  (E18)  together  with  the  following 
equivalent  of  equation  (E42) . 


fh  <y?)  -  o 

XI 

fa  (y?)  *  i 

n 


(E69) 
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III.  Downwash  on  the  Tall  Due  to  Wing  Motion 

For  a  two  dimensional  wing  executing  plunging  and  pitch¬ 
ing  motion,  the  angle  of  attack  governing  the  circulatory 
lift  on  the  wing,  and  hence  the  strength  of  the  vortices 
shed  from  the  trailing  edge,  is  determined  by  the  planform 
downwash  at  the  three-quarter-chord  point  given  by  equa¬ 
tion  (E4)  and  repeated  below. 

w^(t)  -  -[h(t)  +  Va(t)  +  b(%  -  a)a(t)]  (E70) 


In  calculating  the  downwash  at  the  tail  due  to  wing 
motion,  a  spanwise  station  is  chosen  as  a  reference,  and 
the  downwash  at  this  station  is  assumed  to  apply  across 
the  entire  span  of  the  tail.  If  the  tail  mean  aerodynamic 
chord,  y^a  ,  is  chosen  as  the  reference  spanwise  station 

for  calculating  downwash,  equation  (E70)  when  applied  to 
this  particular  station  becomes 

w"<yL c'V  ’  +  Vaw<yLc>t> 

<E71) 


The  downwash  at  the  tail  due  to  wing  motion,  assuming 
instantaneous  downwash  build  up,  is  given  by 


w6(t) 


de 

3a 


w0^ymac’ 


(E72) 


where  tj^  accounts  for  the  time  delay  between  the  wing 
motion  and  the  assumed  start  of  lift  build  tp  on  the  tail. 

Vortices  due  to  flexible  motion  are  assumed  to  be 
dissipated  by  viscous  effects  before  reaching  the  tail; 
therefore,  only  rigid-body  wing  motion  is  considered  in 
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evaluating  the  downwash  at  the  tail  as  given  by  equation  (E71) . 

Hie  vertical  and  angular  deflections  at  the  wing  refer- 

ence  axis  are  expressed  in  terms  of  the  mode  shapes 

and  normal  coordinates  in  a  manner  similar  to  equations  CE12) , 

where  the  summation  is  carried  out  for  the  rigid-body  modes 

only. 


Vic'8)  -J  fh 

n-1  n 

“w^Lc'8)  -  J. 

n-i  n 


(E73) 


Substitution  of  equations  (E73)  into  equation  (E71)  and 
introduction  of  the  non-dimensional  time,  s,  into  equa¬ 
tions  (E71)  and  (E72)  results  in  the  following  relationships. 


w”(ws> 


V  I 

r —  Z 

k  n-1 


fW  (yT  ) 
h^mac' 


[V?Lc> *»<»> 


w6(s) 


3e 

3a 


\ 


+  bRfan<£.c>?„<8> 

(E74) 

n 

w*  <yic-8-V 

(E75) 

The  lift  build  up  on  the  tail  is  evaluated  based  on 
strip  theory,  where  it  is  assumed  that. .the  lift  buildup  on 
each  tail  panel  is  given  by  the  tail  Kussner-type  function 
with  ^(0)  -  0  (see  equation  (BIO)).  The  lift  on  each 
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panel  Is  given  in  terms  of  Duhamel’s  integral  by  the  follow¬ 
ing  expressions. 

l6 (yj, 8)  -  ^  V2t2C(yJ)]£5(8-«(y^))  (E76) 


where 

c(yj)  -  6y(yJ)b<yl)cs  (yj)  (E77) 

s  w.  (6) 

f6(s)  -  /  -^r-  ^(s-6)d6  (E78) 

The  panel  pitching  moment  about  the  reference  axis  of  the 
*.1. 

icn  panel  is  given  in  terms  of  the  i  panel  lift  as  follows; 


t6(yj;,s)  -  “b<yj)[c(yj)  -  a(y£)]i6(y£,s)  (E79) 


The  i^  panel  generalized  force  due  to  downwash  on 
the  tail  caused  by  wing  motion  (r  -  3)  is  given  by 


q3m<yi>s) 


-U6(yJ,s)fh  (yj)  -  t6(y];,s)f  (y|)]  (E80) 


and  the  total  generalized  force  may  now  be  evaluated  from 
equation  (E80)  by  summing  over  the  panels  according  to  the 
following  relationship, 


<53m<s)  *  l  l3m<yi-s> 


(E81) 


A  more  useful  form  of  the  generalized  force  is  obtained 
by  combining  equations  (E76)  and  (E78)  to  (E81) . 
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Q3m<s> 


“ CV[PbR2 


P J(*)]S[ 
m  ^ 


Pm<“> 


f 6 (S-C5 (yj))  ]  ' 


(E82) 


T  T 

where  the  constants  pm(y^,°°)  and  P  (»)  are  evaluated  using 

equations  (C7)  and  (Cll),  and  f^Cs)  is  given  by  equation  (E78) . 

Now,  combining  equations  (E78)  and  (E82)  and  rearranging  the 
terms  results  in  the  following  expressions. 


-CV^Am<s) 


(E83) 


where  the  constant  D 


m 


is  given  by 


V  °R 


£<-> 


(E84) 


M 

and  the  generalized  force  variation-,  a"(s),  is  given  by 

Duhamel's  integral  in  terms  of  the  downwash  due  to  wing 
motion,  w.  (s) ,  and  the  indicial  generalized  force  function, 

£<.). 

m 


(E85) 


M 

where  the  function  6j^(s)  is  given  by  the  following  basic 
expression. 


6M 

m 


(s) 


i=l  P*<») 
m 


6M(s"6(yi)> 


(E86) 


Note  that  <5m(s)  is  given  by  equation  (B21)  . 

Use  of  the  analytical  approximation  to  equation  (E86)  as 
given  by  equation  (C17)  permits  an  evaluation  of  the  integral 
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given  by  equation  (E85)  along  the  same  lines  as  the  evalu¬ 
ation  of  equation  (D8) .  Therefore,  the  final  form  of  the 
generalized  force  variation  due  to  downwash  on  the  tail 
caused  by  wing  motion  will  be  given  by  an  expression  simi¬ 
lar  to  equation  (D17) . 


&m<s)  "  !oi”k<S>  f°r  m  *  L 
A«(s)  »  a“(s)  2  <  m  <  N 


(E87) 


Note,  however,  that  since  the  tail  is  assumed  rigid,  and 
hence  the  tail  indicial  generalized  force  functions  are 
nearly  independent  of  mode  (as  discussed  in  Appendix  C), 
it  only  is  necessary  to  evaluate  the  generalized  force 
for  the  first  mode.  The  component  parts  of  the  first 
mode  generalized  force  are  evaluated  using  the  following 
set  of  differential  equations. 


A**0(s)  -  6m(°)fm(s_sm)’  for  k  -  0 

+ 


for  k  =  1,  2 


Alt(s)  +  2a“kA«k(s)  +  [(c^k)2  +  (^iu)2)A^k(s) 


(E88) 


(E89) 


^ a l!kalk ” b lk^lk ] ^ s ” ^ ^  +  alkC  ^°?k^2+^a?k^2]FM^s"sM^ 


(E90) 


The  expression  for  the  input  to  equations  (E88)  to  (E90) 
is  given  by 


Fm(s)  - 


l  tB«|n(s)  +c«en(s)] 


(E91) 
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where  the  coefficients  in  equation  (E91)  are  given  by 

Bn  “  '  [fhn(yLc)  ‘  bW<yLc> ' a W(yLc> fan^Lc> 


<5-- 


Se 

"Set 


VfW  (yT  ) 
a.  v-ymac/ 
n 


(E92) 

(E93) 
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STRUCTURAL  LOADS 


I.  Aerodynamic  Contributions  to  the  Wing  and  Tail 
Structural  Loads 


A.  Tail  Loads  Due  to  Tail  Motion,  Stabilizer  Rotation, 
Gust  on  the  Tail,  and  Downwash 


Expressions  are  developed  herein  for  the  aerodynamic 
contribution  to  the  structural  loads  at  the  pivot  axis  of 
an  all-movable  stabilizer  at  the  intersection  of  the  sta¬ 
bilizer  and  the  fuselage  side.  Hie  lift  outboard  of  the 
fuselage  (tail  root  chord)  is  evaluated  by  the  force 
summation  method.  Hie  bending  moment  and  torque  on  the 
torque  tube  are  then  obtained  by  multiplying  this  lift  by 
a  bending  moment  arm,  yQ,  and  a  torque  arm,  Xq,  defined 

in  Figure  FI. 


Sf(s)  -  LT(s) 

mTa(s)  *  y0LT<s) 

tJA(s)  -  x0LT(s) 


> 


J 


(FI) 


where  the  function  L^,(s)  is  the  tail  lift  per  side  based 

on  the  surface  outboard  of  the  fuselage  juncture  with  the 
stabilizer  pivot  axis.  This  lift  consists  of  the  follow¬ 
ing  component  parts:  lift  due  to  gust  on  the  tail,  motion 
of  the  tail,  stabilizer  rotation,  and  downwash  on  the  tail 
due  to  gust  on  the  wing  and  motion  of  the  wing.  In  evalua¬ 
ting  these  components,  use  is  made  of  the  expressions 
developed  in  Appendices  D  and  E.  Note  that  the  generalized 


*  Note  that  the  subscript  1  used  with  these  functions 
represents  the  loads  at  the  tail  root  and  is  not  to  be 
associated  with  any  panel  point  index.  Also,  the  letter  T 
represents  loads  on  the  tail,  and  the  letter  A  designates 
the  aerodynamic  contributions  to  the  structural  loads. 


APPENDIX  F 


forces  in  these  appendices  are  based  upon  a  tail  area  de¬ 
fined  between  the  airplane  centerline  and  tail  tip.  To 
evaluate  the  tail  loads,  the  area  outboard  of  the  fuselage 
is  to  be  utilized.  Therefore,  the  final  expression  for 
the  tail  lift  as  developed  from  the  generalized  forces 
will  be  modified  to  conform  to  this  requirement. 

The  total  tail  lift  is  the  negative  of  the  tail  con¬ 
tributions  to  the  generalized  force  for  the  first  mode, 
that  is 


5 

Lt(s)  *  -  2  Qrl(s) 
1  r=2  ri 


(F2) 


The  lift  on  the  tail  due  to  motion  of  the  tail  and  sta¬ 
bilizer  rotation  (r=2)  is  given  by  equation  (E54) ,  which 
is  reproduced  below  for  m  =  1. 

Q21(s)  -  -  Cv  [4>™(s)  +  (F3) 

where  the  non-circulatory  contribution  to  the  lift  is 
given  by  equation  (E55)  and  the  circulatory  contribution 
is  given  by  equation  (E56) .  The  lift  on  the  tail  due  to 
downwash  as  a  result  of  wing  motion  (r*3)  is  given  by 
equation  (E83) ,  also  reproduced  here  for  m  -  1. 


Q31(s)  -  -  Cv  d“  a“(s)  (F4) 

where  the  constant  d”  is  the  first  mode  value  of  Dm  as 

given  by  equation  (E84) ,  and  the  generalized  force  varia- 

tioni  A^(s),  is  given  by  equations  (E87)  to  (E90) .  The 

lift  on  the  tail  due  to  downwash  as  a  result  of  gust  on 
the  wing  (r-4)  is  given  by  equation  (D22)  for  m  *  1. 

Q4l(s)  ■  Cy  A^(s)  (F5) 

where  D^  is  obtained  from  equations  (P23)  and  (D7)  and 
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A^(s)  is  given  by  equations  (D24)  to  (D27) .  Finally,  the 

lift  on  the  tail  due  to  gust  on  the  tail  (r=5)  is  given 
by  equation  (D21)  for  m  =  1. 

Q5l(s)  -  Cy  D*  TT(s)  (f6) 


Note  that  in  equations  (F4)  to  (F6)  the  following  equalities 
are  true  between  the  constants  D-^. 


(F7) 


Equations  (F2)  to  (F7)  are  now  combined  to  provide  the 
final  expression  for  the  lift  on  the  tail  based  on  an 
area  between  the  airplane  centerline  and  tail  tip. 


ys> 


[•?<•> 


•Ft.)]  - 


d^[a“(s) 


A?(s) 


'i'T(s)] 


(F8) 


Expression  (F9)  may  now  be  modified  to  represent  the 
tail  lift  outboard  of  the  fuselage  as  follows.  First, 

T 

the  constant  (-D^)  is  replaced  by  the  symbol  V2  which  is 

defined  by  an  expression  similar  to  equation  (D7)  but  for 
which  P^(«)  as  given  by  equation  (Cll)  is  defined  by  a 

summation  over  only  the  area  outboard  of  the  fuselage. 

If  the  panel  bounded  by  the  airplane  centerline  and  fuse¬ 
lage  side  (see  Figure  Fl)  is  designated  as  the  first 
panel  in  the  development  of  the  tail  generalized  forces, 
the  evaluation  of  V2  will  be  obtained  by  a  summation  start¬ 
ing  from  the  second  panel  and  ending  with  the  maximum 
tail  panel,  IT,  as  follows; 
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where  the  constant  P^Cy^*  00 )  is  given  by  equation  (C7) 
for  the  first  mode  (m=l) ,  in  which  case  f^  (y^)  -  1 

T  i 

and  f  (yr)  =  0.  Next,  the  generalized  force  variations 


M  G  T 

A-^(s),  A^(s),  and  ¥|(s),  which  actually  are  based  upon  a 

tail  area  extending  from  the  airplane  centerline  to  the 
tail  tip,  are  assumed  to  also  be  applicable  to  the  tail 
area  outboard  of  the  fuselage  side.  Finally,  the  non- 
TN  TC 

circulatory,  ^  (s),  and  circulatory,  (s),  functions 

due  to  tail  motion  and  stabilizer  rotation,  as  given  by 
equations  (E55)  to  (E68),  are  modified  to  include  only  the 
area  from  the  fuselage  side  to  the  tail  tip.  The  proce¬ 
dure  will  be  to  express  the  non- circulatory  contributions 

(s),  and  $£q(s)  components  of  the  circulatory  contri¬ 
butions  in  terms  of  the  normal  coordinates,  £n(s)»  and 

the  coefficients  which  will  now  he  based  upon  a  tail  area 
outboard  of  the  fuselage.  The  ^^(s)  components  will  be 

evaluated  as  before,  but  will  be  multiplied  by  the 
constant. 


V,  =  [ 


i=2 


Pi(yr» 


T 

)]/[  s 

i=l 


p^y?-00 


)] 


(F10) 


The  final  expression  for  the  tail  lift  based  on  an 
area  outboard  of  the  fuselage  is 


Lt(s)  =  Cy  JT(s) 


(Fll) 


where 


^T(s) 


N 


rpOO  T9  T 

=  2  JA^n(s)  +  B^(s)  +  C*„(s)] 


n  =1 


n  Jn 


n"n 


Too  t  °  T 

+  7  (s)  +  B*  7(s)  +  C*  7(s) 

2  3 

+  VX  2  *£kCs)  +  V2  Z_[A^k(s)  +  A^k(s)  +  ^(s)] 
k*l  k — 0 

(F12) 
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The  coefficients  A*,  b£,  c£.  Ay  By  and  C*  appearing  in 

equation  (F12>  are  obtained  from  equations  (E61)  to 
(E68)  for  m  »  1,  except  that  the  summations 
are  now  performed  only  over  the  outboard  (I_,-l)  panels. 
The  final  expressions  for  the  coefficients  are 


p  ^ 


i-2 


P^V^n^  +  *T(0)Bta<yi>] 


C  ,  T\ 


P*T(0) 


-  f  1 2  C(yi)b2(yi)a(yj) 


(F13) 


(F14) 


(F15) 


(F16) 


h 

B*  *  $  \  2  2C(y^)b(yJ)[  1  -  20T(O)[a(yJ)-^] 
c7  -  PbR  *T«,)1t2C(yi) 


(F17) 

(F18) 


where 


C(y])  -  Ay(y^)b(yJ)c^(yJ) 


(F19) 
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B.  Wing  Loads  Due  to  Wing  Motion  and  Gust  on  the  Wing 

The  mode  acceleration  method,  discussed  in  References  15 

and  17  ,  is  used  in  evaluating  the  wing  structural  loads 
due  to  wing  motion  and  gust  on  the  wing.  The  total  struc¬ 
tural  loads  are  due  to  the  inertial  contributions  (developed 
subsequently  in  Section  IIA)  and  the  aerodynamic  contribu¬ 
tions  caused  by  wing  motion  and  gust  on  the  wing  which  are 
developed  in  this  section.  The  assumption  is  made  that 
the  wing  airloads  associated  with  the  response  of  the  air¬ 
craft  to  the  gust  may  be  described  in  terms  of  the  rigid- 
wing  airload  distribution  required  to  produce  the  previously 
computed  rigid-body  vertical  acceleration. 

As  the  first  step  in  the  development  of  the  wing 
structural  load  equations  based  on  the  mode  acceleration 
method,  a  wing-induced  rigid-body  acceleration, 

is  defined  in  terms  of  the  airplane  rigid-body  accelera¬ 
tion,  £^(t).  For  this  purpose,  the  acceleration  due  to 

the  tail  lift  is  defined  by  dividing  the  tail  lift  by  the 
airplane  mass.  This  "tail- induced"  vertical  acceleration 
is  then  subtracted  from  the  airplane  acceleration  to  yield 
"wing- induced"  acceleration. 


..  ..  La  (t:) 

Hw^  =  "  I  -  Mj^/g-1  (F2°) 


The  tail  lift  outboard  of  the  center  line,  L.(t),  is  evalu¬ 
ated  from  equations  (F10)  and  (Fll)  as 


LA(t:)  -  VQLT(t)  where  VQ  «  1/V-j^ 


The  negative  sign  within  the  brackets  arises  from  the  fact 
that  lift  is  defined  as  positive  up  whereas  acceleration 
is  defined  as  positive  down.  Finally,  equation  (F20)  is 
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rearranged  as  follows: 


^1W^ 

ga 


^(t) 

S 


LA(t) 


M 


11 


g 

K 


(F21) 


where  the  acceleration  constant,  gQ,  is  introduced  to  non- 

a 

dimensionalize  each  term. 


The  aerodynamic  contributions  to  the  wing  structural 
loads  are  now  defined  in  terms  of  "aerodynamic  loading 
coefficients",  which  are  the  structural  loads  associated 
with  the  rigid-wing  airload  distribution  required  to  pro¬ 
duce  a  one  "g  "  upward  acceleration  of  32.2  feet  per  second 

square.  The  shear  (positive  up)  and  bending  moment., 
(positive  for  compression  in  the  top  skin)  at  the  ic 
load  point  may  be  expressed  in  terms  of  these  coefficients 

and  the  wing-induced  acceleration,  «w<*>  ,  as 


(t) 

M^A(t) 


-cf<V 


®a 

Hw(t) 

8a 


(F22) 

(F23) 


where  the  negative  signs  arise  from  the  opposite  definition 
of  the  aerodynamic  loading  coefficients  and  direction  of 
positive  acceleration.  The  torque  (positive  nose  up)  may 
also  be  expressed  in  a  form  similar  to  equations  (F22) 
and  (F23) .  It  is  assumed  that  the  aerodynamic  loading 
coefficients  for  the  torque  are  given  about  an  arbitrary 
axis  perpendicular  to  the  airplane  centerline,  and  that 
it  is  desired  to  calculate  the  torque  about  the  reference 
axis.  The  torque  about  a  perpendicular  axis  at  the 

i.L 

i  load  point  is  given  in  terms  of  the  known  aerodynamic 

If 

loading  coefficient,  caO^»XpA),  and  the  wing-induced 


l 


I 


} 
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acceleration  as 

T  (hyCt) 

*  -ca<Yi>XPA>  -T—  <F24> 


Note  that  equations  (F22)  to  (F24)  reduce  to  £he  aerody¬ 
namic  loading  coefficients  in  the  case  where  -  gfl. 

To  convert  equation  (F24)  to  the  torque  about  the  reference 
axis,  use  is  made  of  Figure  F2  in  which  the  planform  is 

divided  into  a  number  of  panels.  The  torque  at  the  iL 

load  point  about  the  perpendicular  axis  may  be  expressed 

in  terms  of  the  sum  of  the  panel  airloads,  fA(o,,t),  out- 

th  A  J 

board  of  the  i  load  point. 


(Y.,X 


i’APA 


^j“XPA)fA^j,t:) 


(F25) 


where  the  negative  sign  is  due  to  the  fact  that  an  up  load 
produces  a  nose-down  torque.  Hie  torque  about  a  point  on 
the  reference  axis  may  also  be  obtained  by  summing  up  the 
torques  outboard  of  the  ifc“  load  point. 


CA (^i»Xi> t) 


[  ^j'Xi^fA(oj,t:) 


^j-XPA^"^Xi_XPA^  ]fA^°j,t:^ 


(F26) 


Equations  (F25)  and  (F26)  are  combined  and  rearranged  as 
follows, 


tA^Yi,Xi,t^  ~  ^F27^ 


It  is  noted  that  the  shear  at  the  ic  load  point  is 
given  by 
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qWA 

bi 


(t) 


J 

2 

j-i 


fA^°  j  * 


(F28) 


Equations  (F22),  (F24) ,  (F27),  and  (F28)  are  now  combined. 


tAtti.Xi.t)  -  -[cJ(Y1.3SpA)+(Xi-XpA)cJ(Y1)]- 


go 


(F29) 


For  brevity  define  the  following  identity, 

TiA(t)  5  tA(Yi,Xi*t)  (F30) 

Therefore,  the  final  expression  for  the  torque  about  the 
reference  axis  becomes 

•  • 

T^A(t)  -  -[Ca(Yi»XPA)+(Xi"XPA)ca(Yi)]'J| -  (F31) 

3 


The  equations  for  the  shear,  bending  moment  and  torque 
at  the  elastic  axis  may  be  converted  to  the  non-dimensional 
time,  s,  by  use  of  equations  (Ell).  Hie  functions  L^(t), 

Sf(t),  M^(t),  and  T^A(t)  are  expressed  as  L^(s),  S^A(s), 

WA  WA 

(s)  and  T£  (s)  after  conversion  from  t  to  s.  Equa¬ 
tions  (F22),  (F23),  and  (F31)  become  after  rearranging  terms 


Sf(s) 

mWa(s) 


&_ 

Sa 

£_ 

Sa 

&- 

Sa 


Ef  00 

AiaWs> 

m  oo 

Aia^lW^S^ 


\ 


(F32) 


) 


g  m  ,  pounds (mass) 
ga  *  pounds (force) 


where 
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AIa  '  'CV  cf  <Yi> 

Aia  *  'CV  ca  <Yi> 

Aia  -  -VcI<Yi>)W+<Xi-XPA>cf<Yi>1 


(F33) 


J 


and  the  constant  C„  is  given  by  equation  (E19)  .  The  func- 
oo  V 

tion  ^y(s)  is  similar  in  form  to  equation  (F20)  and  is 
given  below  as 


OO 

^iw(s) 


OO 

^(s) 


.  WS> 

+  'sRir 


(F34) 


*  ^(s) 


V0iT(s) 


Tf 


11 


J 


where  the  first  expression  is  combined  with  equation  (Fll) 
to  obtain  the  second  expression. 
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II.  Inertial  Contributions  to  the  Wing  and  Tall  Structural 
Loads 


A.  Wing  and  Tail  Structural  Loads  Due  to  Motion 


Whether  the  mode-acceleration  or  force- summation 
method  is  used  in  evaluating  the  total  structural  loads, 
the  contribution  of  the  inertial  loads  is  formulated  in 
the  manner  discussed  below.  The  airfoil  surfaces  are  sub¬ 
divided  into  a  number  of  chordwise  strips  or  panels,  which 
is  not  necessarily  the  same  total  number  of  panels  as  that 
used  in  the  development  of  the  generalized  forces.  Hie 
notation  used  in  evaluating  the  inertial  contributions  to 
the  total  structural  loads  is  shown  in  Figure  F3.  The 


vertical  displacement  of  the  j*"*1  panel  center  of  gravity, 
z(yj,t),  positive  down  is  given  in  terms  of  the  vertical, 

h(y.,t),  and  angular,  a(y..,t),  deflections  at  the  reference 

axis  by  the  following  expression. 


z(yj,t)  -  h(y^. ,  t)+(Xj-Xj)a(yj ,  t) 


(F35) 


i-U 

The  acceleration  of  the  j  panel  center  of  gravity  is 
obtained  from  equation  (F35)  by  taking  the  derivative  with 
respect  to  real  time,  t. 


z(yj,t)  *  h(y^ ,  t)+(Xj ~Xj)cc(yj ,  t) 


(F36) 


,  th 


The  incremental  inertial  load,  fj(y^,t),  on  the  j  panel 

due  to  a  positive  acceleration  acting  down  is  equal  to  the 
product  of  the  panel  mass,  m(y.)/g,  and  the  acceleration, 

*  *  M  —  ■  J 

z(y.,t)  where  m(y.)  is  given  in  pound  mass  units. 

J  J 

fjCj^t)  *  -~^[h(yj, t) +(Xj-Xj )a(yj,t)]  (F37) 


If  an  unbalance  about  the  reference  axis  is  defined  as 


s(yj>Xj)  *  m(yj) (xj-Xj) 


(F38) 
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where  the  imbalance  is  positive  for  a  center  of  gravity 
location  aft  of  the  reference  axis,  equation  (F37)  becomes 

fjCyj.t)  -  |[ra(?j)h(yj,t)  +  sG?j,Xj)a(yj,t)]  (F39) 


The  vertical  and  angular  accelerations  may  be  expressed 
in  terms  of  the  contributions  of  the  normal  modes  and  con¬ 
verted  from  the  real  time  independent  variable,  t,  to  the 
non-dimensional  time,  s,  through  the  expressions 

h(y vt)  ■££.  (y,)S  (t)  (F40) 

J  n  n  J 

«<yvt)  -  s  f  (y*)i  (t)  (F41) 

J  n  n  J 

th 

and  equations  (Ell) .  The  j  panel  incremental  inertial 
load,  after  combining  equations  (F39)  through  (F41)  and 
equation  (Ell) ,  becomes 


(F42) 


where  Cy  is  given  by  equation  (E19) . 


th 


Hie  incremental  wing  shear  (positive  up)  at  the  i 
load  point  is  evaluated  by  summing  the  contributions, 

fj(yj,s),  of  all  the  panels  outboard  of  the  i  load  point. 


sf(S)  Cv  E  A™  en(s) 


wS 


(F43) 


where 


— 

rn  w 


>E  *  + 


(FA4) 
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and  the  letter  W  used  as  subscript  or  superscript  desig¬ 
nates  the  wing. 

The  incremental  wing  bending  moment  ^(positive  for 
compression  in  the  top  surface)  at  the  itH  load  point  is 
dependent  upon  the  summation  of  the  incremental  load, 
fjCy^jS),  times  the  moment  arm,  (y^-Y^),  outboard  of  the 

ifck  load  point. 


Hie  incremental  torque  (positive  for  nose-up)  about 
an  axis  perpendicular  to  the  airplane  centerline  which 
intersects  the  elastic  axis  at  the  1“*  load  point  is 
given  by 


Tj(s) 


J 

Z  [-(x.-X. 

j«i  J  1 


)fx(yj»s>  +  tjCyyXj.s)] 


(F47) 


where  tT(y.,x,,s),  the  contribution  of  the  panel  due 
J  J  »•  _ 

to  the  angular  acceleration,  a(yj,s),  is  given  in  terms 

of  the  center-of-gravity  moment  of  inertia  and  the  afore¬ 
mentioned  acceleration  as 

tl(^j’^j,S)  "  -CvHyy^)ZGys)  (F48) 


At  this  time,  the  center-of-gravity  moment  of  inertia  for 

f-U 

the  j  panel  is  expressed  in  terms  of  the  reference-axis 
moment  of  inertia,  i(yj,Xj)>  as  follows. 


135 


APPENDIX  F 


Kyj.xj)  -  1  <y j ■  x j ) '“(yj ) (xj "X j) 2 
Also,  equation  (F47)  is  rewritten  in  the  form. 


(F49) 


T  J  _ 

Ti(s)  *  *  t  ["^j“xi)fI(yj»s)-(xj-Xj)fI(yj»s) 


■+■  tj(yj)Xj,s) 


(F50) 


The  final  expression  for  the  inertial  contribution  to  the 
wing  incremental  torque  at  the  ifc“  load  point  is  obtained 
by  combining  equations  (F38),  (F40)  to  (F42),  and  (F48) 
to  (F50)  and  (Ell)  . 

t”X(s)  -CySjf  fn(s)  (F51) 


where 


.WT 

'in 


-  2 
j-i' 


[s(yj,Xj)£hn(yj>  +  KyJ.x^G")] 


+  s(5«.x«)£an(y«)Hx«-^)(F52) 


Note  that  in  the  case  where  the  reference  axis  is  perpen- 

TJ  JJ 

dicular  to  the  airplane  centerline,  that  is  xj"Xi  for 
i£j£Jw,  the  second  term  in  equation  (F52)  drops  out. 

For  the  tail  loads,  which  are  calculated  at  only  a 
single  point,  the  torque- tube  fuselage  juncture,  the  expressions 
are 
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5  ”n<8> 

(F53) 

■2  “n<8> 

(F54) 

2  ~n<8> 

(F55) 

The  coefficients  appearing  in  equations  (F53)  to  (F55)  are 
given  by  expressions  similar  to  equations  (F44),  (F46)  and 
(F52),  where  in  equation  (F52)  the  second  term  drops  out 
if  the  tail  reference  axis  (torque  tube)  is  perpendicular 
to  the  airplane  center  line  as  is  assumed  herein. 


2  [m(yT)fh  (yT) 
j=“2  J  n  J 


Z  [m(y?)f. 
j-2  1  hn 


+  s(yj»xj)fa  (yj)]  (F56) 

_ m  ni  ip 

+  s(yj >xpfajy j)  ]  (yj-Y{)  (F57) 


T 

where  represents  the  fuselage  station  of  the  inter¬ 
section  of  the  torque  tube. 

— 

AS  -  "^2[S(yj’Xj)£hn(yj)  +  i<yj'Xj)£an<yj)1  <M8) 

B.  Tail  Loads  Due  to  Stabilizer  Rotation 

The  expressions  for  the  inertial  contributions  to  the 
tail  structural  loads  due  to  stabilizer  rotation  will  be 
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similar  to  equations  (F53)  to  (F55)  with  the  variable 


oo  o  P 

£n(s)  being  replaced  by  y  (s) . 

8“<s>  -  Cv  A”  7(«) 

(F59) 

M^(s)  -  Cv  A™  *(s) 

(F60) 

T^(s)  -  Cv  A^  y(s) 

(F61) 

where  the  coefficients  in  the  above  equations  may  be  evalu¬ 
ated  from  equations  (F56)  to  (F58)  with  the  following 
values  for  the  mode  shapes, 


V7J>  ■ 0 
V#  ■ 1 


*•  rp  m  m  m 


(F62) 


(F63) 


(F64) 


(F65) 
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III.  Total  Wing  and  Tail  Loads 

The  total  structural  loads  at  the  1  wing  load  point 
are  given  by  ™"’ 


Si<3>  "  SiX<s>  +  SiA<s> 

Mi<s>  -  s)  +  M*JA(s) 

tJ(s)  *  T^Cs)  +  T^A(s) 


\ 

(F66) 


J 


where  the  component  parts  in  equations  (F66)  have  been 
developed  in  Sections  IIA  and  IB,  respectively. 


Similarly,  the  total  structural  loads  at  the  tail 
load  point  (defined  by  the  intersection  of  the  fuselage 
side  and  torque  tube)  are  given  by 


S*(s)  «  S*A(s)  +  sf  (s)  +  S^(s) 
M^(s)  =  H^A(s)  +  M^Cs)  +  M^Cs) 
T^(s)  -  T^A(s)  +  T^X(s)  +  T^Cs) 


(F67) 


y 


where  the  component  parts  in  equations  (F67)  are  developed 
in  Sections  IA,  IIA,  and  IIB,  respectively. 
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DEFINITION  OF  MOMENT  ARMS  FOR  AERODYNAMIC  CONTRIBUTIONS 
TO  TAIL  BENDING  MOMENT  AND  TORQUE 
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SYMBOLS 


All  symbols  used  in  this  volume  are  defined  below. 

The  symbols  for  the  more  important  or  fundamental  quan¬ 
tities  are  defined  verbally,  and  in  addition,  where  appro¬ 
priate,  reference  is  made  to  equations  and  supporting  fig¬ 
ures.  Other  symbols  which  are  of  lesser  importance  or 
which  represent  intermediate  quantities  with  less  physical 
significance,  are  defined  in  terms  of  their  analytical 
expressions  by  simply  listing  the  corresponding  equation 
or  figure  numbers. 


I.  Unit  Symbols 


Symbol 

Units 

F 

pounds 

Force 

I 

inches 

Linear  dimension 

L 

feet 

Linear  dimension 

P 

pounds 

Mass 

R 

radians 

Angular  dimension 

S 

semichords 

Non-dimensional  time 

T 

seconds 

Real  time 
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II.  Roman  Alphabet  Symbols 


Symbol  Units 


a  -  Chordwise  location  of  the 

reference  axis  of  a  two  dimen- 
sional  planform  measured  as  a 
fraction  of  the  semichord, 
positive  for  reference  axis 
aft  of  midchord. 


a 


Subscript  denoting  amplitude 
of  sinusoidal  function  (ex¬ 
cept  in  the  coefficients 


Aws  awm 

Aia>  Aia’ 


AWT 

Aia* 


the  constant 

ga,  and  the  parameters  in 
Appendix  A,  Section  II) . 


a(y.)  -  Chordwise  location  of  the 

1  reference  axis  at  the  span- 

wise  center  of  the  planform 

i*"*1  panel  measured  as  a 
fraction  of  the  planform 
local  semichord,  positive  for 
a  reference  axis  aft  of 
midchord,  (Figure  Cl) 

t~Vi 

a^  ---  Coefficient  of  k  exponential 

term  in  approximation  for  plan 
form  Kussner-type  function, 
(Equation  (B3)) 


a  .  -  Constant  associated  with  the 

analytical  approximation  of 
the  indicial  generalized  force 
functions,  (Equation  (C45)) 

9 

a„y  T“  Equation  (A24) 
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x  q-4 


A 

Subscript  or  superscript  de¬ 
noting  aerodynamic  contribution 
(except  in  the  variable  L^(s)). 

am 

Aia 

9  9 

fV/p 

Equation  (F33) 

A^ 

Aia 

9  9 

F  S  /PL 

Equation  (F33) 

at 

Aia 

9  9 

fV/P 

Equation  (F33) 

am 

Ain 

PL 

Equations  (F46)  and  (F57) 

A^ 

Ain 

P 

Equations  (F44)  and  (F56) 

at 

Ain 

PL 

Equations  (F52)  and  (F58) 

Al7 

PL2 

Equation  (F64) 

Al7 

PL 

Equation  (F63) 

A  TT 

Al7 

PL2 

Equation  (F65) 

Ak 

— 

Equation  (C41) 

A3n<*> 

L3/R 

Equation  (E15) 

Amn 

P 

Equations  (74),  (E24),  and 
(E61) 

Anry 

PL/R 

Equation  (75) 

at 

An 

P/R 

Equation  (F13) 

1 

ff 
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I 

I  -L<*> 

LA/R 

Equation  (E18) 

I 

J  A7 

PL/R 

Equation  (F16) 

] 

AR 

— 

Planform  aspect  ratio 

1 

b 

-  as 

L 

Semi chord  of  two  dimensional 
planform 

si 

b(yi> 

m 

L 

Semichord  at  the  spanwise  center 

of  the  planform  itb  panel, 

(Figure  Cl) . 

Pi 

bk 

w 

— 

Coefficient  of  ktb  exponential 
term  in  approximation  for  plan- 
form  Wagner- type  function, 

(Equation  (B14) ) . 

1  b  *. 

-  mk 

j 

— 

Constant  associated  with  the 
analytical  approximation  of  the 
indicial  generalized  force 
functions,  (Equation  (C45)) . 

J  \ 

L/S 

Airplane  reference  semichord, 

(Figure  B3) . 

bT 

1 

L/S 

Tail  reference  semichord,  half 
of  tail  mean  aerodynamic  chord, 
(Figures  Bl  and  B3) . 

1  bw 

1 

L/S 

Wing  reference  semichord,  half  of 
wing  mean  aerodynamic  chord. 

1 

1 
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N 

T 

Equation  (A25) 

bB 

L 

Reference  semichord  for  an  arbi¬ 
trary  planform,  (Appendix  B, 
Section  IA) . 

B 

— 

Scale  parameter  in  the  probability 
distribution  of  root-mean- square 
gust  velocity,  (Equation  (113)). 

\ 

-  -  - 

Equation  (C49) 

B<n<5,i) 

l3/rs 

Equation  (E15) 

B 

ran 

P/S 

Equations  (74),  (E25),  (E28), 

(E62),  and  (E63) . 

B 

nry 

PL/S 

Equations  (75),  (E66),  and  (E67) . 

Bn 

S/T 

Equation  (E92) 

Bn 

P/RS 

Equation  (F14) 

Btn(yi) 

l4/rs 

Equation  (E18) 

bt 

7 

PL/RS 

Equation  (F17) 

c(yi> 

— 

Chordwise  location  of  the  aero¬ 
dynamic  center  at  the  spanwise 

,  center  of  the  planform  i*"*1 

panel  measured  as  a  fraction  of 
the  planform  local  semichord, 
positive  for  aerodynamic  center 
aft  of  midchord,  (Figure  Cl) . 

I 

■ 
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I 

J  ca<V 

T 

FL 

fch 

Bending  moment  at  the  i  wing 
load  point  due  to  wing  airloads 
required  to  produce  incremental 
one  gfl  upward  normal  acceleration. 

J 

T 

positive  for  compression  in  the 
top  skin. 

t 

F 

Shear  at  the  l  load  point  due 

to  wing  airloads  required  to 
produce  incremental  one  gfl 

upward  normal  acceleration,  posi¬ 
tive  up. 

CI<V  XPA> 

FL 

Torque  at  the  l  load  point, 

about  the  perpendicular  axis 

XpA»  due  to  wing  airloads  re¬ 

*■« 

quired  to  produce  incremental 
one  gfl  upward  normal  acceleration. 

J 

positive  torque  tends  to  twist 
the  leading  edge  up. 

i 

J 

1/R 

Section  rigid  lift  curve  slope 
at  the  spanwise  center  of  the 
til 

planform  i  panel. 

J 

C 

I 

— 

Gust  power  spectrum  modifying 
constant,  (Equation  (125)). 

I 

— 

Superscript  designating  circula¬ 
tory  component. 

1  c(yi> 

l2/r 

Equations  (E8)  and  (E45) 

|  C(k) 

1 

“  “  “ 

Theodorsen's  function. 

1 

1 
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<4><v 

l3/rs2 

Equation  (E15) 

c 

mn 

P/S2 

Equations  (50),  (E29),  and  (E64) 

C 

nry 

PL/S2 

Equations  (50)  and  (E68) 

% 

1/T 

Equation  (E93) 

< 

P/RS2 

Equation  (F15) 

l4/RS2 

Equation  (E18) 

CV 

fs2/pl 

(V/bR)2/g 

CT 

7 

PL/RS2 

Equation  (F18) 

S 

LF/(R/T) 

Angular  damping  constant  of 
sprashpot  system,  (Figure  Al) . 

Dm 

PT/RS2 

Equations  (D7) ,  (D23) ,  and  (E84) 

^mn 

p/s2 

Equation  (77) 

f(s) 

Variable 

General  notation  for  a  dependent 
variable  (e.g.,  if  the  auto¬ 
pilot  were  being  considered, 
f(s)  would  be  replaced  by 

P-  (s))  . 

f*(s) 

Variable 

l-f(s) 

f  (p) 

Variable 

Laplace  transform  of  f(s). 

f(o>) 

or 

f(n) 

Variable 

General  notation  for  a  sinusoidal 
function. 

*.<“> 

Variable 

Amplitude  of  general  sinusoidal 
function. 

fh  <**> 

n 

—  —  - 

Relative  vertical  deflection  in 

the  nL  mode  at  the  aft  bob- 
weight  . 

\<*) 

Relative  vertical  deflection 
ttl 

in  the  n  mode  at  the  forward 
bobweight . 

£hm(yi) 

or  j 

\(yi> 

| 

Relative  vertical  deflection 

in  the  mt4x  or  n^  modes  of 
the  planform  reference  axis 
at  the  spanwise  center  of  the 

it*1  panel. 

- - 

Relative  vertical  displacement 

in  the  nt41  mode  at  the  it41 
airplane  location  (either 
fuselage,  wing,  or  tail). 

F 

Equation  (F39) 

fM(yi,s) 

L4/S2 

Equation  (E50) 

l4/s2 

Equation  (E14) 
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f  te<yI»  s>  l5/s2 


Equation  (E51) 


ft0(yi,s)  l5/s2 


Equation  (E17) 


Equations  (F25)  to  (F28) 


R 


Phase  angle  of  general  sin¬ 
usoidal  function. 


fc  <xa>  R/L 


Relative  pitching  rotation  in 

the  n  mode  at  the  aft  bob- 
weight  . 


fa  (xf )  R/L 


Relative  pitching  rotation  in 

the  n  mode  at  the  forward 
bobweight . 


£c  <*s>  R/L 


Relative  pitching  rotation  in 

the  n^  mode  at  the  stick  pivot 
point. 


R/L 


Relative  pitching  rotation  in 

the  n  mode  at  the  autopilot 
pitch  angle  sensing  device. 


Relative  pitching  rotation  in 
ttl 

the  n  mode  at  the  autopilot 
pitch  rate  sensing  device. 


fa<yi> 

m 

or 

\(yi) 


R/L 


Relative  pitching  rotation  in 

till 

the  m  or  n  modes  about  an 
axis  perpendicular  to  the  air¬ 
plane  center  line  at  the  inter¬ 
section  of  the  planform  refer¬ 
ence  axis  and  the  spanwise 
position  of  the  center  of  the 

ifck  panel. 


I 

1 

I 

I 


1 

1 
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1 

r 

f 

R/L 

Relative  pitching  rotation  in 

4*U 

the  n  mode  at  the  airplane 

location  (either  fuselage,  wing, 
or  tail) . 

I  *.<-> 

— 

Equation  (E78) 

J  f9<yi-t) 

L/T 

Equation  (E46) 

f0(t) 

~?t 

L/T 

Equations  (E3)  and  (E5) 

L/T 

Equation  (E9) 

w 

Vs) 

«  «  •» 

Equation  (D2) 

MK 

Fl-  F2 

1/S 

Equation  (77) 

F3 

1 

“  “  ~ 

ta; 

7  F(j)’  Fj 

— 

Equations  (C27)  and  (C47) 

! 

j 

J  F(S) 

— 

Equations  (C22)  and  (C42) 

I  F*(s) 

PL/S2 

Equations  (E27)  and  (E60) 

I  Fm(s) 

1 

L/T 

Equation  (E91) 

1 

I 
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F 

mn 

P/S3 

Equation  (77) 

g 

PL/FT2 

Gravitational  constant  which 
includes  conversion  from  mass 
to  force  units. 

sa 

L/T2 

Acceleration  constant  32.2  feet/ 

second  squared. 

* 

G 

-  — 

Subscript  or  superscript  de¬ 
noting  gust  or  downwash  due 
to  gust. 

»  T 

'1 

h(t) 

L 

Vertical  displacement  of 
reference  axis  of  two-dimensional 
airfoil  section,  positive  down. 

i 

l 

1 

h  (xfl,t) 

L/T2 

Vertical  acceleration  at  the  aft 
bobweight,  positive  down. 

I 

i 

i 

h  (xf,t) 

L/T2 

Vertical  acceleration  of  the 
forward  bobweight,  positive  down. 

j 

h(yt,t) 

L 

Vertical  displacement  of  reference 
axis  at  the  spanwise  center  of 

the  ifck  panel,  positive  down. 

yi 

m 

\ 

wW 

hGFj.t) 

L 

Vertical  displacement  of  reference 
axis  at  the  spanwise  center-of- 
gravity  position  of  the 

j  panel,  positive  down. 

I 

I 

I 

I 

I 
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^(s) 

• 

i 

L 

Vertical  displacement,  velocity 
and  acceleration  at  the 

hi<8) 

{^(s) 

L/T 

L/T2 

i  airplane  location,  posi¬ 
tive  down. 

*V^mac* 

L 

Vertical  displacement  of  the 
wing  reference  axis  at  the 
spanwise  position  of  the 
tail  mean  aerodynamic  chord, 
positive  down. 

H20 

1/S 

Equation  (A27) 

H21 

1/S2 

Equation  (A28) 

H22 

1/S2 

Equation  (A29) 

H30 

1/s2 

Equation  (A19) 

H3n 

1/L 

Equation  (A20) 

H3. 

1/S 

Equation  (A19) 

H40 

1/S 

Equation  (A21) 

1 

— 

General  index. 

i 

_  -  — 
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Moment  of  inertia  of  planform 

j  panel  about  an  axis  perpen¬ 
dicular  to  the  airplane  center- 
line  passing  through  the  panel 
center  of  gravity. 


Moment  of  inertia  of  planform 

panel  about  an  axis  perpen¬ 
dicular  to  the  airplane  center- 
line  passing  through  a  point 
defined  by  the  intersection  of 
the  spanwise  location  of  the 

panel  center  of  gravity  and 
the  reference  axis,  (Figure  F3) . 


Subscript  or  superscript  de¬ 
noting  inertial  contribution. 


Number  of  planform  panels  in 
generalized  force  calculations. 


Aft  bobweight  inertia  about 
its  pivot  point. 


Forward  bobweight  inertia 
about  its  pivot  point. 


Equation  (79) 


Equation  (79) 


Control  system  inertia  referred 
to  stick  pivot  point,  (Equation 
(A10)). 


1 

1 
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1 

isO 

1 

PL2 

Inertia  of  the  control  system 
about  the  stick  pivot  point 
without  the  bobweight  contri¬ 
bution,  (Equation  (A10) ) . 

I 

V 

PL2 

Inertia  of  the  bobweight  system 
about  the  stick  pivot  point, 

(Equation  (A9)) . 

1 

i 

Variable 

Integral  defined  by  Equation 
(114)  . 

i  1 

— 

General  index. 

i 

— 

Number  of  planform  panels  in 
load  calculations . 

1  Jl*  J2 

1/S 

Equation  (79) 

I 

i 

J  J“ 

P/S3 

Equation  (79) 

J 

PL  As  3 

Equation  (79) 

I 

— 

General  index. 

I 

k7 

R/  (R/T) 

Pitch  rate  autopilot  gain. 

1 

I 

rA 

Pitch  displacement  autopilot  gain. 

1 

1 
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LF/R 

kt  LF/R 


Stabilizer  to  control  stick 
linkage  ratio. 


Angular  spring  constant  of 
sprashpot,  (Figure  Al) . 


Angular  spring  constant  of 
stick  feel  spring,  (Figure  Al) . 


Number  of  exponential 
terms  in  approximation  to 
indicial  lift  function  for 
planform  motion,  (Equation  (B14)). 


Number  of  exponential 
terms  in  approximation  to 
indicial  lift  function  for 
gust  on  a  planform,  (Equation 
(B3)) . 


3 

] 

1 

1 

I 

I 

I 

I 

r 


Vs) 

PL/RS2 

Equation  (F12) 

Vy£>s) 

F 

Equation  (E76) 

Vyi»s) 

F 

Equation  (E48) 

VO 

F/L 

Equation  (El) 

Vyi»s) 

F 

Equation  (E13) 

Vyi»S> 

F 

Equation  (Dl) 

I 

I 

I 

I 

I 

I 

S 

I 
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L 

L 

Scale  of  turbulence  parameter. 

L10 

1/S 

Equation  (A6) 

Lln 

S/L 

Equation  (A6) 

c 

CN 

hJ 

1/L 

Equation  (A6) 

L3n 

1/LS 

Equation  (A6) 

CO 

F 

Tail  lift  per  side  based  on 
area  from  airplane  center  line 
to  tail  tip,  positive  up, 
(Equation  (F8)). 

lT(s) 

F 

Tail  lift  per  side  based  on  area 
from  fuselage  juncture  to  tail 
tip,  positive  up,  (Equation 
(Fll)). 

m 

— 

Normal  mode  index. 

-(yj) 

P 

Mass  of  planform  j panel. 

ma 

P 

Mass  of  aft  bobweight. 

mac 

— 

Mean  aerodynamic  chord. 

mf 

P 

Mass  of  forward  bobweight. 

M 

-  -  - 

Airplane  Mach  number. 
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Subscript  or  superscript  de¬ 
noting  motion  or  downwash  due 
to  motion. 


H 


Superscript  denoting  bending 
moment. 


M 


M?(s) 


FL 


Gust  front  Mach  number. 


Total  incremental  bending  moment 
ttl 

at  the  i  wing  load  point  posi¬ 
tive  for  compression  in  the  top 
skin,  (Equation  (F66)). 


Mmm 


Normal-mode  generalized  masses 
based  on  half  the  airplane. 


M 


nry 


PL 


Normal-mode  stabilizer-rotation 
generalized  mass  coupling  terms 
based  on  half  the  airplane. 


M 


77 


PL 


Stabilizer-rotation  on-diagonal 
generalized  mass  term  based  on 
half  the  airplane  (stabilizer 
pitching  inertia  about  reference 
axis) . 


MfX(s) 


FL 


Inertial  contribution  to  M^(s) 

due  to  stabilizer  rotation, 
(Equation  (F60)). 


M^(s) 


FL 


Total  tail  incremental  bending 
moment  at  the  intersection  of 
the  fuselage  and  the  torque 
tube,  positive  for  compression 
in  the  top  skin,  (Equation 
(F67». 
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MlT(t)  FL 

M2T(t)  FL 

n,  n 

N 

N 


Equation  (All) 

Equation  (A13) 

Normal  mode  index  (1  £  n  <  N, 
l£n£N  +  l='7* 

Number  of  symmetric  normal  modes. 

Superscript  designating  non- 
circulatory  component. 


N„  1/T 

y 


Characteristic  frequency,  average 
number  of  times  per  second  that 
the  variable  y  crosses  the  value 
zero  with  a  positive  slope, 
(Equation  (111)). 


N(y)  1/T  Average  number  of  times  per 

second  that  the  variable  crosses 
a  given  value  of  y  with  a  posi¬ 
tive  slope;  an  approximation  for 
the  average  number  of  positive 
peaks  per  second  exceeding  a 
given  value  of  y,  (Equation  (113)). 


Ratio  of  aft  bobweight  angular 
motion  to  control  stick  angular 
motion. 


Ratio  of  forward  bobweight 
angular  motion  to  control  stick 
angular  motion. 


1/S  I  i-l-5,  8 

1/S2,  i-6 


Equation  (72) 
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0 


i 


1/S,  i=l-5,  8  Equation  (67) 
1/S2,  i-6 


P 


Laplace  transform  differential 
operator . 


pm(^i’  °°>  l2/R 


Equation  (C7) 


?m(oo)  L2/R 


Equation  (Cll) 


P 


1 


Percentage  of  total  flight  time 
or  distance  in  turbulence. 


ttl 

q (y. ,s)  F  Generalized  force  on  the  ic 

’rm'-'x 

panel  in  the  in  mode  due  to 

r  i.1- 

the  r  contribution,  (Equa¬ 
tions  (C5) ,  (C6) ,  (D3) ,  and 
(E20)) . 


1/S*  i*l-5,  8 
1/S2,  i-6 


Equation  (62) 


Vs> 

<Ws> 


Total  generalized  force  in  the 


mfc^  mode,  (Equations  (36)  and 
(57)). 


til  •*— 

Generalized  force  in  the  m 

mode  due  to  the  r^  contribu-  T 

tion,  (Equations  (C8)  and  (D5)).  «L 

I 

I 

I 

I 


1 

1 
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I 

I 

I 

J 

J 

«•  mm  m 

Index  denoting  generalized 
force  contributions : 
r=l,  wing  motion 

2,  tail  motion  and  stabilizer 
rotation 

3,  downwash  on  tail  due  to 
wing  motion 

4,  downwash  on  tail  due  to 
gust  on  the  wing 

5,  gust  on  the  tail 

6,  gust  on  the  wing 

J 

L 

Moment  arm  of  the  aft  bobweight 
from  its  pivot  point,  (Figure  Al) . 

J 

T 

L 

Moment  arm  of  the  forward  bob- 
weight  from  its  pivot  point, 

(Figure  Al)  . 

J 

1 

— ,  i=7 

1/S,  1=1-5,  8 

1/S2,  1=6 

Equation  (56) 

I 

I 

S 

Non-dimensional  time  based  on 
airplane  reference  semichord, 

(Equation  (B5)). 

I 

J 

S 

Transport  time  delay  between 
airplane  zero  time  reference 
point  and  assumed  start  of 
stabilizer  lift  build  up  due  to 
downwash  resulting  from  gust  on 
the  wing. 

I  “m 

1 

1 

S 

Transport  time  delay  between 
wing  motion  and  assumed  start 
of  resulting  stabilizer  lift 
build  up  due  to  downwash. 

1 

1 
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Transport  time  delay  between 
airplane  zero  time  reference 
point  and  start  of  tail  lift 
build  up  due  to  gust  on  the 
tail. 


Non-dimensional  horizontal  dis¬ 
tance  from  the  midchord  of  the 
wing  mean  aerodynamic  chord  to 
the  leading  edge  of  the  tail 
mean  aerodynamic  chord,  (Figures 
B1  and  B6) . 


Non-dimensional  vertical  dis¬ 
tance  between  the  planes  con¬ 
taining  the  wing  and  the  tail, 
(Figure  Bl) . 


General  expression  for  transport 
time  delay. 


Transport  time  delay  between 
airplane  zero  time  reference 
point  and  assumed  start  of  wing 
lift  build  up  due  to  gust  on 
the  wing. 


i.L 

Unbalance  of  planform  panel 
about  an  axis  perpendicular  to 
the  airplane  centerline  passing 
through  a  point  defined  by  the 
intersection  of  the  spanwise 

4.1. 

location  of  the  j  panel  center 
of  gravity  and  the  reference 
axis,  positive  for  center  of 
gravity  aft  of  reference  axis, 
(Equation  (F38)). 


Superscript  denoting  shear. 


Superscript  denoting  stabilizer 
rotation. 
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S“(s) 

F 

Total  incremental  shear  at  the 

i  wing  load  point,  positive 
up,  (Equation  (F66)). 

sfI(s) 

F 

T 

Inertial  contribution  to  S-^(s) 

due  to  stabilizer  rotation, 
(Equation  (F59)). 

S^(s) 

F 

Total  tail  incremental  shear 
at  the  inersection  of  the  fuse¬ 
lage  and  the  torque  tube,  posi¬ 
tive  up,  (Equation  (F67)). 

t 

T 

Real  time. 

fcA  (Yi » Xi  > 

FL 

Equation  (F26) 

fcA  ^Yi  ’  XPA  * 

1  FL 

Equation  (F25) 

tlCYj.x.,8) 

FL 

Equation  (F48) 

FL 

Equation  (E79) 

t9(yj,t) 

FL 

Equation  (E44) 

V<o 

F 

Equation  (E2) 

Vyi,l:) 

F 

Equation  (E7) 

FL 

Equation  (C4) 

T 

— 

Subscript  or  superscript  de¬ 
noting  tail. 

GRUMMAN  AIRCRAFT  ENGINEERING  CORPORATION 


REPORT 

DATE 


ADR  06-14-63.1 
OCTOBER  1963 


164 


T 


Superscript  denoting  torque. 


T(s)  FL  Kinetic  energy  in  terms  of 

normal  velocity  coordinates 
and  generalized  masses,  (Equa¬ 
tion  (34))  . 


tV(s)  FL  Total  incremental  torque  about 

a  wing  axis  perpendicular  to 
the  airplane  centerline  passing 
through  a  point  defined  by  the 
intersection  of  the  dpanwise 

location  of  the  ic  load  point 
and  the  reference  axis,  posi¬ 
tive  nose  up,  (Equation  (F66)) . 


(s) 


T 

FL 


Equations  (A22)  to  (A25) 

Inertial  contributions  to 
T 

T^(s)  due  to  stabilizer  rota 
tion,  (Equation  (F61)). 


T 

T-^(s)  FL  Total  tail  incremental  torque 

about  the  torque  tube  at  the 
intersection  of  the  fuselage 
and  torque  tube,  positive  nose 
up,  (Equation  (F67)) . 


T  (<x>)  Variable  Frequency-response  function  of 

y  arbitrary  variable,  (Equation 

(HO)). 


V 


L/T 


Airplane  forward  velocity. 


V 


0 
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V1 

“  -  “ 

Equation  (FlO) 

V2 

PT/RS2 

Equation  (F9) 

V(s) 

FL 

Potential  energy  in  terms  of 
normal  coordinates  and  general¬ 
ized  spring  terms,  (Equation 
(35)). 

WG 

L/T 

Vertical  gust  velocity. 

wO^.M.A^, 

ST*  ST* 

— 

Indicial  downwash  function, 
normalized  unsteady  downwash 
build  up  on  the  tail  due  to 
either  a  step  change  in  wing 
motion  or  a  step  change  in 
gust  velocity  on  the  wing, 
(Equation  (26)). 

w6(s) 

L/T 

Equation  (E75) 

we (yj.t) 

L/T 

Equation  (E47) 

Vc) 

L/T 

Downwash  at  the  three-quarter 
chord  of  a  two-dimensional 
planform,  (Equation  (E4)) . 

44 

a 

*t4 

•©- 

£ 

'L/T 

Downwash  at  the  three-quarter 
chord  at  the  spanwise  center 

of  the  i^  panel  of  a  three- 
dimensional  planform,  (Equation 
(E10)). 
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L/T 


L  or  I 

L 


L  or  I 


L  or  I 


L  or  I 


L 


L  or  I 


Equation  (E71) 


Subscript  or  superscript  de¬ 
noting  wing. 


Fuselage  station,  positive 
a£t  of  nose. 


Tail  torque  arm, longitudinal 
distance  between  tail  aero¬ 
dynamic  center  and  tail  refer¬ 
ence  axis,  positive  for  aero¬ 
dynamic  center  forward  of  refer 
ence  axis,  (Figure  FI) . 


Location  of  the  aft  bobweight 
pivot  point,  (Figure  Al) . 


Location  of  airplane  center  of 
gravity,  (Figure  Al) . 


Location  of  the  forward  bob- 
weight  pivot  point,  (Figure  Al) 


f*h 

Longitudinal  location  of  j 
panel  center  of  gravity, 
(Figure  F3)  . 


Location  of  stick  pivot  point, 
(Figure  Al)  . 


Fuselage  station  of  stabilizer 
torque  tube. 


L  or  I 


1.67 


,a 


U 


L  or  I 


Location  of  the  autopilot 
pitch  angle  sensing  device. 


L  or  I 


Location  of  the  autopilot 
pitch  rate  sensing  device. 


X,  L  Longitudinal  location  of  the 

reference  axis  at  the 

i^  load  point,  (Figures 
F2  and  F3) . 


Xp.  L  Longitudinal  location  of  an 

axis  perpendicular  to  the 
fuselage  centerline  about  which 
the  torque  components  of  the 
data  for  the  aerodynamic  load¬ 
ing  coefficients  used  in  the 
mode-acceleration  method  are 
provided,  (Figure  F2) . 


y  L  or  I  Spanwise  coordinate  measured 

from  the  airplane  centerline 
perpendicular  to  the  plane  of 
symmetry,  positive  outboard, 
(Figure  Cl)  . 


?i 


L  or  I  Spanwise  location  of  the  center 

of  the  ith  panel,  (Figure  Cl) . 

th 

L  or  I  Spanwise  location  of  the  j 

panel  center  of  gravity, 

(Figure  F3) . 


L  or  I  Spanwise  location  of  the  plan- 

form  mean  aerodynamic  chord, 
(Figure  B3) . 
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^max 

L  or  I 

Planform  semi span. 

yo 

L 

Tail  bending  moment  arm, span- 
wise  distance  between  tail 
aerodynamic  center  and  nominal 
fuselage  side,  positive  for 
aerodynamic  center  outboard 
of  fuselage,  (Figure  FI) . 

Y. 

l 

L  or  I 

Spanwise  location  of  the  ifc^  load 
point,  (Figures  F2  and  F3) . 

z(7j,  t) 

L 

Vertical  displacement  of  the 

j  ^  panel  center  of  gravity, 
(Equation  (F35)). 

! 


169 


III.  Greek  Alphabet  Symbols 


Symbol 

Units 

a 

— 

Subscript  denoting  phase  angle 
of  sinusoidal  function. 

a(t) 

R 

Angular  displacement  about  the 
reference  axis  of  two-dimen¬ 
sional  airfoil  section,  positive 
nose  up. 

a(xa,t) 

CM 

Pt 

Angular  acceleration  at  the  aft 
bobweight  pivot  point,  positive 
nose  up. 

a  (xf>t) 

R/T2 

Angular  acceleration  at  the 
forward  bobweight  pivot  point, 
positive  nose  up. 

«(xg,t) 

R/T2 

Angular  acceleration  at  the 
stick  pivot  point,  positive 
nose  up. 

R 

Angular  displacement  at  the 
location  of  the  autopilot  pitch 
attack  sensing  device,  positive 
nose  up. 

• 

*  ,  a  - 
a(x^,t) 

R/T 

Angular  velocity  at  the  location 
of  the  autopilot  pitch  rate 
sensing  device,  positive  nose  up 
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Angular  displacement  about  an 
axis  perpendicular  to  the  air¬ 
plane  centerline  passing  through 
a  point  defined  by  the  inter¬ 
section  of  the  planform  refer¬ 
ence  axis  and  the  s panwise 
position  of  the  center  of  the 
th 

ic  panel,  positive  nose  up. 


Angular  displacement  about  an 
axis  perpendicular  to  the 
airplane  centerline  passing 
through  a  point  defined  by  the 
intersection  of  the  planform 
reference  axis  and  the  spanwise 
center-of-gravity  position  of 

the  j**1  panel,  positive  nose  up. 


Angular  displacement,  velocity 
and  acceleration  about  an 
axis  perpendicular  to  the 
airplane  centerline  at  the 

i  airplane  location,  posi¬ 
tive  nose  up. 


Exponent  of  the  k*"  exponential 
term  in  approximation  for  plan- 
form  Kiissner-type  function 
based  upon  the  non-dimensional 
time,  s,  (Equations  (B9)  and 
(B12) ) . 


th 

Exponent  of  the  k  exponential 
term  in  approximation  for 
planform  Kiissner-type  function 
based  upon  the  non-dimensional 
time,r| ,  (Equations  (B9)  and 
(B12)) . 


1 

1 
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1  . 

I  amk 

I 

1/S 

Constant  associated  with  the 
analytical  approximation  of 
the  indicial  generalized  force 
functions,  (Equations  (C45) 
and  (C52)). 

I 

Ak 

1/S 

Equation  (C37) 

T  i 

°V^tnac* 

J 

I 

f 

s)  R 

Angular  displacement  about  an 
axis  perpendicular  to  the  air¬ 
plane  centerline  passing  through 
a  point  defined  by  the  inter¬ 
section  of  the  wing  reference 
axis  and  the  spanwise  position 
of  the  tail  mean  aerodynamic 
chord,  positive  nose  up. 

J 

pk 

I 

I 

1/S 

Exponent  of  the  exponential 

term  in  approximation  for  plan- 
form  Wagner-type  function  based 
upon  the  non-dimensional  time, 
s,  (Equations  (B16)  and  (B18)). 

I  * 

I 

1/S 

Exponent  of  k*"*1  exponential 
term  in  approximation  for  plan- 
form  Wagner-type  function  based 
upon  the  non-dimensional  time. 

T),  (Equations  (B16)  and  (B18)). 

I 

I 

— 

Subscript  denoting  stabilizer 
rotation. 

J§ 

|  7(s) 

1 

R 

Stabilizer  rotation,  positive 
for  trailing  edge  down. 

1 

1 
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7S (t)  R 


Total  signal  input  to 
stabilizer. 


6(s) 

6(Mg,M,ARw,  , 
ART,  Sljl,  s) 


Tail  indicial  lift  function 
due  to  downwash  caused  by 
wing  motion  (M  =  “)  or 

gust  on  the  wing  (Mg  =0  ), 

(Equations  (26),  (B19),  (B20), 
and  (B21)). 


(s) 


Indicial  generalized  force 

th 

function  for  the  m  mode 
due  to  downwash  on  the  tail 
caused  by  gust  on  the  wing. 
(Equations  (C18)  and  (C19)) 


6M 

m 


(s) 


Indicial  generalized  force 

function  for  the  m  mode 
due  to  downwash  on  the  tail 
caused  by  wing  motion, 
(Equations  (C16)  and  (C17)). 


r 

Am(s)  L/T  Generalized  force  variation 

1  due  to  downwash  on  the  tail 

caused  by  gust  on  the  wing, 
(Equation  (D24)). 


A^(s)  L/T 


Q 

Component  parts  of  Am(s), 
(Equations  (D25)  to  (D27)). 


A«(s)  L/T 


Generalized  force  variation 
due  to  downwash  on  the  tail 
caused  by  wing  motion,  (Equa¬ 
tion  (E8/)). 
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ASk<*> 


L/T 


M 

Component  parts  of  A^(s), 
(Equations  (E88)  to  (E90) ) . 


As 


Non-dimensional  time  increment. 


Ay(yt) 


Width  of  the  i  panel  whose 
center  is  located  at  y^, 

(Figure  Cl). 


de/da 


Steady-state  downwash  slope. 


Exponent  defining  rate  of  decay 
of  frequency-response  function 
at  high  frequencies,  (Equations 
(119)  to  (121)). 


Coefficients  in  equation  (C30) . 


Longitudinal  location  of  the 
aerodynamic  center  of  the 

j*"*1  panel,  (Figure  F2) . 


Non-dimensional  time  based  on 
planform  reference  semichord, 
(Equation  (Bl)). 


T)(t) 


Motion  of  spring-damper  attach¬ 
ment  point  in  sprashpot,  positive 
for  positive  stick  rotation, 

(Figure  Al) . 


GRUMMAN  AIRCRAFT  ENGINEERING  CORPORATION 


REPORT 

DATE 


ADR  06-14-63.1 
OCTOBER  1963 


174 


Subscript  denoting  Imaginary 
component. 


K 


R/S 


Reduced  frequency, 
(3)). 


(Equation 


\ 


Equations  (C28)  and  (C34) 


Equation  (C24) 


^  Sweep  angle  of  planform  leading 

edge,  (Figure  Cl) . 


^(s)  R  Autopilot  signal  to  stabilizer, 

positive  for  nose  up  pitch 
attitude  and  pitch  rate, 
(Equation  (A5)). 


v 


i 

k 


Equation  (C50) 


£  (s) 
mv  7 

or  L  Normal  coordinates. 

e„(s) 


2 

^iy(s)  L/S  Wing-induced  rigid-body  verti¬ 

cal  acceleration,  (Equation 
(F34)  ) . 


p 


Subscript  denoting  real 
component . 


P 


Air  density. 


6 


Dummy  variable  of  integration 
in  Duhamel ' s  integral . 


<3(y.)  S  Transport  time  delay  in  terms 

x  of  airplane  reference  semi¬ 

chords  to  the  leading  edge  of 

the  midspan  of  the  i*"*1  panel 
measured  from  the  planform 
apex,  (Figure  Cl) . 


6 (y^)  S  Transport  time  delay  in  terms 

of  planform  reference  semichords 
to  the  leading  edge  of  the  mid- 

span  of  the  i  panel  measured 
from  the  planform  apex, 

(Equation  (C2)). 


6max  S  Transport  time  delay  in  terms 

of  airplane  reference  semi¬ 
chords  to  the  tip  chord  lead¬ 
ing  edge  measured  from  the 
planform  apex. 


L/T 


Root-mean- square  gust  velocity. 


Variable  Root-mean- square  value  of 

variable  y,  (Equation  (112)). 
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t(s) 


’j 

0(s) 


0(M,AR,t!) 


0(ti) 


\(.S) 


R  Stick  rotation,  positive  for 

stick  movement  towards  the 
pilot,  (Equation  (A17)). 


L  or  I  Spanwise  location  of  the  center 

of  the  j  panel,  (Figure  F2) . 


Exponential  approximation  to 
the  planform  Wagner-type 
function,  including  compress¬ 
ibility  and  aspect  ratio  effects 
and  based  upon  the  airplane 
reference  semichord,  (Equations 
(B15)  and  (B17)). 


Planform  indicial  lift  function 
due  to  motion,  including  com¬ 
pressibility  and  aspect  i?atio 
effects,  and  based  upon  the  plan- 
form  reference  semichord;  Wagner f 
type  function,  (Equation  (B13) ) . 


Exponential  approximation  to 
the  planform  Wagner-type 
function,  including  compress¬ 
ibility  and  aspect  ratio 
effects,  and  based  upon  the 
planform  reference  semichord, 
(Equation  (B14)). 


PL/S2 


Unsteady  contribution  to  the 
circulatory  generalized  force 
variation  due  to  motion, 
(Equations  (E39)  and  (E57)). 
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1 

I  *£<*> 

I 

PL/S2 

Circulatory  contribution  to 
the  total  generalized  force 
variation  due  to  motion, 

(Equations  (E26),  (E38),  and 
(E56)) . 

I  *•> 

I 

PL/S2 

Non-circulatory  contribution 
to  the  total  generalized  force 
variation  due  to  motion, 

(Equations  (E23)  and  (E55)). 

I 

J  *nk<s> 

J 

PL/S2 

Component  parts  of  <I>m(s) , 
(Equations  (E41)  and  (E59)). 

J  V* 

J 

l2/rt 

Power-spectral-density  function 
of  vertical  gust  velocity  in 
terms  of  circular  frequency, 
(Equations  (107)  and  (108)). 

J  <D  (co) 

G 

J 

f 

l2At 

Modified  power-spectral-density 
function  of  vertical  gust 
velocity  in  terms  of  circular 
frequency,  (Equation  (125)). 

J 

_  (n) 

I 

1 

l3At2 

Power-spectral-density  function 
of  vertical  gust  velocity  in 
terms  of  reduced  frequency, 
(Equation  (106)). 

■ 

|  *y<“> 

1 

Variable 

Power-spectral-density  function 
of  response  variable  in  terms 
of  circular  frequency, 

(Equation  (110)). 

1 

1 
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4>  (ft)  Variable  Power- spectral-density  function 

y  of  response  variable  in  terms 

of  reduced  frequency, 

(Equation  (109)). 


X(s)  L/T 


Vertical  gust  velocity. 


xj 


L 


Longitudinal  location  of 
reference  axis  at  a  spanwise 
position  corresponding  to  the 
center  of  gravity  of  the 

panel,  (Figure  F3)  . 


^(M,AR,q)  -  Planform  indicial  lift  function 

due  to  gust,  including  com¬ 
pressibility  and  aspect  ratio 
effects,  and  based  upon  the 
planform  reference  semichord; 
Kiissner-type  function, 

(Equation  (B2)). 


^(s)  ---  Exponential  approximation  to 

the  planform  Kiissner-type 
function,  including  compress¬ 
ibility  and  aspect  ratib  effects 
and  based  upon  the  airplane 
reference  semichord,  (Equations 
(B7)  and  (BIO)). 


^0l)  -  Exponential  approximation  to 

the  planform  Kiissner-type 
function,  including  compress¬ 
ibility  and  aspect  ratio 
effects,  and  based  upon  the 
planform  reference  semichord, 
(Equation  (B3) ) . 
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I 

("  I 

J 

J  fm(s)  L/T 

J 


Indicial  generalized  force 

A.1. 

function  for  the  nr  mode 
due  to  gust,  (Equations  (C12), 
(C13),  (C14) ,  and  (C15)). 


Generalized  force  variation 
due  to  gust,  (Equations  (D8) 
and  (D17)). 


Ws>  L/T 


Component  parts  of 
Tm(s),  (Equations  (D18)  to 

(D20)). 


co 


R/T 


Circular  frequency. 


o^(L)  L/T 


Breakpoint  frequency  in  approxi¬ 
mation  of  power- spectral-density 
function  of  vertical  gust  vel¬ 
ocity,  (Equation  (117)) . 


J 

J 

I 

I 

I 

I 

I 


CD 


m 


R/T 


Natural  frequency  of  the  m 
mode. 


cd  ,  1/S  Constant  associated  with  the 

analytical  approximation  of 
the  indicial  generalized  force 
functions,  (Equations  (C45) 
and  (C52)) . 


rA 


Reduced  frequency,  (Equation 
(105)). 
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IV.  Mathematical  Symbols 


Symbols 


(•), 

(“) 

Dots  above  symbols  denote  differentiation 
with  respect  to  real  time,  t. 

(#), 

(°°) 

Circles  above  symbols  denote  differentiation 
with  respect  to  non-dimensional  time,  s, 
(Equations  (9),  (10),  and  (Ell)). 

( )' 

Prime  after  symbol  denotes  differentiation 
with  respect  to  independent  variable. 

1  1 

Absolute  value. 

a 

3a 

Partial  derivative  with  respect  to  a. 

£[ 

] 

Laplace  transform  of  [  ] . 
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